Changeset 11930


Ignore:
Timestamp:
2019-11-19T17:41:22+01:00 (11 months ago)
Author:
laurent
Message:

New functional "stp_ctl" for STATION_ASF, test-case launch script and python diagnostics…

Location:
NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF
Files:
3 added
8 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/file_def_nemo-oce.xml

    r11861 r11930  
    1616      <field field_ref="dt_skin"      name="dt_skin"  /> 
    1717      <!--    --> 
    18       <field field_ref="qt_oce"       name="qt_oce"   /> 
    19       <field field_ref="qsr_oce"      name="qsr_oce"  /> 
    20       <field field_ref="qns_oce"      name="qns_oce"  /> 
    21       <field field_ref="qla_oce"      name="qla_oce"  /> 
    22       <field field_ref="evap_oce"     name="evap_oce" /> 
    23       <field field_ref="qsb_oce"      name="qsb_oce"  /> 
    24       <field field_ref="qlw_oce"      name="qlw_oce" /> 
     18      <field field_ref="qt_oce"       name="qt"   /> 
     19      <field field_ref="qsr_oce"      name="qsr"  /> 
     20      <field field_ref="qns_oce"      name="qns"  /> 
     21      <field field_ref="qla_oce"      name="qla"  /> 
     22      <field field_ref="qsb_oce"      name="qsb" /> 
     23      <field field_ref="qlw_oce"      name="qlw"  /> 
     24      <field field_ref="qemp_oce"     name="qemp" /> 
    2525      <!-- --> 
    26       <field field_ref="qemp_oce"     name="qemp_oce" /> 
    27       <field field_ref="precip"       name="precip"   /> 
     26      <field field_ref="evap_oce"     name="evap" /> 
     27      <field field_ref="precip"       name="precip" /> 
     28      <field field_ref="empmr"        name="empmr" /> 
    2829      <!-- --> 
    2930      <field field_ref="taum"         name="taum"     /> 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_coare3p6_cfg

    r11905 r11930  
    277277&namctl        !   Control prints                                       (default: OFF) 
    278278!----------------------------------------------------------------------- 
    279 ln_ctl = .false. !lolo 
     279ln_ctl = .false. 
    280280/ 
    281281!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_ecmwf_cfg

    r11905 r11930  
    277277&namctl        !   Control prints                                       (default: OFF) 
    278278!----------------------------------------------------------------------- 
    279 ln_ctl = .false. !lolo 
     279ln_ctl = .false. 
    280280/ 
    281281!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_ncar_cfg

    r11905 r11930  
    277277&namctl        !   Control prints                                       (default: OFF) 
    278278!----------------------------------------------------------------------- 
    279 ln_ctl = .false. !lolo 
     279ln_ctl = .false. 
    280280/ 
    281281!----------------------------------------------------------------------- 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/step_c1d.F90

    r11637 r11930  
    6767      ! Control 
    6868      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    69       !LOLO:                        
    70                             ! LOLO: should write a stp_ctl, who checks on flux realism (current SAS version checks ice fields...) 
    71                             ! CALL stp_ctl( kstp, indic ) 
    72       !IF( indic < 0  )  THEN 
    73       !                       CALL ctl_stop( 'step: indic < 0' ) 
    74       !                       CALL dia_wri_state( 'output.abort' ) 
    75       !ENDIF 
    76       !LOLO. 
     69                             CALL stp_ctl( kstp, indic ) 
     70      IF( indic < 0  )  THEN 
     71                             CALL ctl_stop( 'step: indic < 0' ) 
     72                             CALL dia_wri_state( 'output.abort' ) 
     73      ENDIF 
    7774                              
    7875      IF( kstp == nit000   ) CALL iom_close( numror )     ! close input  ocean restart file 
    7976 
    80       !LOLO: from C1D: 
     77      !#LB: from C1D: 
    8178      IF( lrst_oce       )   CALL rst_write( kstp )        ! write output ocean restart file 
    8279      ! 
     
    8582#endif 
    8683 
    87       !LOLO: from SAS" 
     84      !#LB: from SAS" 
    8885!#if defined key_iomput 
    8986!      IF( kstp == nitrst ) THEN 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/usrdef_hgr.F90

    r11637 r11930  
    6767      REAL(wp), DIMENSION(:,:), INTENT(out) ::   pe1e2u, pe1e2v               ! u- & v-surfaces (if reduction in strait)   [m2] 
    6868      ! 
    69       INTEGER  ::   ji, jj               ! dummy loop indices 
    70       REAL(wp) ::   zlam1, zlam0, zcos_alpha, zim1 , zjm1 , ze1  , ze1deg, zf0   ! local scalars 
    71       REAL(wp) ::   zphi1, zphi0, zsin_alpha, zim05, zjm05, zbeta, znorme, zfact !   -      - 
     69      REAL(wp) :: zf0   ! local scalar 
    7270      !!------------------------------------------------------------------------------- 
    7371      ! 
    7472      IF(lwp) WRITE(numout,*) 
    75       IF(lwp) WRITE(numout,*) 'usr_def_hgr : C1D configuration' 
     73      IF(lwp) WRITE(numout,*) 'usr_def_hgr : STATION_ASF configuration' 
    7674      IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   uniform grid spacing WITH Coriolis force' 
    77       IF(lwp) WRITE(numout,*) ' LOLO => rn_lon1d, rn_lat1d =', rn_lon1d, rn_lat1d 
     75      IF(lwp) WRITE(numout,*) '              coordinates of station =>', rn_lon1d, rn_lat1d 
    7876      ! 
    7977      !                       !==  grid point position  ==!   (Special case of 1D config: all lon/lat are the same) 
     
    8886      pphiv(:,:) = rn_lat1d 
    8987      pphif(:,:) = rn_lat1d 
     88 
    9089      !                       !== Horizontal scale factors ==! (in meters) 
    91       !                      
    9290      !                                         ! constant grid spacing 
    9391      pe1t(:,:) = 100.  ;   pe2t(:,:) = 100. 
     
    103101      ! 
    104102      !                       !==  Coriolis parameter  ==! 
    105       kff = 0                 !  indicate to compute Coriolis parameter afterward 
     103      zf0   = 2._wp * omega * SIN( rad * rn_lat1d )       
     104      pff_f(:,:) = zf0 
     105      pff_t(:,:) = zf0       
     106      kff = 1                 !  indicate to skip computing Coriolis parameter afterward 
    106107      ! 
    107108   END SUBROUTINE usr_def_hgr 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/usrdef_nam.F90

    r11831 r11930  
    88   !!====================================================================== 
    99   !! History :  4.0  ! 2016-03  (S. Flavoni, G. Madec)  Original code 
    10    !! History :  trunk  ! 2019-10 (L. Brodeau) for STATION_ASF (C1D meets SAS) 
     10   !! History :  4.x  ! 2019-10 (L. Brodeau) for STATION_ASF (C1D meets SAS) 
    1111   !!---------------------------------------------------------------------- 
    1212 
  • NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/README.md

    r11892 r11930  
    6464 
    6565 
     66## Playing with STATION_ASF 
    6667 
     68First compile the test-case as follows (compile with xios-2.5 support → check your ARCH file): 
    6769 
     70```./makenemo -m <your_arch> -n STATION_ASF -j 4 -a STATION_ASF``` 
    6871 
     72Then you can use the script ``launch_sasf.sh`` found in  ```EXPREF/``` to launch 3 simulations (one for each bulk parameterization available). You need to adapt the following variable to your environment in the script: 
    6973 
     74- ```NEMO_DIR``` : NEMO root directory where to fetch compiled STATION_ASF ```nemo.exe``` + setup (such as ```${NEMO_DIR}/tests/STATION_ASF```) 
    7075 
     76- ```WORK_DIR``` :  Directory where to run the simulation 
    7177 
     78- ```FORC_DIR```  Directory containing sea-surface + atmospheric forcings (get it there https://drive.google.com/file/d/1MxNvjhRHmMrL54y6RX7WIaM9-LGl--ZP/) 
     79 
Note: See TracChangeset for help on using the changeset viewer.