Changeset 146 for branches/iLoveclim/SOURCES/Hemin40_files/output_hemin40_mod.f90

Timestamp:

10/20/17 09:31:39 (7 years ago)

Author:

aquiquet

Message:

Grisli-iLoveclim branch: merged to trunk at revision 145

Location:

branches/iLoveclim

Files:

• . (modified) (1 prop)
• SOURCES/Hemin40_files/output_hemin40_mod.f90 (modified) (3 diffs)

branches/iLoveclim/SOURCES/Hemin40_files/output_hemin40_mod.f90

@@ -30 +30 @@
 ! afq for CONSEAU REAL, dimension(nx,ny) :: old_H_dtt      ! Epaisseur de glace au pas de temps precedent
 integer, dimension(nx,ny) :: write_mask
+real :: ilc_units                     !afq, unit change from Grisli to iLoveclim
 
 CONTAINS
@@ -184 +185 @@
 ! afq for CONSEAU      REAL, dimension(nx,ny) :: delta_H_dtt
 
+      ilc_units = dx*dy/DICE
 
 !     open(unit=4145,file='reg_output_nord.dat')
@@ -330 +332 @@
 
 ! afq -- iLOVECLIM water conservation: 
-      trendWAC        = diff_H*dx*dy !or water_bilan
-      smbWAC(:,:)     = (ablbord_dtt(:,:)+ Bm_dtt(:,:))*dx*dy/dtt
-      bmeltWAC(:,:)   = Bmelt_dtt(:,:)*dx*dy/dtt
-      calvingWAC(:,:) = Calv_dtt(:,:)*dx*dy/dtt
+!        units: GRISLI is in m i.e. / yr
+!               iLOVECLIM expects m3 w.e. /yr
+      trendWAC        = diff_H*ilc_units !or water_bilan
+      smbWAC(:,:)     = (-ablbord_dtt(:,:)+ Bm_dtt(:,:))*ilc_units/dtt
+      bmeltWAC(:,:)   = -Bmelt_dtt(:,:)*ilc_units/dtt
+      calvingWAC(:,:) = Calv_dtt(:,:)*ilc_units/dtt
 !check->
       if (abs(1-diff_H/(water_bilan+1e-20)).gt.0.01) then
-         write(*,*) "Water is not conserved in GRISLI", diff_H, water_bilan, abs(1-diff_H/(water_bilan+1e-20))
+         write(*,*) "Water not conserved in GRISLI", diff_H, water_bilan, abs(1-diff_H/(water_bilan+1e-20))
       endif