Changeset 1711 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

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

Location:

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Files:

• lect_start_archive.F (modified) (9 diffs)
• newstart.F (modified) (4 diffs)
• start2archive.F (modified) (6 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F

@@ -2 +2 @@
      &     date,tsurf,tsoil,emis,q2,
      &     t,ucov,vcov,ps,co2ice,h,phisold_newgrid,
-     &     q,qsurf,tauscaling,surfith,nid)
+     &     q,qsurf,tauscaling,totcloudfrac,surfith,nid)
 c=======================================================================
 c
@@ -49 +49 @@
       REAL,INTENT(OUT) :: q2(ngrid,nlayer+1),qsurf(ngrid,nqtot)
       REAL,INTENT(OUT) :: tauscaling(ngrid) ! dust conversion factor
+      REAL,INTENT(OUT) :: totcloudfrac(ngrid) ! sub grid cloud fraction
       REAL,INTENT(OUT) :: phisold_newgrid(iip1,jjp1)
       REAL,INTENT(OUT) :: t(iip1,jjp1,llm)
@@ -130 +131 @@
       REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqtot)
       real tauscalingS(iip1,jjp1)
+      real totcloudfracS(iip1,jjp1)
 
       real ptotal, co2icetotal
@@ -158 +160 @@
       real, dimension(:,:,:,:), allocatable :: qold
       real, dimension(:,:), allocatable :: tauscalingold
+      real, dimension(:,:), allocatable :: totcloudfracold
 
       real tab_cntrl(100)
@@ -358 +361 @@
       allocate(qsurfold(imold+1,jmold+1,nqtot))
       allocate(tauscalingold(imold+1,jmold+1))
+      allocate(totcloudfracold(imold+1,jmold+1))
 
       allocate(var (imold+1,jmold+1,llm))
@@ -683 +687 @@
 
 c-----------------------------------------------------------------------
-c 5.1 Lecture des champs 2D (co2ice, emis,ps,tsurf,Tg[10], q2surf, tauscaling)
+c 5.1 Lecture des champs 2D (co2ice, emis,ps,tsurf,Tg[10], q2surf, tauscaling,totcloudfrac)
 c-----------------------------------------------------------------------
  
@@ -777 +781 @@
         IF (ierr .NE. NF_NOERR) THEN
            PRINT*, "lect_start_archive: Failed loading <tauscaling>"
+           PRINT*, NF_STRERROR(ierr)
+           CALL abort
+        ENDIF
+      ENDIF
+
+      ierr = NF_INQ_VARID (nid, "totcloudfrac", nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: <totcloudfrac> not in file"
+         totcloudfracold(:,:) = 1
+      ELSE
+#ifdef NC_DOUBLE
+        ierr = NF_GET_VARA_DOUBLE(nid, nvarid,start,count,
+     &                            totcloudfracold)
+#else
+        ierr = NF_GET_VARA_REAL(nid, nvarid,start,count,
+     &                          totcloudfracold)
+#endif
+        IF (ierr .NE. NF_NOERR) THEN
+           PRINT*, "lect_start_archive: Failed loading <totcloudfrac>"
            PRINT*, NF_STRERROR(ierr)
            CALL abort
@@ -1063 +1086 @@
      &                   1,rlonuold,rlatvold,rlonu,rlatv)
       call gr_dyn_fi (1,iim+1,jjm+1,ngrid,tauscalings,tauscaling)
+      
+c Sub grid cloud fraction
+      call interp_horiz (totcloudfracold,totcloudfracs,imold,jmold,iim,
+     &                   jjm,1,rlonuold,rlatvold,rlonu,rlatv)
+      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,totcloudfracs,totcloudfrac)
+      
 c     write(46,*) 'emis',emis
 c-----------------------------------------------------------------------
@@ -1422 +1451 @@
       deallocate(qsurfold)
       deallocate(tauscalingold)
+      deallocate(totcloudfracold)
       deallocate(var,varp1)
 

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -170 +170 @@
       real,allocatable :: coefvmr(:)  ! Correction coefficient when changing composition
       integer :: iloc(1), iqmax
+! sub-grid cloud fraction
+      real :: totcloudfrac(ngridmx)
 
 c sortie visu pour les champs dynamiques
@@ -402 +404 @@
      &   date,tsurf,tsoil,emis,q2,
      &   t,ucov,vcov,ps,co2ice,teta,phisold_newgrid,q,qsurf,
-     &   tauscaling,surfith,nid)
+     &   tauscaling,totcloudfrac,surfith,nid)
         write(*,*) "OK, read start_archive file"
         ! copy soil thermal inertia
@@ -423 +425 @@
         CALL phyetat0 (fichnom,tab0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
      .        day_ini,time,
-     .        tsurf,tsoil,emis,q2,qsurf,co2ice,tauscaling)
+     .        tsurf,tsoil,emis,q2,qsurf,co2ice,tauscaling,totcloudfrac)
         
         ! copy albedo and soil thermal inertia
@@ -1640 +1642 @@
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
      .              dtphys,hour_ini,
-     .              tsurf,tsoil,co2ice,emis,q2,qsurf,tauscaling)
+     .              tsurf,tsoil,co2ice,emis,q2,qsurf,tauscaling,
+     .              totcloudfrac)
 
 c=======================================================================

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

@@ -64 +64 @@
       REAL co2ice(ngridmx)        ! CO2 ice layer
       REAL tauscaling(ngridmx) ! dust conversion factor
+      REAL totcloudfrac(ngridmx) ! sub-grid cloud fraction
       REAL q2(ngridmx,llm+1)
       REAL emis(ngridmx)
@@ -79 +80 @@
       REAL co2iceS(ip1jmp1)
       REAL tauscalingS(ip1jmp1)
+      REAL totcloudfracS(ip1jmp1)
       REAL q2S(ip1jmp1,llm+1)
       REAL,ALLOCATABLE :: qsurfS(:,:)
@@ -156 +158 @@
       CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
      &      day_ini_fi,timefi,tsurf,tsoil,emis,q2,qsurf,co2ice,
-     &      tauscaling)
+     &      tauscaling,totcloudfrac)
 
        ierr = NF_OPEN (fichnom, NF_NOWRITE,nid1)
@@ -244 +246 @@
 c qsurf --> qsurfS
 c tauscaling --> tauscalingS
+c totcloudfrac --> totcloudfracS
 c
 c-----------------------------------------------------------------------
@@ -256 +259 @@
       call gr_fi_dyn(nqtot,ngridmx,iip1,jjp1,qsurf,qsurfS)
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,tauscaling,tauscalingS)
+      call gr_fi_dyn(1,ngridmx,iip1,jjp1,totcloudfrac,totcloudfracS)
 
 c=======================================================================
@@ -346 +350 @@
       call write_archive(nid,ntime,'tauscaling',
      &  'dust conversion factor',' ',2,tauscalingS)
+      call write_archive(nid,ntime,'totcloudfrac',
+     &  'sub grid cloud fraction',' ',2,totcloudfracS)
       call write_archive(nid,ntime,'emis','grd emis',' ',2,emisS)
       call write_archive(nid,ntime,'ps','Psurf','Pa',2,ps)