Changeset 12590 for NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/traisf.F90

Timestamp:

2020-03-23T22:16:19+01:00 (4 years ago)

Author:

techene

Message:

all: add e3 substitute, OCE/DOM/domzgr_substitute.h90: correct a bug for e3f

File:

• NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/traisf.F90 (modified) (6 diffs)

NEMO/branches/2020/dev_r12377_KERNEL-06_techene_e3/src/OCE/TRA/traisf.F90

@@ -33 +33 @@
       !!----------------------------------------------------------------------
       !!                  ***  ROUTINE tra_isf  ***
-      !!                   
+      !!
       !! ** Purpose :  Compute the temperature trend due to the ice shelf melting (qhoce + qhc)
       !!
@@ -61 +61 @@
          !
          ! Dynamical stability at start up after change in under ice shelf cavity geometry is achieve by correcting the divergence.
-         ! This is achieved by applying a volume flux in order to keep the horizontal divergence after remapping 
+         ! This is achieved by applying a volume flux in order to keep the horizontal divergence after remapping
          ! the same as at the end of the latest time step. So correction need to be apply at nit000 (euler time step) and
          ! half of it at nit000+1 (leap frog time step).
@@ -89 +89 @@
       !! *** Purpose :  Compute the temperature trend due to the ice shelf melting (qhoce + qhc) for cav or par case
       !!
-      !! *** Action :: Update pts(:,:,:,:,Krhs) with the surface boundary condition trend 
+      !! *** Action :: Update pts(:,:,:,:,Krhs) with the surface boundary condition trend
       !!
       !!----------------------------------------------------------------------
@@ -98 +98 @@
       REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in   ) :: ptsc , ptsc_b
       !!----------------------------------------------------------------------
-      INTEGER                      :: ji,jj,jk  ! loop index   
+      INTEGER                      :: ji,jj,jk  ! loop index
       INTEGER                      :: ikt, ikb  ! top and bottom level of the tbl
       REAL(wp), DIMENSION(jpi,jpj) :: ztc       ! total ice shelf tracer trend
@@ -117 +117 @@
          END DO
          !
-         ! level partially include in ice shelf boundary layer 
+         ! level partially include in ice shelf boundary layer
          pts(ji,jj,ikb,jp_tem) = pts(ji,jj,ikb,jp_tem) + ztc(ji,jj) * pfrac(ji,jj)
          !
@@ -128 +128 @@
       !!                  ***  ROUTINE tra_isf_cpl  ***
       !!
-      !! *** Action :: Update pts(:,:,:,:,Krhs) with the ice shelf coupling trend 
+      !! *** Action :: Update pts(:,:,:,:,Krhs) with the ice shelf coupling trend
       !!
       !!----------------------------------------------------------------------