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

Timestamp:

Jan 25, 2024, 6:06:27 PM (2 years ago)

Author:

jbclement

Message:

Mars PCM:
Some changes to prepare the introduction of slopes in 1D: in particular, the subroutine "getslopes.F90" and "param_slope.F90" are now inside the module "slope_mod.F90" + Few small cleanings throughout the code.
JBC

Location:

trunk/LMDZ.MARS/libf/phymars/dyn1d

Files:

• init_testphys1d_mod.F90 (modified) (6 diffs)
• testphys1d.F90 (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -41 +41 @@
 use conf_phys_mod,            only: conf_phys
 use paleoclimate_mod,         only: paleoclimate
+use comslope_mod,             only: nslope, subslope_dist, ini_comslope_h, end_comslope_h
 ! Mostly for XIOS outputs:
 use mod_const_mpi,            only: COMM_LMDZ
@@ -208 +209 @@
 allocate(nqfils(nqtot))
 nqperes = 0
-nqfils(:) = 0
+nqfils = 0
 do iq = 1,nqtot
     if (tnom_transp(iq) == 'air') then
@@ -373 +374 @@
 longitude = longitude*pi/180.
 
-! Some initializations (some of which have already been
-! done above!) and loads parameters set in callphys.def
-! and allocates some arrays
+! Some initializations (some of which have already been done above!)
+! and loads parameters set in callphys.def and allocates some arrays
 ! Mars possible matter with dttestphys in input and include!!!
 ! Initializations below should mimick what is done in iniphysiq for 3D GCM
@@ -383 +383 @@
 call init_geometry_cell_area_for_outputs(1,cell_area)
 ! Ehouarn: init_vertial_layers called later (because disvert not called yet)
-!      call init_vertical_layers(nlayer,preff,scaleheight,
-!     &                      ap,bp,aps,bps,presnivs,pseudoalt)
+! call init_vertical_layers(nlayer,preff,scaleheight,ap,bp,aps,bps,presnivs,pseudoalt)
 call init_dimphy(1,nlayer) ! Initialize dimphy module
-call phys_state_var_init(1,llm,nq,tname,day0,dayn,time,daysec,dttestphys,rad,g,r,cpp,nqperes,nqfils)! MVals: variables isotopes
+call phys_state_var_init(1,llm,nq,tname,day0,dayn,time,daysec,dttestphys,rad,g,r,cpp,nqperes,nqfils) ! MVals: variables isotopes
 call ini_fillgeom(1,latitude,longitude,(/1.0/))
 call conf_phys(1,llm,nq)
@@ -461 +460 @@
 
 ! For the non-orographic gravity waves scheme
-du_nonoro_gwd(1,:)=0
-dv_nonoro_gwd(1,:)=0
+du_nonoro_gwd(1,:) = 0
+dv_nonoro_gwd(1,:) = 0
 
 ! For the slope wind scheme
@@ -559 +558 @@
 if (.not. therestart1D) then
     tsurf(:,1) = tmp2(0)
-    temp(:) = tmp2(1:)
+    temp = tmp2(1:)
 else
     read(3,*) header, (tsurf(1,:), j = 1,size(tsurf,2)), (temp(ilayer), ilayer = 1,nlayer)

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

@@ -255 +255 @@
     ! Compute winds and temperature for next time step
     ! ------------------------------------------------
-    u(:) = u(:) + dttestphys*du(:)
-    v(:) = v(:) + dttestphys*dv(:)
-    temp(:) = temp(:) + dttestphys*dtemp(:)
+    u = u + dttestphys*du
+    v = v + dttestphys*dv
+    temp = temp + dttestphys*dtemp
 
     ! Compute pressure for next time step