Changeset 3051 for trunk/LMDZ.MARS/libf/phymars/dyn1d/writerestart1D.F90

Timestamp:

Sep 26, 2023, 6:10:40 PM (16 months ago)

Author:

jbclement

Message:

Mars PCM 1D:
Addition of the file "start1D.txt" which mimics a "start.nc" file to be able to run chained simulations in 1D. Like this, key variables and profiles can be got from one run to the next. As a consequence, the subroutine "write_profile.F90" is replaced by "writerestart1D.F90".
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/writerestart1D.F90 (moved) (moved from trunk/LMDZ.MARS/libf/phymars/dyn1d/write_profile.F90) (1 diff)

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

@@ -1 @@
-SUBROUTINE writeprofile(nlev,profilename,surfdata,profiledata)
+SUBROUTINE writerestart1D(psurf,tsurf,nlayer,temp,u,v,nq,qnames,qsurf,q)
 
 implicit none
 
-! arguments
-integer, intent(in)               :: nlev
-real, intent(in)                  :: surfdata
-real, dimension(nlev), intent(in) :: profiledata
-character(len = 30), intent(in)   :: profilename
+! Arguments
+integer, intent(in)                           :: nlayer, nq
+real, intent(in)                              :: psurf, tsurf
+real, dimension(nlayer), intent(in)           :: temp, u, v
+real, dimension(nlayer,nq), intent(in)        :: q
+real, dimension(nq), intent(in)               :: qsurf
+character(len = *), dimension(nq), intent(in) :: qnames
 
-! Local
-integer :: il
+! Local variables
+integer :: il, iq
 
-! Write the data
-open(1,file = 'profile_out_'//trim(profilename),form = 'formatted')
-write(1,*) surfdata
-do il = 1,nlev
-    write(1,*) profiledata(il)
+! Write the data needed for a restart in "restart1D.txt"
+open(1,file = 'restart1D.txt',status = "replace",action = "write")
+do iq = 1,nq
+    write(1,*) qnames(iq), qsurf(iq), (q(il,iq), il = 1,nlayer)
 enddo
+write(1,*) 'ps', psurf
+write(1,*) 'u', (u(il), il = 1,nlayer)
+write(1,*) 'v', (v(il), il = 1,nlayer)
+write(1,*) 'teta', tsurf, (temp(il), il = 1,nlayer)
 close(1)
-return
 
-END
+END SUBROUTINE writerestart1D