Changeset 2160 for LMDZ5/branches/testing/libf/dyn3dmem/iniacademic_loc.F90

Timestamp:

Nov 28, 2014, 4:36:29 PM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes -r2070:2158 into testing branch. Compilation problems introduced by revision r2155 have been corrected by hand

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/dyn3dmem/iniacademic_loc.F90 (modified) (3 diffs)

LMDZ5/branches/testing/libf/dyn3dmem/iniacademic_loc.F90

@@ -4 +4 @@
 SUBROUTINE iniacademic_loc(vcov,ucov,teta,q,masse,ps,phis,time_0)
 
+  USE filtreg_mod, ONLY: inifilr
   use exner_hyb_m, only: exner_hyb
   use exner_milieu_m, only: exner_milieu
-  USE filtreg_mod
   USE infotrac, ONLY : nqtot
   USE control_mod, ONLY: day_step,planet_type
-  USE parallel_lmdz
+  USE parallel_lmdz, ONLY: ijb_u, ije_u, ijb_v, ije_v
 #ifdef CPP_IOIPSL
-  USE IOIPSL
+  USE IOIPSL, ONLY: getin
 #else
   ! if not using IOIPSL, we still need to use (a local version of) getin
-  USE ioipsl_getincom
+  USE ioipsl_getincom, ONLY: getin
 #endif
   USE Write_Field
@@ -41 +41 @@
   !   ----------
 
-  real time_0
-
-  !   variables dynamiques
-  REAL vcov(ijb_v:ije_v,llm),ucov(ijb_u:ije_u,llm) ! vents covariants
-  REAL teta(ijb_u:ije_u,llm)                 ! temperature potentielle
-  REAL q(ijb_u:ije_u,llm,nqtot)               ! champs advectes
-  REAL ps(ijb_u:ije_u)                       ! pression  au sol
-  REAL masse(ijb_u:ije_u,llm)                ! masse d'air
-  REAL phis(ijb_u:ije_u)                     ! geopotentiel au sol
+  REAL,INTENT(OUT) :: time_0
+
+  !   fields
+  REAL,INTENT(OUT) :: vcov(ijb_v:ije_v,llm) ! meridional covariant wind
+  REAL,INTENT(OUT) :: ucov(ijb_u:ije_u,llm) ! zonal covariant wind
+  REAL,INTENT(OUT) :: teta(ijb_u:ije_u,llm) ! potential temperature (K)
+  REAL,INTENT(OUT) :: q(ijb_u:ije_u,llm,nqtot) ! advected tracers (.../kg_of_air)
+  REAL,INTENT(OUT) :: ps(ijb_u:ije_u) ! surface pressure (Pa)
+  REAL,INTENT(OUT) :: masse(ijb_u:ije_u,llm) ! air mass in grid cell (kg)
+  REAL,INTENT(OUT) :: phis(ijb_u:ije_u) ! surface geopotential
 
   !   Local:
@@ -80 +81 @@
   character(len=80) :: abort_message
 
+  ! Sanity check: verify that options selected by user are not incompatible
+  if ((iflag_phys==1).and. .not. read_start) then
+    write(lunout,*) trim(modname)," error: if read_start is set to ", &
+    " false then iflag_phys should not be 1"
+    write(lunout,*) "You most likely want an aquaplanet initialisation", &
+    " (iflag_phys >= 100)"
+    call abort_gcm(modname,"incompatible iflag_phys==1 and read_start==.false.",1)
+  endif
+  
   !-----------------------------------------------------------------------
   ! 1. Initializations for Earth-like case