Changeset 2311 for LMDZ5/trunk/libf/phylmd/inifis_mod.F90

Timestamp:

Jun 25, 2015, 9:45:24 AM (9 years ago)

Author:

Ehouarn Millour

Message:

Further modifications to enforce physics/dynamics separation:

• moved iniprint.h and misc_mod back to dyn3d_common, as these should only be used by dynamics.
• created print_control_mod in the physics to store flags prt_level, lunout, debug to be local to physics (should be used rather than iniprint.h)
• created abort_physic.F90 , which does the same job as abort_gcm() did, but should be used instead when in physics.
• reactivated inifis (turned it into a module, inifis_mod.F90) to initialize physical constants and print_control_mod flags.

EM

File:

• LMDZ5/trunk/libf/phylmd/inifis_mod.F90 (moved) (moved from LMDZ5/trunk/libf/phylmd/inifis.F90) (1 diff)

LMDZ5/trunk/libf/phylmd/inifis_mod.F90

@@ +1 @@
+! $Id$
+MODULE inifis_mod
 
-! $Id$
+CONTAINS
 
-SUBROUTINE inifis(ngrid, nlayer, punjours, pdayref, ptimestep, plat, plon, &
-    parea, prad, pg, pr, pcpp)
-  USE dimphy
+  SUBROUTINE inifis(punjours, prad, pg, pr, pcpp)
+  ! Initialize some physical constants and settings
+  USE print_control_mod, ONLY: init_print_control, lunout
   IMPLICIT NONE
 
-  ! =======================================================================
+  include "YOMCST.h"
+  REAL,INTENT(IN) :: prad, pg, pr, pcpp, punjours
 
-  ! subject:
-  ! --------
-
-  ! Initialisation for the physical parametrisations of the LMD
-  ! martian atmospheric general circulation modele.
-
-  ! author: Frederic Hourdin 15 / 10 /93
-  ! -------
-
-  ! arguments:
-  ! ----------
-
-  ! input:
-  ! ------
-
-  ! ngrid                 Size of the horizontal grid.
-  ! All internal loops are performed on that grid.
-  ! nlayer                Number of vertical layers.
-  ! pdayref               Day of reference for the simulation
-  ! firstcall             True at the first call
-  ! lastcall              True at the last call
-  ! pday                  Number of days counted from the North. Spring
-  ! equinoxe.
-
-  ! =======================================================================
-
-  ! -----------------------------------------------------------------------
-  ! declarations:
-  ! -------------
-
-  ! ym#include "dimensions.h"
-  ! ym#include "dimphy.h"
-
-  include 'iniprint.h'
-  REAL prad, pg, pr, pcpp, punjours
-
-  INTEGER ngrid, nlayer
-  REAL plat(ngrid), plon(ngrid), parea(klon)
-  INTEGER pdayref
-
-  REAL ptimestep
   CHARACTER (LEN=20) :: modname = 'inifis'
   CHARACTER (LEN=80) :: abort_message
 
+  ! Initialize flags lunout, prt_level, debug
+  CALL init_print_control
 
-  IF (nlayer/=klev) THEN
-    PRINT *, 'STOP in inifis'
-    PRINT *, 'Probleme de dimensions :'
-    PRINT *, 'nlayer     = ', nlayer
-    PRINT *, 'klev   = ', klev
-    abort_message = ''
-    CALL abort_gcm(modname, abort_message, 1)
+  ! suphel => initialize some physical constants (orbital parameters,
+  !           geoid, gravity, thermodynamical constants, etc.) in the
+  !           physics
+  CALL suphel
+
+  ! check that physical constants set in 'suphel' are coherent
+  ! with values set in the dynamics:
+  IF (rday/=punjours) THEN
+    WRITE (lunout, *) 'inifis: length of day discrepancy!!!'
+    WRITE (lunout, *) '  in the dynamics punjours=', punjours
+    WRITE (lunout, *) '   but in the physics RDAY=', rday
+    IF (abs(rday-punjours)>0.01*punjours) THEN
+        ! stop here if the relative difference is more than 1%
+      abort_message = 'length of day discrepancy'
+      CALL abort_physic(modname, abort_message, 1)
+    END IF
+  END IF
+  IF (rg/=pg) THEN
+    WRITE (lunout, *) 'inifis: gravity discrepancy !!!'
+    WRITE (lunout, *) '     in the dynamics pg=', pg
+    WRITE (lunout, *) '  but in the physics RG=', rg
+    IF (abs(rg-pg)>0.01*pg) THEN
+        ! stop here if the relative difference is more than 1%
+      abort_message = 'gravity discrepancy'
+      CALL abort_physic(modname, abort_message, 1)
+    END IF
+  END IF
+  IF (ra/=prad) THEN
+    WRITE (lunout, *) 'inifis: planet radius discrepancy !!!'
+    WRITE (lunout, *) '   in the dynamics prad=', prad
+    WRITE (lunout, *) '  but in the physics RA=', ra
+    IF (abs(ra-prad)>0.01*prad) THEN
+        ! stop here if the relative difference is more than 1%
+      abort_message = 'planet radius discrepancy'
+      CALL abort_physic(modname, abort_message, 1)
+    END IF
+  END IF
+  IF (rd/=pr) THEN
+    WRITE (lunout, *) 'inifis: reduced gas constant discrepancy !!!'
+    WRITE (lunout, *) '     in the dynamics pr=', pr
+    WRITE (lunout, *) '  but in the physics RD=', rd
+    IF (abs(rd-pr)>0.01*pr) THEN
+        ! stop here if the relative difference is more than 1%
+      abort_message = 'reduced gas constant discrepancy'
+      CALL abort_physic(modname, abort_message, 1)
+    END IF
+  END IF
+  IF (rcpd/=pcpp) THEN
+    WRITE (lunout, *) 'inifis: specific heat discrepancy !!!'
+    WRITE (lunout, *) '     in the dynamics pcpp=', pcpp
+    WRITE (lunout, *) '  but in the physics RCPD=', rcpd
+    IF (abs(rcpd-pcpp)>0.01*pcpp) THEN
+        ! stop here if the relative difference is more than 1%
+      abort_message = 'specific heat discrepancy'
+      CALL abort_physic(modname, abort_message, 1)
+    END IF
   END IF
 
-  IF (ngrid/=klon) THEN
-    PRINT *, 'STOP in inifis'
-    PRINT *, 'Probleme de dimensions :'
-    PRINT *, 'ngrid     = ', ngrid
-    PRINT *, 'klon   = ', klon
-    abort_message = ''
-    CALL abort_gcm(modname, abort_message, 1)
-  END IF
-
-  RETURN
-  abort_message = 'Cette version demande les fichier rnatur.dat &
-    &                                                         &
-    &        et surf.def'
-  CALL abort_gcm(modname, abort_message, 1)
-
-END SUBROUTINE inifis
+  END SUBROUTINE inifis
+  
+END MODULE inifis_mod