Changeset 84 for trunk/NEMO/OPA_SRC/SBC/ocesbc.F90

Timestamp:

2004-04-22T15:32:41+02:00 (20 years ago)

Author:

opalod

Message:

CT : UPDATE057 : # General syntax, alignement, comments corrections

# l_ctl alone replace the set (l_ctl .AND. lwp)
# Add of diagnostics which are activated when using l_ctl logical

File:

• trunk/NEMO/OPA_SRC/SBC/ocesbc.F90 (modified) (9 diffs)

trunk/NEMO/OPA_SRC/SBC/ocesbc.F90

@@ -6 +6 @@
 
    !!----------------------------------------------------------------------
-   !!    oce_sbc : initialization and namelist read
+   !!   oce_sbc     : ???
+   !!   oce_sbc_dmp : ???
    !!----------------------------------------------------------------------
    !! * Modules used
@@ -34 +35 @@
 
    !! * Shared module variables
-   REAL(wp), PUBLIC ::   &
+   REAL(wp), PUBLIC ::   &  !:
       aplus, aminus,     &  !:
       empold = 0.e0         !: current year freshwater budget correction
-   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &
+   REAL(wp), PUBLIC, DIMENSION(jpi,jpj) ::   &  !:
       qt  ,         &  !: total surface heat flux (w/m2)
       q   ,         &  !: surface heat flux (w/m2)
@@ -158 +159 @@
 
          ! Re-initialization of fluxes
-         sst_io(:,:) = 0.0
-         sss_io(:,:) = 0.0
-         u_io  (:,:) = 0.0
-         v_io  (:,:) = 0.0
-         gtaux (:,:) = 0.
-         gtauy (:,:) = 0.
+         sst_io(:,:) = 0.e0
+         sss_io(:,:) = 0.e0
+         u_io  (:,:) = 0.e0
+         v_io  (:,:) = 0.e0
+         gtaux (:,:) = 0.e0
+         gtauy (:,:) = 0.e0
 
       ENDIF
@@ -257 +258 @@
 
          ! Re-initialization of fluxes
-         sst_io(:,:) = 0.0
-         sss_io(:,:) = 0.0
-         u_io  (:,:) = 0.0
-         v_io  (:,:) = 0.0
+         sst_io(:,:) = 0.e0
+         sss_io(:,:) = 0.e0
+         u_io  (:,:) = 0.e0
+         v_io  (:,:) = 0.e0
 
       ENDIF
@@ -651 +652 @@
       DO jj = 1, jpj
          DO ji = 1, jpi
-            freezn(ji,jj) = MAX(0., SIGN(1., freeze(ji,jj)-rsmall))
+            freezn(ji,jj) = MAX(0., SIGN(1., freeze(ji,jj)-rsmall) )
          END DO
       END DO
@@ -683 +684 @@
       END DO
       ! volume flux associated to internal damping to climatology
-      dmp(:,:) = zstrdmp(:,:) * rauw / ( zsss(:,:) + rsmall )
+!!ibu dmp(:,:) = zstrdmp(:,:) * rauw / ( zsss(:,:) + rsmall )
+      dmp(:,:) = zstrdmp(:,:) * rauw / ( zsss(:,:) + 1.e-20 )
 # else
       dmp(:,:) = 0.e0            ! No internal damping
@@ -698 +700 @@
             zerp = ( 1. - 2.*upsrnfh(ji,jj) ) * zsrp   &
                & * ( zsss(ji,jj) - s_dta(ji,jj,1) )     &
-               & / ( zsss(ji,jj) + rsmall        )
+               & / ( zsss(ji,jj) + 1.e-20        )
+!ib            & / ( zsss(ji,jj) + rsmall        )
             
             zerp = MIN( zerp, zplus  )
@@ -719 +722 @@
       IF( lk_mpp )   CALL mpp_sum( aplus  )   ! sums over the global domain
       IF( lk_mpp )   CALL mpp_sum( aminus )
-      IF( l_ctl .AND. lwp ) WRITE(numout,*) ' oce_sbc_dmp : a+ = ', aplus, ' a- = ', aminus
+      IF(l_ctl)   WRITE(numout,*) ' oce_sbc_dmp : a+ = ', aplus, ' a- = ', aminus
 
       zadefi = MIN( aplus, aminus )
-      IF( zadefi == 0.0 ) THEN 
+      IF( zadefi == 0.e0 ) THEN 
          erp(:,:) = 0.e0
       ELSE
@@ -732 +735 @@
       erp(:,:) = ( 1. - zfreeze(:,:) ) * zsrp    &   ! surface restoring term
          &     * ( zsss(:,:) - s_dta(:,:,1) )     &
-         &     / ( zsss(:,:) + rsmall      )
+         &     / ( zsss(:,:) + 1.e-20      )
+!ib      &     / ( zsss(:,:) + rsmall      )
 #endif