Changeset 479

Timestamp:

2006-06-19T15:19:23+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_bugfix_048 : CT : correction of the ssh initial state computation

File:

• trunk/NEMO/OPA_SRC/istate.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/istate.F90

@@ -238 +238 @@
       INTEGER  ::   inum              ! temporary logical unit
       INTEGER  ::   ji, jj, jk        ! dummy loop indices
+      INTEGER  ::   ijloc
       REAL(wp) ::   &
-         zh1, zh2, zslope, zcst       ! temporary scalars
+         zh1, zh2, zslope, zcst,   &  ! temporary scalars
+         zfcor
       REAL(wp) ::   &
          zt1  = 12._wp,            &  ! surface temperature value (EEL R5)
          zt2  =  2._wp,            &  ! bottom  temperature value (EEL R5)
-         zsal = 35.5_wp                ! constant salinity (EEL R2, R5 and R6)
+         zsal = 35.5_wp,           &  ! constant salinity (EEL R2, R5 and R6)
+         zueel = 0.1_wp               ! constant uniform zonal velocity (EEL R5)
 # if ! defined key_dynspg_rl
       REAL(wp), DIMENSION(jpiglo,jpjglo) ::   &
@@ -291 +294 @@
             ! Start EEL5 configuration with barotropic geostrophic velocities 
             ! according the sshb and sshn SSH imposed.
-            ub(:,:,:) = 0.1 * umask(:,:,:)
+            ! we assume a uniform grid (hence the use of e1t(1,1) for delta_y)
+            ! we use the Coriolis frequency at mid-channel.   
+    
+            ub(:,:,:) = zueel * umask(:,:,:)
             un(:,:,:) = ub(:,:,:)
+            ijloc = mj0(INT(jpjglo-1)/2)
+            zfcor = ff(1,ijloc)
 
             DO jj = 1, jpjglo
-               zssh(:,jj) = ( .22 - ( FLOAT(jj-3) * (0.44) ) / 99. )
+               zssh(:,jj) = - (FLOAT(jj)- FLOAT(jpjglo-1)/2.)*zueel*e1t(1,1)*zfcor/grav 
             END DO
+
+            IF(lwp) THEN
+               WRITE(numout,*) ' Uniform zonal velocity for EEL R5:',zueel
+               WRITE(numout,*) ' Geostrophic SSH profile as a function of y:'
+               WRITE(numout,'(12(1x,f6.2))') zssh(1,:)
+            ENDIF
+
             DO jj = 1, nlcj
                DO ji = 1, nlci