Changeset 1944 for trunk/LMDZ.MARS/libf/phymars/dyn1d

Timestamp:

Jun 5, 2018, 8:38:22 AM (7 years ago)

Author:

emillour

Message:

Mars GCM:
A first step towards 1+1=2 (for now only works without tracers):

• store and load "albedo" (surface albedo) and wstar (thermals' max vertical velocity) in physics (re)start file.
• turn phyetat0 into a module in the process.

EM

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F (modified) (5 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F

@@ -21 +21 @@
      &                            ecritphy, iphysiq
       use dimradmars_mod, only: tauscaling,tauvis,totcloudfrac
+      use co2cloud_mod, only: mem_Mccn_co2, mem_Mh2o_co2,
+     &                        mem_Nccn_co2
       USE comvert_mod, ONLY: ap,bp,aps,bps,pa,preff,sig,
      &                       presnivs,pseudoalt,scaleheight
@@ -92 +94 @@
       REAL co2ice(1)        ! co2ice layer (kg.m-2)
       REAL emis(1)          ! surface layer
+      REAL albedo(1,1)      ! surface albedo
+      REAL :: wstar(1)=0.    ! Thermals vertical velocity
       REAL q2(nlayer+1)   ! Turbulent Kinetic Energy
       REAL zlay(nlayer)   ! altitude estimee dans les couches (km)
@@ -522 +526 @@
       call getin("albedo",albedodat(1))
       write(*,*) " albedo = ",albedodat(1)
+      albedo(1,1)=albedodat(1)
 
       inertiedat(1,1)=400 ! default value for inertiedat
@@ -593 +598 @@
       call getin("co2ice",co2ice)
       write(*,*) " co2ice = ",co2ice
-
+! Initialization for CO2 clouds (could be improved to read initial profiles)
+      mem_Mccn_co2(:,:)=0
+      mem_Mh2o_co2(:,:)=0
+      mem_Nccn_co2(:,:)=0
 c
 c  emissivity
@@ -710 +718 @@
      &              albedodat,inertiedat,zmea,zstd,zsig,zgam,zthe)
       call physdem1("startfi.nc",nsoilmx,ngrid,llm,nq,
-     .              dtphys,time,
-     .              tsurf,tsoil,co2ice,emis,q2,qsurf,tauscaling,
-     .              totcloudfrac)
+     &              dtphys,time,
+     &              tsurf,tsoil,co2ice,albedo,emis,q2,qsurf,tauscaling,
+     &              totcloudfrac,wstar,
+     &              mem_Mccn_co2,mem_Nccn_co2,
+     &              mem_Mh2o_co2)
 
 c=======================================================================