Changeset 3045 for trunk/LMDZ.MARS/libf/phymars/dyn1d/write_profile.F90

Timestamp:

Sep 21, 2023, 10:19:40 AM (14 months ago)

Author:

jbclement

Message:

Mars PCM 1D:
Addition of the pressure profile to be able to run chained simulations in 1D. Like this, the surface pressure can be got from one run to the next.
Rework of "write_profile.F90" + correction of arguments in its call in "testphys1d.F": 'qsurf' was wrongly given!
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/write_profile.F90 (modified) (1 diff)

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

@@ -1 @@
-SUBROUTINE writeprofile(nlev,data,name2,qsurf)
+SUBROUTINE writeprofile(nlev,profilename,surfdata,profiledata)
 
+implicit none
 
-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 nlev
-REAL data(nlev)
-REAL qsurf
-CHARACTER(len=30) :: name2
+! Local
+integer :: il
 
-!local
-INTEGER il
+! Write the data
+open(1,file = 'profile_out_'//trim(profilename),form = 'formatted')
+write(1,*) surfdata
+do il = 1,nlev
+    write(1,*) profiledata(il)
+enddo
+close(1)
+return
 
-!write the data
-OPEN(1,file='profile_out_'//trim(name2),form='formatted')
-write(1,*) qsurf
-DO il=1,nlev
-    write(1,*) data(il)
-ENDDO
-CLOSE(1)
-RETURN
 END