New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
user/flavoni/libIGCM/iom_put (diff) – NEMO

Changes between Version 49 and Version 50 of user/flavoni/libIGCM/iom_put


Ignore:
Timestamp:
2009-11-27T15:24:59+01:00 (14 years ago)
Author:
flavoni
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • user/flavoni/libIGCM/iom_put

    v49 v50  
    1919[[BR]] 
    2020 
    21  '''today 30/10/2009 if you download ./model NEMO_IGCM you can't use iom_put'''. [[BR]] 
     2127/11/2009: downloading NEMO_DEV you have iomput by default (tag v3_2) 
    2222 
    23 === If you want activate iom_put === 
     23=== By default iom_put is activated === 
    2424  
    25  If you want to activate iom_put you have to: 
    26  
    27  * download model NEMO_IGCM: [[BR]] 
    28  {{{ 
    29  cd modipsl/util  
    30  ./model NEMO_IGCM  
    31  }}} 
    32  
    33  
    34  * comment lines in dommsk: [[BR]] 
    35  {{{ 
    36         ij0 = 115   ;   ij1 = 115  
    37          ii0 = 145   ;   ii1 = 146   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 4.0e0  
    38          ij0 = 116   ;   ij1 = 116  
    39          ii0 = 145   ;   ii1 = 146   ;   fmask( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 1:jpk ) = 4.0e0  
    40   }}} 
    41  
    42  * put in namelist 
    43  {{{ 
    44   rn_bfri1    =    0.e-4  
    45  }}} 
    46  
    47   * Add key in modipsl/config/ORCA2_LIM/scripts/BB_make.ldef 
     25  *  key in modipsl/config/ORCA2_LIM/scripts/BB_make.ldef 
    4826 {{{ 
    4927 key_iomput 
    5028 }}}  
    5129 
    52  * Do an ins_make : [[BR]] 
    53  {{{ 
    54  ./modipsl/util/ins_make 
    55  }}} 
    56   
    57  * Go in ORCA2_LIM directory : [[BR]] 
    58  {{{ 
    59  cd modipsl/config/ORCA2_LIM 
    60  }}} 
    6130 
    62  * Compile : [[BR]] 
    63  {{{ 
    64  gmake 
    65  }}} 
     31 * namelist in IGCM00 is the same that namelist in EXP00 (nothing to change if you want to use iodef.xml by default, like an example) 
    6632 
    67  * you can copy modified files (opa9.driver, opa9.card, lim2.card, POST directory) that you find in: [[BR]] 
    68  {{{ 
    69  brodie: /workdir/rech/eee/reee605/LIBIGCM_NEMO/LIKE-COUPLED/test_claire/   
    70  mercure: /work/cont003/flavoni/LIBIGCM_NEMO/LIKE-COUPLED/test_iom_lib14_rev1661_WR/modipsl/config/ORCA2_LIM/IGCM_lim14_io/COMP   
    71  }}} 
     33 * iodef.xml in IGCM00 is the same that namelist in EXP00 (nothing to change if you want to use iodef.xml by default, like an example) 
    7234 
    73  * Change in opa9.card: [[BR]] 
    74  {{{ 
    75  io_xml='''y'''   
    76  }}} 
    77   
    78 ---- 
    79  
    80  === My job: modifications done till 30/10/2009 to prepare files for using iom_put: === 
    81  
    82  Modifications done till now were necessary if not it was not possible to have good file for iodef.xml 
    83  
    84  1. Changed revision for nemo in modispl/util/mod.def in NEMO_IGCM: 
    85  
    86  #-H- NEMO_IGCM  NEMO trunk revision '''1661''' [[BR]]  (revision without bottom friction bug) 
    87   
    88  
    89  in '''lim2.card'''   
    90  
    91  2. changed PREFIX name for output .nc file  [[BR]] 
    92  [OutputFiles]  [[BR]] 
    93  List= ('''${PREFIX_NEMO}'''_icemod.nc, ${R_OUT_ICE_NWRITE}/${PREFIX}_${WF1}_icemod.nc, Post_1M_icemod) 
    94   
    95  3. removed deptht (for monitoring, no good to do ncrcat in 3 dimensions for ice)   [[BR]] 
    96  [Post_1M_icemod] [[BR]] 
    97  GatherWithInternal = (nav_lon, nav_lat, time_counter) 
    98  
    99  4. removed ileadfra '''for the moment'''  [[BR]] 
    100  TimeSeriesVars = (iicethic, isnowthi, iicevelu, iicevelv) 
    101  
    102  NOTE: for this moment ileadfra removed because in iomdef.xml there are new variables' names (like IPPC requested) and ileadfra there is no more in output from nemo's code, waiting for that nemo put in output 1-soicecov (=ileadfra) 
     35 * To have post traitement different by default you need to create POST directory in IGCM00 [[BR]] 
    10336 
    10437 
    105  in '''opa9.card''' 
    106  
    107  5. changed PREFIX name for output .nc file  [[BR]] 
    108  [OutputFiles]  [[BR]] 
    109  List= ('''${PREFIX_NEMO}''''_grid_T.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_T.nc, Post_1M_grid_T),\ [[BR]] 
    110  ('''${PREFIX_NEMO}'''_grid_U.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_U.nc, Post_1M_grid_U),\ [BR]] 
    111  ('''${PREFIX_NEMO}'''_grid_V.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_V.nc, Post_1M_grid_V),\ [[BR]] 
    112  ('''${PREFIX_NEMO}'''_grid_W.nc, ${R_OUT_OCE_NWRITE}/${PREFIX}_${WF1}_grid_W.nc, Post_1M_grid_W),\ [[BR]] 
    113  
    114  (in opa9.driver ${PREFIX_NEMO} is:  [[BR]] 
    115  ${PREFIX_NEMO}=${PREFIX}_{WF1}  if we use iomput [[BR]] 
    116  ${PREFIX_NEMO}=${PREFIX_NWRITE}_${DATE_OPA} if not )   [[BR]] 
    117  
    118  in '''opa9.driver'''  
    119  
    120  6. possibility to use io_server (but not testing now) [[BR]] 
    121  ##--Variables used for XMLIO_SERVER configuration file  [[BR]] 
    122  # using_server: .TRUE.(.FALSE.) to use(bypass) the io_server  [[BR]] 
    123  PAT_SERVER=$( supergrep using_server ${SUBMIT_DIR}/PARAM/xmlio_server.def )  [[BR]] 
    124  
    125  
    126  7. Added output_frequency in days, if not xml file is not correct  [[BR]] 
    127  (( output_frequency_oce = $WriteFrequencyInDays * 86400 )) [[BR]] 
    128  (( output_frequency_ice = $WriteFrequencyInDays * 86400 ))  
    129  
    130  8. replaced output filenames in xml file [[BR]] 
    131  sed -e "s/<file_name_grid_T>/ \"${PREFIX_NEMO}_grid_T.nc\"/" \ [[BR]] 
    132  -e "s/<file_name_grid_U>/ \"${PREFIX_NEMO}_grid_U.nc\"/" \ [[BR]] 
    133  -e "s/<file_name_grid_V>/ \"${PREFIX_NEMO}_grid_V.nc\"/" \ [[BR]] 
    134  -e "s/<file_name_grid_W>/ \"${PREFIX_NEMO}_grid_W.nc\"/" \ [[BR]] 
    135  -e "s/<file_name_icemod>/ \"${PREFIX_NEMO}_icemod.nc\"/" \ [[BR]] 
    136  -e "s/<output_frequency_oce>/ \"${output_frequency_oce}\"/g" \ [[BR]] 
    137  -e "s/<output_frequency_ice>/ \"${output_frequency_ice}\"/g" \ [[BR]] 
    138  iodef.xml > iodef.xml.tmp  
    139  
    140  IGCM_sys_Mv iodef.xml.tmp iodef.xml [[BR]] 
    141  
    142  
    143  9. determine number of day(s) in WriteFrequency :  [[BR]] 
    144    case ${config_UserChoices_OCE_WriteFrquency} in 
    145             1Y|1y) 
    146                 (( WriteFrequencyInDays = $( IGCM_date_DaysInYear $year ) )) ;; 
    147             1M|1m) 
    148                 (( WriteFrequencyInDays = $( IGCM_date_DaysInMonth $year $month ) )) ;; 
    149             5D|5d) 
    150                 (( WriteFrequencyInDays=5 )) ;; 
    151             1D|1d) 
    152                 (( WriteFrequencyInDays=1 )) ;; 
    153             *) 
    154                 IGCM_debug_Exit "IGCM_config_OCE_WriteFrquency " ${config_UserChoices_OCE_WriteFrquency} " invalid write frquency length : choose in 1Y, 1M, 5D, 1D." 
    155                 IGCM_debug_Verif_Exit ;; 
    156    esac 
    157  
    158  10. Created POST directory under $SUBMIT_DIR:  [[BR]] 
    159  In POST directrory there are:   [[BR]] 
     38 Create POST directory under $SUBMIT_DIR:  [[BR]] 
     39 In POST directrory you can put:   [[BR]] 
    16040 monitoring01_lim2_ORCA2.cfg  [[BR]] 
    16141 monitoring01_opa9_ORCA2.cfg  [[BR]] 
    16242   
    163  In monitoring01_lim2_ORCA2.cfg I removed ileadfra because now this variable does not exists, and if I don't use "dynamic monitoring" I can't have ICE monitoring (monitoring don't create monitoring images if there is one or more variables that does not exist).  
    164  
     43 (In monitoring01_lim2_ORCA2.cfg I removed ileadfra because now this variable does not exists)  
    16544 see https://forge.ipsl.jussieu.fr/nemo/wiki/private/libIGCM/postprocessingpage 
    16645 
     
    16847 '''open problems''': 
    16948 1. Don't have ATLAS for "ice cover": because new variable is named soicecov and in atlas is written ileadfra 
    170  2. to do: calcul for 2 frequencies (waiting for a new iodef.xml) 
    171  
    172  
    173 === Tests on vargas === 
    174  
    175  * Directory : 
    176  1- /workgpfs/rech/psl/rpsl003/VARGAS/NEMO/CLIM_IOM/modipsl/config/ORCA2_LIM 
    177  2- /workgpfs/rech/psl/rpsl003/VARGAS/NEMO/CLIM_IOM2/modipsl/config/ORCA2_LIM 
    178  * Main idea :  
    179   * use xml.def to list fields and file for differents types of output (1M, 1D, 1Y, ...) 
    180     * Usefull for grid_[TUVW] and icemod files 
    181     * old fashion for diaptr 
    182   * use -1 for 1Y length simulation to have monthly (variable length) outputs 
    183    * OK for xml files 
    184    * nead adjustment for old fashion (diaptr) files 
    185  * Modifications : 
    186   * BB_make.ldef (key_iomput added) 
    187   * Sources : commented lines in dommsk 
    188   * COMP directory : 
    189    * opa9.driver : 
    190      * use Vxx_enable variable to set information in xml.def file (1M enable for example) 
    191    * opa9.card : 
    192      * add a list of possible files : ORCA2_1m_grid_T.nc ... 
    193      * add a Patch for time counter (to rename the first time axis time_counter, t_ave_01month allowable, works when no tule_counter rename is required) 
    194        * See : /workgpfs/rech/psl/rpsl003/VARGAS/NEMO/CLIM_IOM/modipsl/libIGCM/libIGCM_post/IGCM_Patch_20091030_histcom_time_axis.ksh 
    195      * delete sophtldf variable in monitoring list 
    196    * lim2.card : 
    197      * add a list of possible files : ORCA2_1m_icemod.nc ... 
    198      * add a Patch for time counter (to rename the first time axis time_counter, t_ave_01month allowable, works when no tule_counter rename is required) 
    199      * delete deptht variable in gather list 
    200   * PARAM directory : 
    201    * add iodef.xml 
    202    * xmlio_server.def 
    203   * POST directory : 
    204    * add monitoring01_lim2.cfg file : 
    205     * a copy (take care : without _ORCA2) of /home/rech/psl/rpsl035/atlas/monitoring01_lim2_ORCA2.cfg 
    206     * without ileadfra_south and ileadfra_north lines 
    207  
    208  * 4 tests : All tests have been done with 1Y atlas/monitoring : 
    209   * Modifications in config.card for all tests : 
    210  {{{ 
    211 JobNumProcTot=16 
    212 RebuildFrequency=1Y 
    213 DateEnd=2002-12-31 
    214 TimeSeriesFrequency=1Y 
    215 SeasonalFrequency=1Y 
    216 WriteFrequency="1M 1D" 
    217 WriteFrequency="1M 1D" 
    218 }}} 
    219   * Name of tests 
    220  
    221  || PeriodLength || 1M || 1Y || 
    222  || 360d || NLIO5 / CLIM_IOM/.../IGCM360MO '''OK''' || NLIO7 / CLIM_IOM/.../IGCM360 '''OK''' || 
    223  || leap || NLIO20 /IGCM365MO6 '''KO monitoring''' || NLIO15 / IGCM3656 '''KO impossible for diaptr''' || 
    224  || leap || NLIOD5 / CLIM_IOM2/.../IGCM365MO2 '''OK''' || NLIOD4 '''KO impossible yet''' false for xml files || 
    225  
    226 ==== NLIO5 ==== 
    227  * Status : OK 
    228  * Directory 1 
    229  * Modifications in IGCM360MO/config.card 
    230  {{{ 
    231 JobName=NLIO5 
    232 PeriodLength=1M 
    233 }}} 
    234  * Issues already known : 
    235    * No mask in atlas for continents 
    236    * ileadfra maps empty 
    237  * See : http://dodsp.idris.fr/rpsl003/ORCA2_LIM2/NLIO5 
    238  
    239 ==== NLIO7 ==== 
    240  * Status : OK 
    241  * Directory 1 
    242  * Modifications in IGCM360/config.card 
    243  {{{ 
    244 JobName=NLIO7 
    245 }}} 
    246  * Issues already known : 
    247    * No mask in atlas for continents 
    248    * ileadfra maps empty 
    249  * See : http://dodsp.idris.fr/rpsl003/ORCA2_LIM2/NLIO7 
    250  
    251 ==== NLIO15 ==== 
    252  * Status : KO 
    253  * Directory 1 
    254  * Modifications in IGCM3656/config.card 
    255  {{{ 
    256 JobName=NLIO15 
    257 CalendarType=leap 
    258 }}} 
    259  * Issues : 
    260   * diaptr not saved name : NLIO15_CU*_20010101_20011231_diaptr_0000.nc 
    261    * ZNL atlas empty 
    262    * Monitoring : only ICE variables, only 1 year,  
    263    * will work better when diaptr files will exist 
    264  * Issues already known : 
    265    * No mask in atlas for continents 
    266    * ileadfra maps empty 
    267  * See : http://dodsp.idris.fr/rpsl003/ORCA2_LIM2/NLIO15 
    268  
    269 ==== NLIO20 ==== 
    270  * Status : KO 
    271  * Directory 1 
    272  * Modifications in IGCM365MO6/config.card 
    273  {{{ 
    274 JobName=NLIO20 
    275 CalendarType=leap 
    276 PeriodLength=1M 
    277 }}} 
    278  * Issues : 
    279     * Monitoring : only 1Y for some variables, 2Y for others. ????? 
    280  * Issues already known : 
    281    * No mask in atlas for continents 
    282    * ileadfra maps empty 
    283  * See : http://dodsp.idris.fr/rpsl003/ORCA2_LIM2/NLIO20 
    284  
    285  
    286 ==== NLIOD4 ==== 
    287  * Status : KO 
    288  * Directory 2  ..../NEMO/CLIM_IOM2/modipsl/... 
    289  * libIGCM 195 
    290  * Modifications in IGCM3654/COMP/opa9.driver 
    291  * Modifications in IGCM3654/config.card 
    292  {{{ 
    293 JobName=NLIOD4 
    294 CalendarType=leap 
    295 }}} 
    296  * Issues : 
    297     * Monitoring KO : 
    298       * 22 values in Analyse/TS_MO/* files 
    299       * nb second per month (instead of -1) for diaptr and xml files : grid_[TUVW] and icemod 
    300  * See : http://dodsp.idris.fr/rpsl003/ORCA2_LIM2/NLIOD4 
    301  
    302 ==== NLIOD5 ==== 
    303  * Status : OK 
    304  * Directory 2  ..../NEMO/CLIM_IOM2/modipsl/... 
    305  * libIGCM 195 
    306  * Modifications in IGCM3654/COMP/opa9.driver 
    307  * Modifications in IGCM3654/config.card 
    308  {{{ 
    309 JobName=NLIOD5 
    310 CalendarType=leap 
    311 PeriodLength=1M 
    312 }}} 
    313  * Issues : 
    314     * Monitoring OK : 
    315       * nb second per month (instead of -1) for diaptr and xml files : grid_[TUVW] and icemod 
    316       * Patch OK 
    317  * Issues already known : 
    318    * No mask in atlas for continents 
    319    * ileadfra maps empty 
    320  * See : http://dodsp.idris.fr/rpsl003/ORCA2_LIM2/NLIOD5 
    321  
    322  
    323 ==== Tests without iomput ==== 
    324  
    325 =====  leap calendar ===== 
    326  
    327  * NL365 with 1M Period Length  
    328  * To verify the monitoring produced without iomput. 
    329  * Sources : /workgpfs/rech/psl/rpsl003/VARGAS/NEMO/CLIM/modipsl/config/ORCA2_LIM/NL365  
    330  * See : http://dodsp.idris.fr/rpsl003/ORCA2_LIM2/NL365 
    331  * Same problem on some variables. Only the first period of monitoring is shown. 
    332  
    333 =====  360d calendar ===== 
    334  
    335  * with and without ln_mskland, to verify monitoring 
    336  * NL360M (with ln_mskland=.true.) 
    337  * CLIMV16 (with ln_msklnd=.false.)  
    338  
    339 ==== Proposition 1 ==== 
    340  
    341 Add a directory IGCM01 on svn with key_iomput activated and with this files COMP, PARAM, POST files here : https://forge.ipsl.jussieu.fr/nemo/browser/branches/libIGCM/ORCA2_LIM 
    342  
    343 ==== Proposition 2 ==== 
    344  
    345 Organisation used for IPSLCM5_v2 configuration : 
    346  * Use same directory with this modifications : 
    347   * Add a README file to explain how to switch to noiomput 
    348   * Copy opa9.card/opa9.driver lim2.card xxx.noiomput 
    349   * Adapt opa9.cardopa9.driver lim2.card to iomput 
    350  
    351 === Tests in coupled model IPSLCM5 === 
    352  
    353 More information : http://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM5iomputNEMO 
    35449