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

Timestamp:

Jan 22, 2019, 5:02:23 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:

• Update of rocketduststorm_mod.F90 :

We want to separate both parametrizations related to the formation of the detached dust layers. Therefore, rocketduststorm_mod.F90 now only comprises the rocket dust storm scheme, whereas it contained
also before the calculation of the vertical velocity induced by the presence of the sub-grid scale topography. This latter part is under development and will be integrated as a fully independant
parametrization: the aim is to simulate the entrainment of dust by slope winds, from the boundary layer up to the top of the sub-grid scale topography.

• Addition of initial surface parameters "summit" and "base" to prepare the previously described slope wind parametrization

MV

Location:

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Files:

• datareadnc.F (modified) (5 diffs)
• newstart.F (modified) (5 diffs)

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

@@ -2 +2 @@
       SUBROUTINE datareadnc(relief,phisinit,alb,ith,z0,
      &                    zmea,zstd,zsig,zgam,zthe,
-     &                    hmons,summit,zavg)
+     &                    hmons,summit,base,zavg)
 c=======================================================================
 c
@@ -79 +79 @@
       REAL,intent(out) :: zthe(imdp1*jmdp1)
       REAL,intent(out) :: hmons(imdp1*jmdp1) !CW17,361*180 hmons
-      REAL,intent(out) :: summit(imdp1*jmdp1) !CW17,361*180 hmons
-      REAL,intent(out) :: zavg(imdp1*jmdp1) !CW17,361*180 hmons
+      REAL,intent(out) :: summit(imdp1*jmdp1)
+      REAL,intent(out) :: base(imdp1*jmdp1)
+      REAL,intent(out) :: zavg(imdp1*jmdp1)
       
 ! Local variables:
@@ -327 +328 @@
 c     we don't need to put them into restartfi.nc
       string(6) = 'summit' !subgrid summit
-c     string(7) = 'base'   !base of summit (not subgrid)
+      string(7) = 'base'   !base of summit (not subgrid)
 c     string(8) = 'bottom' !subgrid base
-      do k=5,6
+      do k=5,7
           write(*,*) 'string',k,string(k)  
 c-----------------------------------------------------------------------
@@ -373 +374 @@
           enddo
 
-          if (k .eq. 5 .or. k .eq. 6) then 
+          if (k .eq. 5 .or. k .eq. 6 .or. k .eq. 7) then 
               ! hmons, summit, keep the maximum of subgrids 
               call mvc_horiz(imd,jmdp1,iim,jjm,longitude,
@@ -404 +405 @@
           if (k.eq.6) then  ! summit
              do i=1,iimp1*jjp1
+                  if (pfield(i) .ne. -999999.) then             
                     summit(i) = pfield(i)
+                  else
+                    summit(i)=0.
+                  endif
+             enddo
+          endif
+
+          if (k.eq.7) then  ! base
+             do i=1,iimp1*jjp1
+                  if (pfield(i) .ne. -999999.) then             
+                    base(i) = pfield(i)
+                  else
+                    base(i)=0.
+                  endif
              enddo
           endif

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

@@ -27 +27 @@
       use surfdat_h, only: phisfi, z0, zmea, zstd, zsig, zgam, zthe,
      &                     albedodat, z0_default, qsurf, tsurf,
-     &                     co2ice, emis, hmons
+     &                     co2ice, emis, hmons, summit, base
       use comsoil_h, only: inertiedat, layer, mlayer, nsoilmx, tsoil
       use control_mod, only: day_step, iphysiq, anneeref, planet_type
@@ -100 +100 @@
       real hmonsS(iip1,jjp1)
       real summitS(iip1,jjp1)
+      real baseS(iip1,jjp1)
       real zavgS(iip1,jjp1)
       real z0S(iip1,jjp1)
@@ -385 +386 @@
       CALL datareadnc(relief,phis,alb,surfith,z0S,
      &          zmeaS,zstdS,zsigS,zgamS,ztheS,
-     &          hmonsS,summitS,zavgS)
+     &          hmonsS,summitS,baseS,zavgS)
 
       CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phis,phisfi)
@@ -398 +399 @@
       CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,ztheS,zthe)
       CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,hmonsS,hmons)
-!      CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,summitS,summit)
+      CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,summitS,summit)
+      CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,baseS,base)
 !      CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,zavgS,zavg)
 
@@ -1651 +1653 @@
      &              nqtot,dtphys,real(day_ini),0.0,cell_area,
      &              albfi,ithfi,zmea,zstd,zsig,zgam,zthe,
-     &              hmons)
+     &              hmons,summit,base)
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
      &              dtphys,hour_ini,