THE・つろく

2017年4月26日水曜日

Generating PIOMAS_initialization.nc

 1014  vi namelist.piomas2naa 
 1015  ../bin/sosie.x -f namelist.piomas2naa 
 1016  vi namelist.piomas2naa 
 1017  ../bin/sosie.x -f namelist.piomas2naa 
 1019  ncks -h -A SIC_PIOMAS_1979_2013-naa_arlan.nc4 SIT_PIOMAS_1979_2013-naa_arlan.nc4 
 1021  mv SIT_PIOMAS_1979_2013-naa_arlan.nc4 PIOMAS_NAA_1979_2013.nc4
 1022  mv PIOMAS_NAA_1979_2013.nc4 ~/scratch/input_naa/hhayashi/piomas/
 1024  cdo seltimestep,1/1 ~/scratch/input_naa/hhayashi/piomas/PIOMAS_NAA_1979_2013.nc4 PIOMAS_initialization.nc
 1031  cdo chname,SIT,hicif PIOMAS_initialization.nc PIOMAS_initialization.nc2 
 1033  cdo chname,SIC,frld PIOMAS_initialization.nc2 PIOMAS_initialization.nc3

 1041  ncap -s "frld2=1.-frld" PIOMAS_initialization.nc3 PIOMAS_initialization.nc4
 1043  ncks -x -v frld PIOMAS_initialization.nc4 PIOMAS_initialization.nc5
 1046  cdo chname,frld2,frld PIOMAS_initialization.nc5 PIOMAS_initialization.nc6
 1050  mv PIOMAS_initialization.nc6 ~/scratch/input_naa/hhayashi/piomas/PIOMAS_initialization.nc

References:

https://code.zmaw.de/boards/1/topics/996

http://research.jisao.washington.edu/data_sets/nco/#example5

2017年4月24日月曜日

ncdump anaconda

Solved: Not being able to use "ncdump"

Solution: Execute the following two commands:

conda remove hdf4 hdf5 libnetcdf netcdf4

conda install -c conda-forge netcdf4

Ref: https://github.com/spencerahill/aospy/issues/58

2017年4月19日水曜日

前々から気になっていた目の症状があって目のスペシャリストに見てもらったら

マイグレイン・オーラ

と診断された。

徹夜で仕事をし始めた学士論文ごろから起きるようになって

寝不足・コーヒー・チョコレートなども要因らしい。

目の病気だと思って心配していたのでよかった。

昨日の夜も９：３０から３０分ぐらい症状が出た。

起きるのは年に２〜３度。

日本語では 閃輝暗点（せんきあんてん）

というんだね。

2017年4月17日月曜日

LIM2 read air temperature to decide if snow or rain

CORE2で行くのでもう必要ないけどもしいつか必要になった時のために残しておくことにする。

2017年4月12日水曜日

random notes

one way to make existing "time" dimension to be "unlimited":

cdo seltimestep,1/12 ifile ofile

where 1/12 is the first and last indices of your time record.

Ref: https://code.zmaw.de/boards/1/topics/798

 1057  cdo seltimestep,1/1 sic_G10010_SIBT1850_v1.1-naa_icecon.nc4 sic.nc
 1058  cdo seltimestep,1/1 sit_G10010_SIBT1850_v1.1-naa_icecon.nc4 sic.nc
 1059  cdo seltimestep,1/1 sit_PIOMAS_1979_2013-naa_arlan.nc4 sit.nc

Add all variables in sic.nc to sit.nc
 1062  ncks -h -A sic.nc sit.nc

Remove time_counter dimension from the lists of dimensions and variables. 
 1065  ncwa -a time_counter sit.nc sit2.nc

Same names (labels) but different IDs for nodes in PyGraphviz:

add_node(1,label=A)

add_node(2,label=B)

add_node(3,label=A)

will create three nodes consisting of two As and one B.

Ref: http://stackoverflow.com/questions/15711909/how-to-create-duplicate-nodes-in-pygraphviz