Changeset 11930 for NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk

Timestamp:

2019-11-19T17:41:22+01:00 (4 years 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:

• EXPREF/file_def_nemo-oce.xml (modified) (1 diff)
• EXPREF/launch_sasf.sh (added)
• EXPREF/namelist_coare3p6_cfg (modified) (1 diff)
• EXPREF/namelist_ecmwf_cfg (modified) (1 diff)
• EXPREF/namelist_ncar_cfg (modified) (1 diff)
• EXPREF/plot_station_asf.py (added)
• MY_SRC/step_c1d.F90 (modified) (2 diffs)
• MY_SRC/stpctl.F90 (added)
• MY_SRC/usrdef_hgr.F90 (modified) (3 diffs)
• MY_SRC/usrdef_nam.F90 (modified) (1 diff)
• README.md (modified) (1 diff)

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/file_def_nemo-oce.xml

@@ -16 +16 @@
       <field field_ref="dt_skin"      name="dt_skin"  />
       <!--    -->
-      <field field_ref="qt_oce"       name="qt_oce"   />
-      <field field_ref="qsr_oce"      name="qsr_oce"  />
-      <field field_ref="qns_oce"      name="qns_oce"  />
-      <field field_ref="qla_oce"      name="qla_oce"  />
-      <field field_ref="evap_oce"     name="evap_oce" />
-      <field field_ref="qsb_oce"      name="qsb_oce"  />
-      <field field_ref="qlw_oce"      name="qlw_oce"  />
+      <field field_ref="qt_oce"       name="qt"   />
+      <field field_ref="qsr_oce"      name="qsr"  />
+      <field field_ref="qns_oce"      name="qns"  />
+      <field field_ref="qla_oce"      name="qla"  />
+      <field field_ref="qsb_oce"      name="qsb"  />
+      <field field_ref="qlw_oce"      name="qlw"  />
+      <field field_ref="qemp_oce"     name="qemp" />
       <!-- -->
-      <field field_ref="qemp_oce"     name="qemp_oce" />
-      <field field_ref="precip"       name="precip"   />
+      <field field_ref="evap_oce"     name="evap" />
+      <field field_ref="precip"       name="precip" />
+      <field field_ref="empmr"        name="empmr" />
       <!-- -->
       <field field_ref="taum"         name="taum"     />

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_coare3p6_cfg

@@ -277 +277 @@
 &namctl        !   Control prints                                       (default: OFF)
 !-----------------------------------------------------------------------
-ln_ctl = .false. !lolo
+ln_ctl = .false.
 /
 !-----------------------------------------------------------------------

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_ecmwf_cfg

@@ -277 +277 @@
 &namctl        !   Control prints                                       (default: OFF)
 !-----------------------------------------------------------------------
-ln_ctl = .false. !lolo
+ln_ctl = .false.
 /
 !-----------------------------------------------------------------------

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/EXPREF/namelist_ncar_cfg

@@ -277 +277 @@
 &namctl        !   Control prints                                       (default: OFF)
 !-----------------------------------------------------------------------
-ln_ctl = .false. !lolo
+ln_ctl = .false.
 /
 !-----------------------------------------------------------------------

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/step_c1d.F90

@@ -67 +67 @@
       ! Control
       !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
-      !LOLO:                       
-                            ! LOLO: should write a stp_ctl, who checks on flux realism (current SAS version checks ice fields...)
-                            ! CALL stp_ctl( kstp, indic )
-      !IF( indic < 0  )  THEN
-      !                       CALL ctl_stop( 'step: indic < 0' )
-      !                       CALL dia_wri_state( 'output.abort' )
-      !ENDIF
-      !LOLO.
+                             CALL stp_ctl( kstp, indic )
+      IF( indic < 0  )  THEN
+                             CALL ctl_stop( 'step: indic < 0' )
+                             CALL dia_wri_state( 'output.abort' )
+      ENDIF
                              
       IF( kstp == nit000   ) CALL iom_close( numror )     ! close input  ocean restart file
 
-      !LOLO: from C1D:
+      !#LB: from C1D:
       IF( lrst_oce       )   CALL rst_write( kstp )        ! write output ocean restart file
       !
@@ -85 +82 @@
 #endif
 
-      !LOLO: from SAS"
+      !#LB: from SAS"
 !#if defined key_iomput
 !      IF( kstp == nitrst ) THEN

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/usrdef_hgr.F90

@@ -67 +67 @@
       REAL(wp), DIMENSION(:,:), INTENT(out) ::   pe1e2u, pe1e2v               ! u- & v-surfaces (if reduction in strait)   [m2]
       !
-      INTEGER  ::   ji, jj               ! dummy loop indices
-      REAL(wp) ::   zlam1, zlam0, zcos_alpha, zim1 , zjm1 , ze1  , ze1deg, zf0   ! local scalars
-      REAL(wp) ::   zphi1, zphi0, zsin_alpha, zim05, zjm05, zbeta, znorme, zfact !   -      -
+      REAL(wp) :: zf0   ! local scalar
       !!-------------------------------------------------------------------------------
       !
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) 'usr_def_hgr : C1D configuration'
+      IF(lwp) WRITE(numout,*) 'usr_def_hgr : STATION_ASF configuration'
       IF(lwp) WRITE(numout,*) '~~~~~~~~~~~   uniform grid spacing WITH Coriolis force'
