Changeset 1711 for trunk/LMDZ.MARS/libf/phymars/phyetat0.F90

Timestamp:

May 19, 2017, 11:19:17 AM (8 years ago)

Author:

mvals

Message:

Mars GCM

• implementation of a sub-grid water cloud fraction scheme (by A. Pottier): If CLFvarying is set to true in callphys.def, then the sub-grid cloud fraction routine is applied in watercloud.F and aeropacity.F.
• accordingly modified files are: watercloud.F, aeropacity.F, callradite.F, conf_phys.F, phyetat0.F, phyredem.F90, physiq_mod.F, testphys1d.F, callkeys.h, newtsart.F, start2archive.F, lect_start_archive.F
• added file: tcondwater.F90, used by watercloud.F to calculate the condensation temperature of water
• watercloud.F, aeropacity.F, callradite.F are converted to module files as watercloud_mod.F, aeropacity_mod.F, callradite_mod.F

MV

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0.F90 (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/phyetat0.F90

@@ -1 +1 @@
 subroutine phyetat0 (fichnom,tab0,Lmodif,nsoil,ngrid,nlay,nq, &
                      day_ini,time0,tsurf,tsoil,emis,q2,qsurf,co2ice, &
-                     tauscaling)
+                     tauscaling,totcloudfrac)
 !  use netcdf
   use tracer_mod, only: noms ! tracer names
@@ -46 +46 @@
   real,intent(out) :: co2ice(ngrid) ! co2 ice cover
   real,intent(out) :: tauscaling(ngrid) ! dust conversion factor
-
+  real,intent(out) :: totcloudfrac(ngrid) ! total cloud fraction
 !======================================================================
 !  Local variables:
@@ -241 +241 @@
 endif
 
+! Sub-grid cloud fraction
+call get_field("totcloudfrac",totcloudfrac,found,indextime)
+if (.not.found) then
+  write(*,*) "phyetat0: <totcloudfrac> not in file WARNING put to 1"
+  totcloudfrac(:) = 1.0 !valeur par defaut (CLFfixval par defaut)
+else
+  write(*,*) "phyetat0: total cloud fraction <totcloudfrac> range:", &
+             minval(totcloudfrac), maxval(totcloudfrac)
+endif
 
 ! Surface temperature :