Changeset 2683 for trunk/LMDZ.VENUS/libf/phyvenus

Timestamp:

May 23, 2022, 6:25:07 PM (3 years ago)

Author:

flefevre

Message:

Venus: un peu de nettoyage dans l'initialisation de la chimie et de la microphysique

File:

• trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F (modified) (2 diffs)

trunk/LMDZ.VENUS/libf/phyvenus/physiq_mod.F

@@ -700 +700 @@
       ENDDO
 
-c---------
-c       Ecriture fichier initialisation
-c       PRINT*,'Ecriture Initial_State.csv'
-c       OPEN(88,file='Trac_Point.csv',
-c     & form='formatted')
-c---------
-     
-c---------
-c       Initialisation des parametres des nuages
-c===============================================
-     
-c MICROPHY SANS CHIMIE: seulement si full microphy (cl_scheme=2)
-
-      if (ok_chem.and..not.ok_cloud) then
-        print*,"LA CHIMIE A BESOIN DE LA MICROPHYSIQUE"
-        print*,"ok_cloud doit etre = a ok_chem"
-      stop
-      endif
-      if (.not.ok_chem.and.ok_cloud.and.(cl_scheme.eq.1)) then
-        print*,"cl_scheme=1 doesnot work without chemistry"
-      stop
-      endif
-       if (.not.ok_chem.and.ok_cloud.and.(cl_scheme.eq.2)) then
-        print*,"Full microphysics without chemistry"
-c indexation of microphys tracers
-        CALL chemparam_ini() 
-      endif
+!     initialisation of microphysical and chemical parameters
+
+      if (ok_chem .and. .not. ok_cloud) then
+         print*,"chemistry requires clouds"
+         print*,"ok_cloud must be .true."
+         stop
+      end if
+
+      if (.not. ok_chem .and. ok_cloud .and. (cl_scheme == 1)) then
+         print*,"cl_scheme = 1 requires chemistry"
+         print*,"ok_chem must be .true."
+         stop
+      end if
     
-! number of microphysical tracers
+!     number of microphysical tracers
 
       nmicro = 0
@@ -734 +720 @@
       if (ok_cloud .and. (cl_scheme == 2)) nmicro = 12
  
-c CAS 1D POUR MICROPHYS Aurelien
-      if ((nlon .EQ. 1) .AND. ok_cloud .AND. (cl_scheme.eq.1)) then
-        PRINT*,'Open profile_cloud_parameters.csv'
-        OPEN(66,file='profile_cloud_parameters.csv',
-     &   form='formatted')
-      endif
-
-      if ((nlon .EQ. 1) .AND. ok_sedim .AND. (cl_scheme.eq.1)) then
-        PRINT*,'Open profile_cloud_sedim.csv'
-        OPEN(77,file='profile_cloud_sedim.csv',
-     &   form='formatted')
-      endif
-           
-c INIT PHOTOCHEMISTRY ! includes the indexation of microphys tracers
-c     if ((nlon .GT. 1) .AND. ok_chem) then
-c !!! DONC 3D !!!  POURQUOI ???
-      if (ok_chem) then
-        CALL chemparam_ini() 
-      endif
+!     initialise chemical parameters. includes the indexation of microphys tracers
+
+      if (ok_chem .or. cl_scheme == 2) then
+         call chemparam_ini() 
+      end if
          
-c INIT MICROPHYS SCHEME 1 (AURELIEN)  
-      if ((nlon .GT. 1) .AND. ok_cloud .AND. (cl_scheme.eq.1)) then
-c !!! DONC 3D !!!
-        CALL cloud_ini(nlon,nlev)
-      endif
+!     initialise cloud parameters (for cl_scheme = 1)
+
+      if (ok_cloud .and. (cl_scheme == 1)) then
+         call cloud_ini(nlon,nlev)
+      end if
 
 !     initialise mmean
 
-      if(callthermos) then
+      if (callthermos) then
          call concentrations2(pplay,t,qx,nqmax)
-      endif
+      end if
 
 c========================