Changeset 2909 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Mar 9, 2023, 9:38:24 AM (22 months ago)

Author:

romain.vande

Message:

Mars PEM:
New Boolean options for following orbital parameters of ob_ex_lsp.asc: var_obl, var_ex, var_lsp.
If using evol_orbit_pem=true, you can specify which parameter to follow.
True by default: Do you want to change the parameter XXX during the PEM run as prescribed in ob_ex_lsp.asc.
If false, it is set to constant (to the value of the tab_cntrl in the start)
RV

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• comslope_mod.F90 (modified) (3 diffs)
• comsoil_h.F90 (modified) (1 diff)
• dimradmars_mod.F90 (modified) (1 diff)
• phyetat0_mod.F90 (modified) (3 diffs)
• phys_state_var_init_mod.F90 (modified) (2 diffs)
• physiq_mod.F (modified) (1 diff)
• surfdat_h.F90 (modified) (1 diff)

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

@@ -23 +23 @@
 
 CONTAINS 
-  SUBROUTINE ini_comslope_h(ngrid)
+  SUBROUTINE ini_comslope_h(ngrid,nslope_in)
 
   IMPLICIT NONE
   INTEGER, INTENT(IN) :: ngrid
+  INTEGER, INTENT(IN) :: nslope_in
 
-    allocate(def_slope(nslope+1))
-    allocate(def_slope_mean(nslope))
-    allocate(sky_slope(nslope))
-    allocate(subslope_dist(ngrid,nslope))
+    allocate(def_slope(nslope_in+1))
+    allocate(def_slope_mean(nslope_in))
+    allocate(sky_slope(nslope_in))
+    allocate(subslope_dist(ngrid,nslope_in))
     allocate(major_slope(ngrid))
   END SUBROUTINE ini_comslope_h
@@ -70 +71 @@
       qsurf_meshavg(:,:) = 0.
 
+  if(nslope.eq.1) then
+      albedo_meshavg(:,:) = albedo_slope(:,:,1)
+      emis_meshavg(:) = emis_slope(:,1)
+      tsurf_meshavg(:) = tsurf_slope(:,1)
+      qsurf_meshavg(:,:) = qsurf_slope(:,:,1)
+  else
+
   DO ig = 1,ngrid
       DO islope = 1,nslope
@@ -84 +92 @@
   ENDDO
 
+  endif
+
   END SUBROUTINE compute_meshgridavg
 

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

@@ -70 +70 @@
   end subroutine end_comsoil_h
 
+  subroutine ini_comsoil_h_slope_var(ngrid,nslope)
+  
+  implicit none
+  integer,intent(in) :: ngrid ! number of atmospheric columns
+  integer,intent(in) :: nslope ! number of sub grid slopes
+ 
+    allocate(tsoil(ngrid,nsoilmx,nslope)) ! soil temperatures
+    allocate(mthermdiff(ngrid,0:nsoilmx-1,nslope))
+    allocate(thermdiff(ngrid,nsoilmx-1,nslope))
+    allocate(coefd(ngrid,nsoilmx-1,nslope))
+    allocate(alph(ngrid,nsoilmx-1,nslope))
+    allocate(beta(ngrid,nsoilmx-1,nslope))
+ 
+  end subroutine ini_comsoil_h_slope_var
+
+
+  subroutine end_comsoil_h_slope_var
+
+  implicit none
+
+    if (allocated(tsoil)) deallocate(tsoil)
+    if (allocated(mthermdiff)) deallocate(mthermdiff)
+    if (allocated(thermdiff)) deallocate(thermdiff)
+    if (allocated(coefd)) deallocate(coefd)
+    if (allocated(alph)) deallocate(alph)
+    if (allocated(beta)) deallocate(beta)
+
+  end subroutine end_comsoil_h_slope_var
+
 end module comsoil_h

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

@@ -216 +216 @@
 
   end subroutine end_dimradmars_mod
+
+  subroutine ini_dimradmars_mod_slope_var(ngrid,nslope)
+  
+  implicit none
+  
+  integer,intent(in) :: ngrid ! number of atmospheric columns
+  integer,intent(in) :: nslope ! number of subgrid scale slopes
+
+   allocate(albedo(ngrid,2,nslope))
+   allocate(fluxrad_sky(ngrid,nslope))
+   allocate(fluxrad(ngrid,nslope))
+
+  end subroutine ini_dimradmars_mod_slope_var
+
+  subroutine end_dimradmars_mod_slope_var
+
+  implicit none
+
+   if (allocated(albedo))      deallocate(albedo)
+   if (allocated(fluxrad_sky)) deallocate(fluxrad_sky)
+   if (allocated(fluxrad))     deallocate(fluxrad)
+
+  end subroutine end_dimradmars_mod_slope_var
 
  

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

@@ -134 +134 @@
 endif ! of if (startphy_file)
 
-if(nslope.ne.1) then
-  call abort_physic(modname, &
-                "phyetat0: For now, nslope should be 1 (set in comslope_mod)",1)
-endif
-
-allocate(default_def_slope(nslope+1))
-!Sub-grid scale  subslopes
-if (nslope.eq.7) then
-  default_def_slope(1) = -43.
-  default_def_slope(2) = -19.
-  default_def_slope(3) = -9.
-  default_def_slope(4) = -3.
-  default_def_slope(5) = 3.
-  default_def_slope(6) = 9.
-  default_def_slope(7) = 19.
-  default_def_slope(8) = 43.
-elseif (nslope.eq.5) then
-  default_def_slope(1) = -43.
-  default_def_slope(2) = -9.
-  default_def_slope(3) = -3.
-  default_def_slope(4) = 3.
-  default_def_slope(5) = 9.
-  default_def_slope(6) = 43.
-elseif (nslope.eq.1) then
-  default_def_slope(1) = 0.
-  default_def_slope(2) = 0.
-endif
-
 if (startphy_file) then
  call get_var("def_slope",def_slope,found)
