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

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...

File:

• NEMO/branches/2019/dev_r11085_ASINTER-05_Brodeau_Advanced_Bulk/tests/STATION_ASF/MY_SRC/usrdef_hgr.F90 (modified) (3 diffs)

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