Changeset 3975

Timestamp:

Nov 21, 2025, 7:20:50 PM (3 weeks ago)

Author:

lrosset

Message:

TITAN PCM :
Updating and fixing the 1D version of the model : can now run with all couplings activated + XIOS outputs added.
Documentation about 1D version updated in LMDZ.TITAN/DOC, and reference files added in deftank.
LR

Location:

trunk/LMDZ.TITAN

Files:

• DOC/rcm1d.pdf (added)
• DOC/rcm1d.tex (added)
• DOC/rcm1d_old.tex (added)
• deftank/callphys.def (added)
• deftank/rcm1d.def (added)
• deftank/run.def.1d (added)
• libf/chimtitan/titan.h (modified) (1 diff)
• libf/dynphy_lonlat/phytitan/gr_kim_vervack.F90 (modified) (1 diff)
• libf/phytitan/dyn1d/mod_const_mpi.F90 (added)
• libf/phytitan/dyn1d/mod_interface_dyn_phys.F90 (added)
• libf/phytitan/dyn1d/parallel_lmdz.F90 (added)
• libf/phytitan/dyn1d/rcm1d.F (modified) (8 diffs)
• libf/phytitan/profile.F (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/chimtitan/titan.h

@@ -8 +8 @@
 #define R0    (double)(2575.0) /* Titan's radius */
 #define NLEV  (int)(133)  /* Nbre de niv verticaux -> Need to be coherent with the vertical grid used !! */
-#define NLD   (int)(40)   /* Nbre de niv verticaux faits sans diff -> Need to be coherent with the vertical grid used !! */
-//#define NLD   (int)(0)  /* -> This is for 1D, diffusion on all column */
+#ifdef CPP_1D
+    #define NLD   (int)(0)  /* -> This is for 1D, diffusion on all column */
+#else
+    #define NLD   (int)(40)   /* Nbre de niv verticaux faits sans diff -> Need to be coherent with the vertical grid used !! */
+#endif
 #define NLRT  (int)(650)  /* Nbre de niv verticaux dans table fmoy - aussi dans common_mod */
 

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/gr_kim_vervack.F90

@@ -63 +63 @@
   
   ! NB : At the top of the model we are in pure pressure coord. -> ap
-  ! ( except for 1D where we have only bp )
   
-  IF (jjm.GT.1) THEN
-    ptop = 2.0*aps(llm) - ap(llm)
-  ELSE
-    ptop = preff*(2.0*bps(llm) - bp(llm))
-  ENDIF
+  ptop = 2.0*aps(llm) - ap(llm)
   
   ! --------------------------------------------

trunk/LMDZ.TITAN/libf/phytitan/dyn1d/rcm1d.F

@@ -37 +37 @@
       use phys_state_var_mod, only: phys_state_var_init
       use physiq_mod, only: physiq
+! For XIOS outputs:
+      use mod_const_mpi, only: init_const_mpi
+      use parallel_lmdz, only: init_parallel
       implicit none
 
@@ -135 +138 @@
 c=======================================================================
 
+! For XIOS outputs
+#ifdef CPP_XIOS
+      call init_const_mpi
+      call init_parallel
+#endif
+
       ! read nq from traceur.def
       open(90,file='traceur.def',status='old',form='formatted',
@@ -204 +213 @@
         write(*,*) 'Please remove the file and restart the run.'
         write(*,*) 'Runtime parameters are supposed to be in rcm1d.def'
-        stop
+        !stop
       else
-        call system('touch run.def')
-        call system("echo 'INCLUDEDEF=callphys.def' >> run.def")
-        call system("echo 'INCLUDEDEF=rcm1d.def' >> run.def")
+        ! call system('touch run.def')
+        ! call system("echo 'INCLUDEDEF=callphys.def' >> run.def")
+        ! call system("echo 'INCLUDEDEF=rcm1d.def' >> run.def")
       endif
 
@@ -217 +226 @@
       write(*,*) " tracer = ",tracer
 
-! OK. now that run.def has been read once -- any variable is in memory.
-! so we can dump the dummy run.def
-!      call system("rm -rf run.def") ! Ehouarn: delay this to after inifis
-
 ! while we're at it, check if there is a 'traceur.def' file
-! and preocess it, if necessary. Otherwise initialize tracer names
+! and process it, if necessary. Otherwise initialize tracer names
       if (tracer) then
       ! load tracer names from file 'traceur.def'
@@ -449 +454 @@
       iphysiq=1 ! in 1D model physics are called evry time step
       ecritphy=day_step ! default value for ecritphy
-      PRINT *,'Nunber of steps between writediagfi ?'
+      PRINT *,'Number of steps between writediagfi ?'
       call getin("ecritphy",ecritphy)
       write(*,*) " ecritphy = ",ecritphy
@@ -503 +508 @@
 ! At this point, both getin() and getin_p() functions have been used,
 ! and the run.def file can be removed.
-      call system("rm -rf run.def")
+      ! call system("rm -rf run.def")
 
 !!! We check everything is OK.
@@ -941 +946 @@
 c    ========================================================
       if(saveprofile)then
-         OPEN(12,file='profile.out',form='formatted')
+         OPEN(12,file='profile_T.out',form='formatted')
          write(12,*) tsurf
          DO ilayer=1,llm
@@ -956 +961 @@
 c    ========================================================
       end                       !rcm1d
- 
-

trunk/LMDZ.TITAN/libf/phytitan/profile.F

@@ -178 +178 @@
       ! then profile of atmospheric temperature
       ELSE IF(ichoice.GE.8) THEN
-      OPEN(11,file='profile',status='old',form='formatted',err=101)
+      OPEN(11,file='profile_T',status='old',form='formatted',err=101)
       DO il=1,nlev
         READ (11,*) temp(il)
@@ -184 +184 @@
 
       GOTO 201
-101   STOP'fichier profile inexistant'
+101   STOP'fichier profile_T inexistant'
 201   CONTINUE
       CLOSE(10)