Changeset 12631

Timestamp:

2020-03-31T10:55:26+02:00 (4 years ago)

Author:

jchanut

Message:

Solve #2400 and replace at the same time useless area array by e1e2t

File:

• NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diaar5.F90 (modified) (11 diffs)

NEMO/releases/r4.0/r4.0-HEAD/src/OCE/DIA/diaar5.F90

@@ -32 +32 @@
    REAL(wp)                         ::   vol0         ! ocean volume (interior domain)
    REAL(wp)                         ::   area_tot     ! total ocean surface (interior domain)
-   REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:  ) ::   area         ! cell surface (interior domain)
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:  ) ::   thick0       ! ocean thickness (interior domain)
    REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) ::   sn0          ! initial salinity
@@ -54 +53 @@
       !!----------------------------------------------------------------------
       !
-      ALLOCATE( area(jpi,jpj), thick0(jpi,jpj) , sn0(jpi,jpj,jpk) , STAT=dia_ar5_alloc )
+      ALLOCATE( thick0(jpi,jpj) , sn0(jpi,jpj,jpk) , STAT=dia_ar5_alloc )
       !
       CALL mpp_sum ( 'diaar5', dia_ar5_alloc )
@@ -89 +88 @@
          ALLOCATE( zrhd(jpi,jpj,jpk) , zrhop(jpi,jpj,jpk) )
          ALLOCATE( ztsn(jpi,jpj,jpk,jpts) )
-         zarea_ssh(:,:) = area(:,:) * sshn(:,:)
-      ENDIF
-      !
-      CALL iom_put( 'e2u'      , e2u (:,:) )
-      CALL iom_put( 'e1v'      , e1v (:,:) )
-      CALL iom_put( 'areacello', area(:,:) )
+         zarea_ssh(:,:) = e1e2t(:,:) * sshn(:,:)
+      ENDIF
+      !
+      CALL iom_put( 'e2u'      , e2u  (:,:) )
+      CALL iom_put( 'e1v'      , e1v  (:,:) )
+      CALL iom_put( 'areacello', e1e2t(:,:) )
       !
       IF( iom_use( 'volcello' ) .OR. iom_use( 'masscello' )  ) THEN  
          zrhd(:,:,jpk) = 0._wp        ! ocean volume ; rhd is used as workspace
          DO jk = 1, jpkm1
-            zrhd(:,:,jk) = area(:,:) * e3t_n(:,:,jk) * tmask(:,:,jk)
+            zrhd(:,:,jk) = e1e2t(:,:) * e3t_n(:,:,jk) * tmask(:,:,jk)
          END DO
          CALL iom_put( 'volcello'  , zrhd(:,:,:)  )  ! WARNING not consistent with CMIP DR where volcello is at ca. 2000
@@ -152 +151 @@
          END IF
          !                                         
-         zarho = glob_sum( 'diaar5', area(:,:) * zbotpres(:,:) ) 
+         zarho = glob_sum( 'diaar5', e1e2t(:,:) * zbotpres(:,:) ) 
          zssh_steric = - zarho / area_tot
          CALL iom_put( 'sshthster', zssh_steric )
@@ -178 +177 @@
          END IF
          !    
-         zarho = glob_sum( 'diaar5', area(:,:) * zbotpres(:,:) ) 
+         zarho = glob_sum( 'diaar5', e1e2t(:,:) * zbotpres(:,:) ) 
          zssh_steric = - zarho / area_tot
          CALL iom_put( 'sshsteric', zssh_steric )
@@ -194 +193 @@
              DO jj = 1, jpj
                 DO ji = 1, jpi
-                   zztmp = area(ji,jj) * e3t_n(ji,jj,jk)
+                   zztmp = e1e2t(ji,jj) * e3t_n(ji,jj,jk)
                    ztsn(ji,jj,1,jp_tem) = ztsn(ji,jj,1,jp_tem) + zztmp * tsn(ji,jj,jk,jp_tem)
                    ztsn(ji,jj,1,jp_sal) = ztsn(ji,jj,1,jp_sal) + zztmp * tsn(ji,jj,jk,jp_sal)
@@ -242 +241 @@
                z2d(:,:) = 0._wp
                DO jk = 1, jpkm1
-                 z2d(:,:) = z2d(:,:) + area(:,:) * e3t_n(:,:,jk) * ztpot(:,:,jk)
+                 z2d(:,:) = z2d(:,:) + e1e2t(:,:) * e3t_n(:,:,jk) * ztpot(:,:,jk)
                END DO
                ztemp = glob_sum( 'diaar5', z2d(:,:)  ) 
@@ -249 +248 @@
              !
              IF( iom_use( 'ssttot' ) ) THEN   ! Output potential temperature in case we use TEOS-10
-               zsst = glob_sum( 'diaar5',  area(:,:) * ztpot(:,:,1)  ) 
+               zsst = glob_sum( 'diaar5',  e1e2t(:,:) * ztpot(:,:,1)  ) 
                CALL iom_put( 'ssttot', zsst / area_tot )
              ENDIF
@@ -268 +267 @@
       ELSE       
          IF( iom_use('ssttot') ) THEN   ! Output sst in case we use EOS-80
-            zsst  = glob_sum( 'diaar5', area(:,:) * tsn(:,:,1,jp_tem) )
+            zsst  = glob_sum( 'diaar5', e1e2t(:,:) * tsn(:,:,1,jp_tem) )
             CALL iom_put('ssttot', zsst / area_tot )
          ENDIF
@@ -400 +399 @@
          IF( dia_ar5_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'dia_ar5_init : unable to allocate arrays' )
 
-         area(:,:) = e1e2t(:,:)
-         area_tot  = glob_sum( 'diaar5', area(:,:) )
+         area_tot  = glob_sum( 'diaar5', e1e2t(:,:) )
 
          ALLOCATE( zvol0(jpi,jpj) )
@@ -410 +408 @@
                DO ji = 1, jpi
                   idep = tmask(ji,jj,jk) * e3t_0(ji,jj,jk)
-                  zvol0 (ji,jj) = zvol0 (ji,jj) +  idep * area(ji,jj)
+                  zvol0 (ji,jj) = zvol0 (ji,jj) +  idep * e1e2t(ji,jj)
                   thick0(ji,jj) = thick0(ji,jj) +  idep    
                END DO