Changeset 2366 for trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/start2archive.F

Timestamp:

Jun 11, 2020, 7:40:22 PM (5 years ago)

Author:

jvatant

Message:

Titan GCM : Major maintenance catching up commits from the generic including :

• r2356 and 2354 removing obsolete old dynamical core
• various minor addition to physics and gestion of phys_state_var_mode, especially in dyn1d
• adding MESOSCALE CPP keys around chemistry and microphysics (disabled in mesoscale for now)

File:

• trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/start2archive.F (modified) (6 diffs)

trunk/LMDZ.TITAN/libf/dynphy_lonlat/phytitan/start2archive.F

@@ -22 +22 @@
       USE comsoil_h
       USE comchem_h, only : cnames, nkim, nlaykim_up, preskim, ykim_up
-!      USE comgeomfi_h, ONLY: lati, long, area
-!      use control_mod
-!      use comgeomphy, only: initcomgeomphy
-! to use  'getin'
-      USE ioipsl_getincom
+      USE ioipsl_getincom, only: getin
       USE planete_mod, only: year_day
       USE mod_const_mpi, ONLY: COMM_LMDZ
@@ -35 +31 @@
       USE temps_mod, ONLY: day_ini
       USE iniphysiq_mod, ONLY: iniphysiq
+      use phys_state_var_mod, only: phys_state_var_init
       use phyetat0_mod, only: phyetat0
       use tracer_h
+      use exner_hyb_m, only: exner_hyb
       implicit none
 
@@ -44 +42 @@
       include "comdissip.h"
       include "comgeom.h"
-!#include "control.h"
-
-!#include "dimphys.h"
-!#include "planete.h"
-!#include"advtrac.h"
+
       include "netcdf.inc"
 c-----------------------------------------------------------------------
@@ -62 +56 @@
       REAL pk(ip1jmp1,llm)
       REAL pkf(ip1jmp1,llm)
-      REAL beta(iip1,jjp1,llm)
       REAL phis(ip1jmp1)                     ! geopotentiel au sol
       REAL masse(ip1jmp1,llm)                ! masse de l'atmosphere
@@ -262 +255 @@
       Lmodif=0
 
+! Allocate saved arrays (as in firstcall of physiq)
+      call phys_state_var_init(nqtot)
+      
 ! Initialize tracer names, indexes and properties
       CALL initracer2(nqtot,tname)
@@ -312 +308 @@
 
       CALL pression(ip1jmp1, ap, bp, ps, p3d)
-      call exner_hyb(ip1jmp1, ps, p3d, beta, pks, pk, pkf)
+      call exner_hyb(ip1jmp1, ps, p3d, pks, pk, pkf)
 
 c=======================================================================