Changeset 4043 for trunk/LMDZ.MARS/libf

Timestamp:

Feb 5, 2026, 2:47:26 PM (4 weeks ago)

Author:

jbclement

Message:

Mars PCM:
Addition of 'pa' and 'preff' in the information provided by "start1D.txt".
JBC

Location:

trunk/LMDZ.MARS/libf/phymars/dyn1d

Files:

• init_testphys1d_mod.F90 (modified) (3 diffs)
• testphys1d.F90 (modified) (2 diffs)
• writerestart1D_mod.F90 (modified) (2 diffs)

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

@@ -313 +313 @@
 ! ------------------------
 psurf = 610. ! Default value for psurf
+pa = 20.     ! Default value for transition pressure (for hybrid coord.)
+preff = 610. ! Default value for reference surface pressure
 write(*,*) 'Surface pressure (Pa)?'
 if (.not. therestart1D) then
@@ -318 +320 @@
 else
     open(3,file = trim(start1Dname),status = "old",action = "read")
-    read(3,*) header, psurf
-endif
-write(*,*) " psurf = ",psurf
+    read(3,*) header, psurf, pa, preff
+endif
+write(*,*) " psurf = ", psurf
+write(*,*) " pa    = ", pa
+write(*,*) " preff = ", preff
 
 ! Compute pressures and altitudes of atmospheric levels
@@ -330 +334 @@
 call getin("hybrid",hybrid)
 write(*,*) " hybrid = ", hybrid
-
-! Reference pressures
-pa = 20.     ! Transition pressure (for hybrid coord.)
-preff = 610. ! Reference surface pressure
 
 call disvert_noterre

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

@@ -11 +11 @@
 use dimradmars_mod,      only: tauvis, totcloudfrac, albedo
 use dust_param_mod,      only: tauscaling
-use comvert_mod,         only: ap, bp, aps, bps
+use comvert_mod,         only: ap, bp, aps, bps, pa, preff
 use physiq_mod,          only: physiq
 use turb_mod,            only: q2
@@ -261 +261 @@
 
 ! Writing the "restart1D.txt" file for the next run
-if (startfiles_1D) call writerestart1D('restart1D.txt',psurf,tsurf(1,:),nlayer,size(tsurf,2),temp,u,v,nq,noms,qsurf(1,:,:),q)
+if (startfiles_1D) call writerestart1D('restart1D.txt',psurf,pa,preff,tsurf(1,:),nlayer,size(tsurf,2),temp,u,v,nq,noms,qsurf(1,:,:),q)
 
 write(*,*) "testphys1d: everything is cool!"

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

@@ -5 +5 @@
 contains
 
-SUBROUTINE writerestart1D(filename,psurf,tsurf,nlayer,nslope,temp,u,v,nq,qnames,qsurf,q)
+SUBROUTINE writerestart1D(filename,psurf,pa,preff,tsurf,nlayer,nslope,temp,u,v,nq,qnames,qsurf,q)
 
 implicit none
 
 ! Arguments
-character(len = *),                intent(in) :: filename
-integer,                           intent(in) :: nlayer, nq, nslope
-real,                              intent(in) :: psurf
-real, dimension(nslope),           intent(in) :: tsurf
-real, dimension(nlayer),           intent(in) :: temp, u, v
-real, dimension(nlayer,nq),        intent(in) :: q
-real, dimension(nq,nslope),        intent(in) :: qsurf
-character(len = *), dimension(nq), intent(in) :: qnames
+character(*),                intent(in) :: filename
+integer,                     intent(in) :: nlayer, nq, nslope
+real,                        intent(in) :: psurf, pa, preff
+real, dimension(nslope),     intent(in) :: tsurf
+real, dimension(nlayer),     intent(in) :: temp, u, v
+real, dimension(nlayer,nq),  intent(in) :: q
+real, dimension(nq,nslope),  intent(in) :: qsurf
+character(*), dimension(nq), intent(in) :: qnames
 
 ! Local variables
@@ -24 +24 @@
 ! Write the data needed for a restart in "restart1D.txt"
 open(1,file = filename,status = "replace",action = "write")
-write(1,*) 'ps', psurf
+write(1,*) 'ps', psurf, pa, preff
 do i = 1,nq
     write(1,*) qnames(i), (qsurf(i,j), j = 1,nslope), (q(il,i), il = 1,nlayer)