Changeset 5802 for branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90

Timestamp:

2015-10-16T19:19:43+02:00 (9 years ago)

Author:

mathiot

Message:

ice sheet coupling: reproducibility fixed in MISOMIP configuration + contrain bathy to be constant during all the run

File:

• branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90 (modified) (1 diff)

branches/NERC/dev_r5589_is_oce_cpl/NEMOGCM/NEMO/OPA_SRC/DIA/diahsb.F90

@@ -200 +200 @@
 !      ENDIF
 !!gm end
-      
+      IF ( lk_vvl ) THEN
+         IF (lwp) PRINT *, 'cons heat : ', kt, zdiff_hc / zvol_tot, zdiff_sc / zvol_tot 
+         IF (lwp) PRINT *, 'cons volu : ', kt, zdiff_v2 * 1.e-9      
+      ELSE
+         IF (lwp) PRINT *, 'cons heat : ', kt, zdiff_hc1 * 1.e-20 * rau0 * rcp, zdiff_sc1 * 1.e-9
+         IF (lwp) PRINT *, 'cons vol  : ', kt, zdiff_v1 * 1.e-9      
+      END IF
       IF( lk_vvl ) THEN
         CALL iom_put( 'bgtemper' , zdiff_hc / zvol_tot )              ! Temperature variation (C)