Changeset 1493

Timestamp:

Nov 4, 2015, 8:58:27 PM (9 years ago)

Author:

mturbet

Message:

some corrections

Location:

trunk/LMDZ.GENERIC

Files:

• README (modified) (1 diff)
• libf/phystd/phyetat0.F90 (modified) (1 diff)
• libf/phystd/physiq.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/README

@@ -1100 +1100 @@
 == 22/10/2015 == MT
 - Cleaning of condense_co2.F90.
+
+== 04/11/2015 == MT
+- Corrections of bugs in phyetat0.F90 (vap tracer) and in physiq.F90 (ice evolution scheme).

trunk/LMDZ.GENERIC/libf/phystd/phyetat0.F90

@@ -376 +376 @@
   do iq=1,nq
     txt=tname(iq)
-    if (txt.eq."h2o_vap") then
+    
+    !! There was a bug here. MT2015.
+    
+    !if (txt.eq."h2o_vap") then
       ! There is no surface tracer for h2o_vap;
       ! "h2o_ice" should be loaded instead
-      txt="h2o_ice"
-      write(*,*) 'phyetat0: loading surface tracer', &
-                           ' h2o_ice instead of h2o_vap'
-    endif
+     ! txt="h2o_ice"
+     ! write(*,*) 'phyetat0: loading surface tracer', &
+     !                      ' h2o_ice instead of h2o_vap'
+    !endif
+    
     call get_field(txt,qsurf(:,iq),found,indextime)
     if (.not.found) then

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -1831 +1831 @@
 
             ! enforce ice conservation
-            ice_tot= SUM(qsurf_hist(:,igcm_h2o_ice)*area(:) )
+            ice_tot= SUM(qsurf_hist(:,igcm_h2o_ice)*area(:) )/SUM(area(:))
             qsurf_hist(:,igcm_h2o_ice) = qsurf_hist(:,igcm_h2o_ice)*(icesrf/ice_tot)