Changeset 3254 for trunk

Timestamp:

Mar 6, 2024, 5:00:24 PM (9 months ago)

Author:

emillour

Message:

Venus PCM:
Some minor fixes to rcm1d (uninitialized variables) and add possibility to run a
fraction of a Venus day and not just an integer number of days.
EM

File:

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

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

@@ -51 +51 @@
       INTEGER ilayer,ilevel,isoil,idt,iq,i
       LOGICAl firstcall,lastcall
+      REAL :: nb_days ! number of Vdays (and/or fraction thererof) to run
 c
       INTEGER day0          ! date initial (sol ; =0 a Ls=0)
@@ -85 +86 @@
 
 !     initialisation des traceurs
-
+      logical :: file_is_present
       integer :: idummy
       real :: dummy
@@ -195 +196 @@
 
       PRINT *,'nombre de jours simules ?'
-      READ(unit,*) ndt
-      print*,ndt
-
-      ndt=ndt*day_step     
+      READ(unit,*) nb_days
+      print*,nb_days
+
+      ndt=nint(nb_days*day_step)  
+      write(*,*) " => will run ", ndt," timesteps"   
       dtphys=daysec/day_step  
       dtime=dtphys
@@ -265 +267 @@
          u(ilayer)=gru
          v(ilayer)=grv
+         w(ilayer)=0
       ENDDO
 
@@ -346 +349 @@
       ENDDO
 
-      print*,"lecture des profils chimiques"
-      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*,"rcm1d: Loading chemistry profiles from init_1D.txt"
+      ! check if the file is indeed there
+      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)
+        end do
+        close(21)
+      else
+        write(*,*) "Cannot find input file init_1D.txt!"
+        write(*,*) "Might as well stop here"
+        stop
+      endif ! of if(file_is_present)
 
 c    Initialisation des parametres d'oro
@@ -380 +391 @@
       sollw(1)    = 0.
       fder(1)     = 0.
+      dlw(1)      = 0.
+      sollwdown(1)= 0.
       radsol(1)   = 0.
-     
+
+      t_ancien(1,:)=0.
+      q2(1,:)=0.
+
       radpas      = NINT(1.*day_step/nbapp_rad)
       soil_model  = .true.
@@ -402 +418 @@
 c     Ouverture du fichier d'écriture des VMR
 
-      OPEN(5,file='chem.txt',form = 'formatted')
-      WRITE(5,'(2x,100a12)')'hpa',(specname(i), i=1,nqtot)
+      OPEN(15,file='chem.txt',form = 'formatted')
+      WRITE(15,'(2x,100a12)')'hpa',(specname(i), i=1,nqtot)
 
 !     debut de boucle temporelle
 
-      DO idt=1,48000
+      DO idt=1,ndt
         IF (idt.eq.ndt) then
          lastcall=.true.
@@ -516 +532 @@
          if (mod(idt,9600)==0) then   
               DO ilayer=1,nlayer 
-                 write (5,'(100e12.4)')play(ilayer)/100., q(ilayer,:)
+                 write (15,'(100e12.4)')play(ilayer)/100., q(ilayer,:)
      $                 *mmean(1,ilayer)/mmol(:)
               ENDDO   
@@ -522 +538 @@
       ENDDO   ! fin de la boucle temporelle
 
-      close(5)
+      close(15)
 c    ========================================================
 c    GESTION DES SORTIE