Changeset 4059 for trunk/LMDZ.MARS/libf

Timestamp:

Feb 10, 2026, 10:18:01 AM (2 weeks ago)

Author:

jbclement

Message:

Mars PCM:
Deletion of 'qsurf' in the information provided by "start1D.txt" since it is unused and unnecessary until now. Addition of an error-leading check to prevent an unintended situation (which should not happen) regarding the presence of starting files.
JBC

Location:

trunk/LMDZ.MARS/libf/phymars/dyn1d

Files:

• init_testphys1d_mod.F90 (modified) (1 diff)
• testphys1d.F90 (modified) (3 diffs)
• writerestart1D_mod.F90 (modified) (3 diffs)

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

@@ -538 +538 @@
 else
     do iq = 1,nq
-        read(3,*,iostat = ierr) header, (qsurf(1,iq,j), j = 1,size(qsurf,3)), (q(1,ilayer,iq), ilayer = 1,nlayer)
+        read(3,*,iostat = ierr) header, (q(1,ilayer,iq), ilayer = 1,nlayer)
         if (ierr /= 0) error stop 'Not enough atmospheric layers defined in the file "'//trim(start1Dname)//'" for the tracer "'//trim(header)//'"!'
         if (trim(tname(iq)) /= trim(header)) then

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

@@ -129 +129 @@
 if (startfiles_1D) then
     inquire(file = 'start1D.txt',exist = therestart1D)
-    if (startfiles_1D .and. .not. therestart1D) then
+    if (.not. therestart1D) then
         write(*,*) 'There is no "start1D.txt" file!'
         write(*,*) 'Initialization is done with default values.'
@@ -139 +139 @@
     endif
 endif
+if (therestart1D .and. .not. therestartfi) error stop 'There is a "start1D.txt" but no "startfi.nc". There might be problems for the initialization (for example ''qsurf'')!'
 
 call init_testphys1d('start1D.txt','startfi.nc',therestart1D,therestartfi,ngrid,nlayer,odpref,nq,q, &
@@ -261 +262 @@
 
 ! Writing the "restart1D.txt" file for the next run
-if (startfiles_1D) call writerestart1D('restart1D.txt',psurf,pa,preff,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,q)
 
 write(*,*) "testphys1d: everything is cool!"

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

@@ -5 +5 @@
 contains
 
-SUBROUTINE writerestart1D(filename,psurf,pa,preff,tsurf,nlayer,nslope,temp,u,v,nq,qnames,qsurf,q)
+SUBROUTINE writerestart1D(filename,psurf,pa,preff,tsurf,nlayer,nslope,temp,u,v,nq,qnames,q)
 
 implicit none
@@ -16 +16 @@
 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
 
@@ -26 +25 @@
 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)
+    write(1,*) qnames(i), (q(il,i), il = 1,nlayer)
 enddo
 write(1,*) 'u', (u(il), il = 1,nlayer)