Changeset 3993

Timestamp:

Jan 5, 2026, 8:31:25 AM (7 days ago)

Author:

emillour

Message:

Venus PCM:
Improving rcm1d: make reading/loading of tracer mixing ratios from
file init_1D.txt more flexible: identify tracers by their names so
that the ordering in init_1D.txt does not matter any more and that
missing tracers are simply set to 0.
EM

File:

• trunk/LMDZ.VENUS/libf/phyvenus/dyn1d/rcm1d.F (modified) (10 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/dyn1d/rcm1d.F

@@ -1 +1 @@
       PROGRAM rcm1d
       
-      USE infotrac
+      USE infotrac, only: infotrac_init, nqtot, tname
       use control_mod, only: planet_type,day_step,cpofT
       USE phys_state_var_mod
@@ -164 +164 @@
 c    Date (en sols depuis le solstice de printemps) du debut du run
       day0 = 0
-      PRINT *,'date de depart ?'
+      PRINT *,'start date ?'
       READ(unit,*) day0
       day=REAL(day0)
       PRINT *,day0
 c  Heure de demarrage
-      PRINT *,'heure de debut de simulation (entre 0 et 24) ?'
+      PRINT *,'starting hour of the simulation (between 0 et 24) ?'
       READ(unit,*) time
       time=time/24.E+0
@@ -179 +179 @@
       nlevel=nlayer+1
       nsoil=nsoilmx
-      PRINT *,'nombre de pas de temps par jour ?'
+      PRINT *,'nunmber of time steps per day ?'
       READ(unit,*) day_step
       print*,day_step
@@ -189 +189 @@
       nbapp_rad = 24000
 
-      PRINT *,'nombre de jours simules ?'
+      PRINT *,'number of days to run ?'
       READ(unit,*) nb_days
       print*,nb_days
@@ -201 +201 @@
 c ------------------------------------
 c
-      PRINT *,'pression au sol'
+      PRINT *,'surface pressure'
       READ(unit,*) psurf
       PRINT *,psurf
@@ -212 +212 @@
 c  latitude/longitude 
 c  -------------------
-      PRINT *,'latitude en degres ?'
+      PRINT *,'latitude in degrees ?'
       READ(unit,*) lati(1)
       PRINT *,lati(1)
@@ -253 +253 @@
  
 c    vent geostrophique
-      PRINT *,'composante vers l est du vent geostrophique (U) ?'
+      PRINT *,'eastward component of the geostrophic wind (U) ?'
       READ(unit,*) gru
-      PRINT *,'composante vers le nord du vent geostrophique (V) ?'
+      PRINT *,'northward component of the geostrophic wind (V) ?'
       READ(unit,*) grv
 
@@ -313 +313 @@
       CLOSE(unit)
 
-      print*,"               Pression        Altitude     Temperature"
+      print*,"               Pressure        Altitude     Temperature"
       ilayer=1
       ftsol(1)=tmp2(0)
@@ -349 +349 @@
       inquire(file="init_1D.txt",exist=file_is_present)
       if (file_is_present) then
-        open(21, form = 'formatted', file = 'init_1D.txt')
-        read(21,*)
-        do ilayer = nlayer,1,-1
-          read(21,*) idummy, dummy, dummy, (q(ilayer,iq), iq = 1,nqtot)
-!        print*, idummy, q(ilayer,1), q(ilayer,nqtot)
-        end do
-        close(21)
+        call read_init_1D(nlayer,nqtot,tname,q)
       else
         write(*,*) "Cannot find input file init_1D.txt!"
@@ -541 +535 @@
       ENDDO
 
+      CONTAINS
+
+        SUBROUTINE read_init_1D(nlayer,nqtot,tname,q)
+          IMPLICIT NONE
+          
+          integer,intent(in) :: nlayer
+          integer,intent(in) :: nqtot
+          character(len=*),intent(in) :: tname(nqtot)
+          real,intent(out) :: q(nlayer,nqtot)
+        
+          character(len=512) :: line ! to store a line of text
+          character(len=22) :: words(nqtot+2) ! to store items in line of text
+          integer :: words_b(nqtot+2),words_e(nqtot+2)
+          logical :: new_word
+          integer :: iq(nqtot) ! tracer index
+          real :: lineq(nqtot+2) ! values in a line of file
+          integer :: idummy
+          integer :: i,j,ind
+          
+          !1. Open file and process first line
+          open(21,form='formatted',file='init_1D.txt')
+          read(21,'(a)') line
+
+          ! identify words in line
+          ind=0
+          do i=1,len_trim(line)
+            if (line(i:i) == ' ') then
+              new_word=.true.
+            else
+              if (new_word) then
+                ind=ind+1
+                words_b(ind)=i
+                new_word=.false.
+              endif
+              words_e(ind)=i
+            endif
+          enddo ! of do i=1,len_trim(line)
+          ! write(*,*) "found ",ind," words"
+          do j=1,ind
+            write(words(j),"(a)") line(words_b(j):words_e(j))
+          enddo
+          
+          ! identify tracer index (if present in file)
+          iq(:)=0
+          do i=1,nqtot
+            do j=1,ind
+              ! look for match
+              if (trim(tname(i))==words(j)) then
+                iq(i)=j
+                exit ! move on to next tracer
+              endif
+            enddo
+            ! write(*,*) "i=",i,"iq(i)=",iq(i)
+          enddo
+          
+          ! 2. Read following lines and assign tracer mixing ratios
+          do ilayer=nlayer,1,-1
+            read(21,*) idummy,(lineq(j),j=1,ind)
+            ! reorder/store values
+            do i=1,nqtot
+              if (iq(i)/=0) then
+                q(ilayer,i)=lineq(iq(i))
+              else
+                ! set tracer value to zero since it is not in the file
+                q(ilayer,i)=0
+              endif
+            enddo
+          enddo ! of do lay=nlayer,1,-1
+          
+          ! 3. close file
+          close(21)
+          
+        END SUBROUTINE read_init_1D
+
 c    ========================================================
-      END
- 
-c***********************************************************************
-c***********************************************************************
-
-!#include "../dyn3d_common/disvert_noterre.F"
-!#include "../dyn3d/abort_gcm.F"
-
+      END PROGRAM rcm1d