-      IF(lwp) WRITE(numout,*) ' LOLO => rn_lon1d, rn_lat1d =', rn_lon1d, rn_lat1d
+      IF(lwp) WRITE(numout,*) '              coordinates of station =>', rn_lon1d, rn_lat1d
       !
       !                       !==  grid point position  ==!   (Special case of 1D config: all lon/lat are the same)
@@ -88 +86 @@
       pphiv(:,:) = rn_lat1d
       pphif(:,:) = rn_lat1d
+
       !                       !== Horizontal scale factors ==! (in meters)
-      !                     
       !                                         ! constant grid spacing
       pe1t(:,:) = 100.  ;   pe2t(:,:) = 100.
@@ -103 +101 @@
       !
       !                       !==  Coriolis parameter  ==!
-      kff = 0                 !  indicate to compute Coriolis parameter afterward
+      zf0   = 2._wp * omega * SIN( rad * rn_lat1d )      
+      pff_f(:,:) = zf0
+      pff_t(:,:) = zf0      
+      kff = 1                 !  indicate to skip computing Coriolis parameter afterward
       !
    END SUBROUTINE usr_def_hgr

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/usrdef_nam.F90

@@ -8 +8 @@
    !!======================================================================
    !! History :  4.0  ! 2016-03  (S. Flavoni, G. Madec)  Original code
-   !! History :  trunk  ! 2019-10 (L. Brodeau) for STATION_ASF (C1D meets SAS)
+   !! History :  4.x  ! 2019-10  (L. Brodeau) for STATION_ASF (C1D meets SAS)
    !!----------------------------------------------------------------------
 

NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/README.md

@@ -64 +64 @@
 
 
+## Playing with STATION_ASF
 
+First compile the test-case as follows (compile with xios-2.5 support → check your ARCH file):
 
+```./makenemo -m <your_arch> -n STATION_ASF -j 4 -a STATION_ASF```
 
+Then 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:
 
+- ```NEMO_DIR``` : NEMO root directory where to fetch compiled STATION_ASF ```nemo.exe``` + setup (such as ```${NEMO_DIR}/tests/STATION_ASF```)
 
+- ```WORK_DIR``` :  Directory where to run the simulation
 
+- ```FORC_DIR```  Directory containing sea-surface + atmospheric forcings (get it there https://drive.google.com/file/d/1MxNvjhRHmMrL54y6RX7WIaM9-LGl--ZP/)
+