@@ -167 +139 @@
      startphy_slope=.false.
      write(*,*)'slope_settings: Problem while reading <def_slope>'
-     write(*,*)'default def_slope will be used'
-     do islope=1,nslope+1
-       def_slope(islope) = default_def_slope(islope)
-     enddo
-     write(*,*)'computing corresponding distribution <subslope_dist>'
-     write(*,*)'For now, woth nslope=1, subslope_dist is straigforward'
-     write(*,*)'Later this operation will be done by newstart with a specific routine'
-     subslope_dist(:,:)=1.
-     !call subslope_mola(ngrid,nslope,def_slope,subslope_dist)
+     write(*,*)'Checking that nslope=1'
+     if(nslope.eq.1) then
+      write(*,*)'We take default def_slope and subslope dist'
+      allocate(default_def_slope(nslope+1))
+      default_def_slope(1) = 0.
+      default_def_slope(2) = 0.
+      subslope_dist(:,:)=1.
+     else
+       write(*,*)'Variable def_slope is not present in the start and'
+       write(*,*)'you are trying to run with nslope!=1'
+       write(*,*)'This is not possible, you have to go through newstart'
+     endif
    else
      startphy_slope=.true.
@@ -181 +156 @@
      if(.not.found) then 
        write(*,*)'slope_settings: Problem while reading <subslope_dist>'
-       write(*,*)'computing a new distribution'
-       write(*,*)'For now, woth nslope=1, subslope_dist is straigforward'
-       write(*,*)'Later this operation will be done by newstart with a specific routine'
-       subslope_dist(:,:)=1.
-     !call subslope_mola(ngrid,nslope,def_slope,subslope_dist)
+       write(*,*)'We have to abort.'
+       write(*,*)'Please check that nslope is coherent, you can modify this parameter with newstart'
+       call abort_physic(modname, &
+                "phyetat0: Failed loading <subslope_dist>",1)
      endif
     endif
-else ! startphy_file  
- do islope=1,nslope+1
-       def_slope(islope) = default_def_slope(islope)
-     enddo
-     write(*,*)'computing corresponding distribution <subslope_dist>'
-     write(*,*)'For now, woth nslope=1, subslope_dist is straigforward'
-     write(*,*)'Later this operation will be done by newstart with a specific routine'
-     subslope_dist(:,:)=1.
-     !call subslope_mola(ngrid,nslope,def_slope,subslope_dist)
+else ! (no startphy_file)
+     if(nslope.eq.1) then
+      allocate(default_def_slope(nslope+1))
+      default_def_slope(1) = 0.
+      default_def_slope(2) = 0.
+      subslope_dist(:,:)=1.
+    else
+     write(*,*)'Without starfi, lets first run with nslope=1'
+       call abort_physic(modname, &
+                "phyetat0: No startfi and nslope!=1",1)
+   endif
 endif   
 

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

@@ -61 +61 @@
       use dust_rad_adjust_mod, only: ini_dust_rad_adjust_mod, &
                                      end_dust_rad_adjust_mod
-      use comslope_mod, ONLY: nslope
+      use comslope_mod, ONLY: nslope,end_comslope_h,ini_comslope_h
       use netcdf
       USE mod_phys_lmdz_para, ONLY: is_master,bcast
@@ -182 +182 @@
       call ini_dust_rad_adjust_mod(ngrid)
 
+      ! allocate arrays in "comslope_mod"
+      call end_comslope_h
+      call ini_comslope_h(ngrid,nslope)
+
       END SUBROUTINE phys_state_var_init
 

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -577 +577 @@
 c        ~~~~~~~~~~~~
 #ifndef MESOSCALE
-
-         call ini_comslope_h(ngrid)
 
 ! GCM. Read netcdf initial physical parameters.

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

@@ -116 +116 @@
   end subroutine end_surfdat_h
 
+  subroutine ini_surfdat_h_slope_var(ngrid,nq,nslope)
+  
+  implicit none
+  integer,intent(in) :: ngrid ! number of atmospheric columns
+  integer,intent(in) :: nq ! number of tracers  
+  integer,intent(in) :: nslope ! number of sub-grid scale slope  
+    allocate(qsurf(ngrid,nq,nslope))
+    allocate(tsurf(ngrid,nslope))
+    allocate(watercap(ngrid,nslope))
+    allocate(emis(ngrid,nslope))
+    allocate(capcal(ngrid,nslope))
+    allocate(fluxgrd(ngrid,nslope))
+       
+  end subroutine ini_surfdat_h_slope_var
+
+  subroutine end_surfdat_h_slope_var
+
+  implicit none
+
+    if (allocated(qsurf))         deallocate(qsurf)
+    if (allocated(tsurf))         deallocate(tsurf)
+    if (allocated(watercap))      deallocate(watercap)
+    if (allocated(emis))          deallocate(emis)
+    if (allocated(capcal))        deallocate(capcal)
+    if (allocated(fluxgrd))       deallocate(fluxgrd)
+    
+  end subroutine end_surfdat_h_slope_var
+
 end module surfdat_h