Changeset 5835 for branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/iscplhsb.F90

Timestamp:

2015-10-26T15:41:05+01:00 (9 years ago)

Author:

mathiot

Message:

ice sheet coupling branche: cosmetic changes

File:

• branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/iscplhsb.F90 (modified) (3 diffs)

branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DOM/iscplhsb.F90

@@ -82 +82 @@
       pvol_flx(:,:,:  ) = 0.0_wp
       pts_flx (:,:,:,:) = 0.0_wp
+      r1_tiscpl = 1._wp / (prdt_iscpl * rn_rdt) 
 
       ! mask tsn and tsb 
@@ -87 +88 @@
       tsb(:,:,:,jp_sal)=tsb(:,:,:,jp_sal)*ptmask_b(:,:,:); tsn(:,:,:,jp_sal)=tsn(:,:,:,jp_sal)*tmask(:,:,:);
 
-      ! diagnose non conservation of heat, salt and volume 
-      r1_tiscpl = 1._wp / (prdt_iscpl * rn_rdt) 
-
+      !==============================================================================
+      ! diagnose the heat, salt and volume input and compute the correction variable
+      !==============================================================================
+
+      ! 
       zssh0(:,:)        = sshn(:,:) * ssmask(:,:) - sshb(:,:) * psmask_b(:,:)
       IF ( lk_vvl ) zssh0 = 0.0_wp ! already include in the levels by definition
-
+      
       DO jk = 1,jpk-1
          DO ji = 2,jpi-1
@@ -119 +122 @@
                   pts_flx (ji,jj,jk,jp_tem)=   pts_flx (ji,jj,jk,jp_tem) + zdtem(ji,jj) * r1_tiscpl
 
+                  ! case where we close a cell: check if the neighbour cells are wet 
                   IF ( tmask(ji,jj,jk) == 0._wp .AND. ptmask_b(ji,jj,jk) == 1._wp ) THEN
-                     ! case where we close a cell: check if the neighbour cells are wet 
 
                      jip1=ji+1 ; jim1=ji-1 ; jjp1=jj+1 ; jjm1=jj-1 ;