Changeset 12475 for NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceupdate.F90

Timestamp:

2020-02-27T12:27:50+01:00 (4 years ago)

Author:

dancopsey

Message:

• Add more print statements.
• Move away from using snow to ice diagnostics and use a new snow to pond one instead.

File:

• NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceupdate.F90 (modified) (4 diffs)

NEMO/branches/UKMO/NEMO_4.0_add_pond_lids_prints/src/ICE/iceupdate.F90

@@ -154 +154 @@
                &           + wfx_opw(ji,jj) + wfx_dyn(ji,jj) + wfx_res(ji,jj) + wfx_lam(ji,jj) + wfx_pnd(ji,jj)
 
+
+
             ! add the snow melt water to snow mass flux to the ocean
             wfx_snw(ji,jj) = wfx_snw_sni(ji,jj) + wfx_snw_dyn(ji,jj) + wfx_snw_sum(ji,jj)
@@ -161 +163 @@
             emp(ji,jj)    = emp_oce(ji,jj) - wfx_ice(ji,jj) - wfx_snw(ji,jj) - wfx_err_sub(ji,jj)   ! mass flux + F/M mass flux (always ice/ocean mass exchange)
 
+
+            
 
             ! Salt flux at the ocean surface      
@@ -166 +170 @@
             sfx(ji,jj) = sfx_bog(ji,jj) + sfx_bom(ji,jj) + sfx_sum(ji,jj) + sfx_sni(ji,jj) + sfx_opw(ji,jj)   &
                &       + sfx_res(ji,jj) + sfx_dyn(ji,jj) + sfx_bri(ji,jj) + sfx_sub(ji,jj) + sfx_lam(ji,jj) + sfx_pnd(ji,jj)
+
+
             
             ! Mass of snow and ice per unit area   
@@ -177 +183 @@
          END DO
       END DO
+
+            write(numout,*)'iceupdate: MEAN(wfx_pnd) = ',SUM(wfx_pnd)/SIZE(wfx_pnd)
+            write(numout,*)'iceupdate: MEAN(wfx_sum) = ',SUM(wfx_sum)/SIZE(wfx_sum)
+            write(numout,*)'iceupdate: MEAN(wfx_bom) = ',SUM(wfx_bom)/SIZE(wfx_bom)
+            write(numout,*)'iceupdate: MEAN(wfx_ice) = ',SUM(wfx_ice)/SIZE(wfx_ice)
+            write(numout,*)'iceupdate: MEAN(wfx_sni) = ',SUM(wfx_sni)/SIZE(wfx_sni)
+            write(numout,*)'iceupdate: MEAN(wfx_snp) = ',SUM(wfx_snp)/SIZE(wfx_snp)
+            write(numout,*)'iceupdate: MEAN(fmmflx) = ',SUM(fmmflx)/SIZE(fmmflx)
+            write(numout,*)'iceupdate: MEAN(wfx_snw) = ',SUM(wfx_snw)/SIZE(wfx_snw)
+            write(numout,*)'iceupdate: MEAN(emp) = ',SUM(emp)/SIZE(emp)
+            write(numout,*)'iceupdate: MEAN(emp_oce) = ',SUM(emp_oce)/SIZE(emp_oce)
+            write(numout,*)'iceupdate: MEAN(sfx) = ',SUM(sfx)/SIZE(sfx)
+            write(numout,*)'iceupdate: MEAN(sfx_sni) = ',SUM(sfx_sni)/SIZE(sfx_sni)
+            write(numout,*)'iceupdate: MEAN(sfx_pnd) = ',SUM(sfx_pnd)/SIZE(sfx_pnd)
+            write(numout,*)'iceupdate: MEAN(sfx_sum) = ',SUM(sfx_sum)/SIZE(sfx_sum)
+            write(numout,*)'iceupdate: MEAN(sfx_bom) = ',SUM(sfx_sum)/SIZE(sfx_bom)
+            write(numout,*)'iceupdate: MAXVAL(wfx_pnd) = ',MAXVAL(wfx_pnd)
+            write(numout,*)'iceupdate: MINVAL(wfx_pnd) = ',MINVAL(wfx_pnd)
+            write(numout,*)'iceupdate: MAXVAL(wfx_bom) = ',MAXVAL(wfx_bom)
+            write(numout,*)'iceupdate: MINVAL(wfx_bom) = ',MINVAL(wfx_bom)
+            write(numout,*)'iceupdate: MAXVAL(wfx_sni) = ',MAXVAL(wfx_sni)
+            write(numout,*)'iceupdate: MINVAL(wfx_sni) = ',MINVAL(wfx_sni)
+            write(numout,*)'iceupdate: MAXVAL(sfx_pnd) = ',MAXVAL(sfx_pnd)
+            write(numout,*)'iceupdate: MINVAL(sfx_pnd) = ',MINVAL(sfx_pnd)
+            write(numout,*)'iceupdate: MAXVAL(sfx_bom) = ',MAXVAL(sfx_bom)
+            write(numout,*)'iceupdate: MINVAL(sfx_bom) = ',MINVAL(sfx_bom)
+
 
       ! Storing the transmitted variables