Changeset 1707 for LMDZ5/branches/testing/libf/phylmd/iniphysiq.F

Timestamp:

Jan 11, 2013, 10:19:19 AM (11 years ago)

Author:

Laurent Fairhead

Message:

Version testing basée sur la r1706

---

Testing release based on r1706

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/iniphysiq.F (modified) (6 diffs)

LMDZ5/branches/testing/libf/phylmd/iniphysiq.F

@@ -8 +8 @@
      $           pdayref,ptimestep,
      $           plat,plon,parea,pcu,pcv,
-     $           prad,pg,pr,pcpp)
-      USE dimphy
-      USE mod_grid_phy_lmdz
-      USE mod_phys_lmdz_para
-      USE comgeomphy
+     $           prad,pg,pr,pcpp,iflag_phys)
+      USE dimphy, only : klev
+      USE mod_grid_phy_lmdz, only : klon_glo
+      USE mod_phys_lmdz_para, only : klon_omp,klon_omp_begin,
+     &                               klon_omp_end,klon_mpi_begin
+      USE comgeomphy, only : airephy,cuphy,cvphy,rlond,rlatd
 
       IMPLICIT NONE
@@ -18 +19 @@
 c=======================================================================
 c
-c   subject:
-c   --------
+c   Initialisation of the physical constants and some positional and 
+c   geometrical arrays for the physics
 c
-c   Initialisation for the physical parametrisations of the LMD 
-c   martian atmospheric general circulation modele.
-c
-c   author: Frederic Hourdin 15 / 10 /93
-c   -------
-c
-c   arguments:
-c   ----------
-c
-c   input:
-c   ------
 c
 c    ngrid                 Size of the horizontal grid.
@@ -37 +27 @@
 c    nlayer                Number of vertical layers.
 c    pdayref               Day of reference for the simulation
-c    firstcall             True at the first call
-c    lastcall              True at the last call
-c    pday                  Number of days counted from the North. Spring
-c                          equinoxe.
 c
 c=======================================================================
-c
-c-----------------------------------------------------------------------
-c   declarations:
-c   -------------
  
 cym#include "dimensions.h"
@@ -52 +34 @@
 cym#include "comgeomphy.h"
 #include "YOMCST.h"
-      REAL prad,pg,pr,pcpp,punjours
- 
-      INTEGER ngrid,nlayer
-      REAL plat(ngrid),plon(ngrid),parea(klon_glo)
-      REAL pcu(klon_glo),pcv(klon_glo)
-      INTEGER pdayref
-      INTEGER :: ibegin,iend,offset
- 
-      REAL ptimestep
+#include "iniprint.h"
+
+      REAL,INTENT(IN) :: prad ! radius of the planet (m)
+      REAL,INTENT(IN) :: pg ! gravitational acceleration (m/s2)
+      REAL,INTENT(IN) :: pr ! ! reduced gas constant R/mu
+      REAL,INTENT(IN) :: pcpp ! specific heat Cp
+      REAL,INTENT(IN) :: punjours ! length (in s) of a standard day
+      INTEGER,INTENT(IN) :: ngrid ! number of horizontal grid points in the physics
+      INTEGER,INTENT(IN) :: nlayer ! number of atmospheric layers
+      REAL,INTENT(IN) :: plat(ngrid) ! latitudes of the physics grid
+      REAL,INTENT(IN) :: plon(ngrid) ! longitudes of the physics grid
+      REAL,INTENT(IN) :: parea(klon_glo) ! area (m2)
+      REAL,INTENT(IN) :: pcu(klon_glo) ! cu coeff. (u_covariant = cu * u)
+      REAL,INTENT(IN) :: pcv(klon_glo) ! cv coeff. (v_covariant = cv * v)
+      INTEGER,INTENT(IN) :: pdayref ! reference day of for the simulation
+      REAL,INTENT(IN) :: ptimestep !physics time step (s)
+      INTEGER,INTENT(IN) :: iflag_phys ! type of physics to be called
+
+      INTEGER :: ibegin,iend,offset 
       CHARACTER (LEN=20) :: modname='iniphysiq'
       CHARACTER (LEN=80) :: abort_message
  
       IF (nlayer.NE.klev) THEN
-         PRINT*,'STOP in inifis'
-         PRINT*,'Probleme de dimensions :'
-         PRINT*,'nlayer     = ',nlayer
-         PRINT*,'klev   = ',klev
+         write(lunout,*) 'STOP in ',trim(modname)
+         write(lunout,*) 'Problem with dimensions :'
+         write(lunout,*) 'nlayer     = ',nlayer
+         write(lunout,*) 'klev   = ',klev
          abort_message = ''
          CALL abort_gcm (modname,abort_message,1)
@@ -74 +66 @@
 
       IF (ngrid.NE.klon_glo) THEN
-         PRINT*,'STOP in inifis'
-         PRINT*,'Probleme de dimensions :'
-         PRINT*,'ngrid     = ',ngrid
-         PRINT*,'klon   = ',klon_glo
+         write(lunout,*) 'STOP in ',trim(modname)
+         write(lunout,*) 'Problem with dimensions :'
+         write(lunout,*) 'ngrid     = ',ngrid
+         write(lunout,*) 'klon   = ',klon_glo
          abort_message = ''
          CALL abort_gcm (modname,abort_message,1)
       ENDIF
-c$OMP PARALLEL PRIVATE(ibegin,iend) 
-c$OMP+         SHARED(parea,pcu,pcv,plon,plat)
+
+!$OMP PARALLEL PRIVATE(ibegin,iend) 
+!$OMP+         SHARED(parea,pcu,pcv,plon,plat)
       
       offset=klon_mpi_begin-1
@@ -92 +85 @@
       rlatd(1:klon_omp)=plat(offset+klon_omp_begin:offset+klon_omp_end)
 
+      ! suphel => initialize some physical constants (orbital parameters,
+      !           geoid, gravity, thermodynamical constants, etc.) in the
+      !           physics
       call suphel
+      
+!$OMP END PARALLEL
 
-c$OMP END PARALLEL
+      ! check that physical constants set in 'suphel' are coherent
+      ! with values set in the dynamics:
+      if (RDAY.ne.punjours) then
+        write(lunout,*) "iniphysiq: length of day discrepancy!!!"
+        write(lunout,*) "  in the dynamics punjours=",punjours
+        write(lunout,*) "   but in the physics RDAY=",RDAY
+        if (abs(RDAY-punjours).gt.0.01) then
+          ! stop here if the relative difference is more than 1%
+          abort_message = 'length of day discrepancy'
+          CALL abort_gcm (modname,abort_message,1)
+        endif
+      endif
+      if (RG.ne.pg) then
+        write(lunout,*) "iniphysiq: gravity discrepancy !!!"
+        write(lunout,*) "     in the dynamics pg=",pg
+        write(lunout,*) "  but in the physics RG=",RG
+        if (abs(RG-pg).gt.0.01) then
+          ! stop here if the relative difference is more than 1%
+          abort_message = 'gravity discrepancy'
+          CALL abort_gcm (modname,abort_message,1)
+        endif
+      endif
+      if (RA.ne.prad) then
+        write(lunout,*) "iniphysiq: planet radius discrepancy !!!"
+        write(lunout,*) "   in the dynamics prad=",prad
+        write(lunout,*) "  but in the physics RA=",RA
+        if (abs(RA-prad).gt.0.01) then
+          ! stop here if the relative difference is more than 1%
+          abort_message = 'planet radius discrepancy'
+          CALL abort_gcm (modname,abort_message,1)
+        endif
+      endif
+      if (RD.ne.pr) then
+        write(lunout,*)"iniphysiq: reduced gas constant discrepancy !!!"
+        write(lunout,*)"     in the dynamics pr=",pr
+        write(lunout,*)"  but in the physics RD=",RD
+        if (abs(RD-pr).gt.0.01) then
+          ! stop here if the relative difference is more than 1%
+          abort_message = 'reduced gas constant discrepancy'
+          CALL abort_gcm (modname,abort_message,1)
+        endif
+      endif
+      if (RCPD.ne.pcpp) then
+        write(lunout,*)"iniphysiq: specific heat discrepancy !!!"
+        write(lunout,*)"     in the dynamics pcpp=",pcpp
+        write(lunout,*)"  but in the physics RCPD=",RCPD
+        if (abs(RCPD-pcpp).gt.0.01) then
+          ! stop here if the relative difference is more than 1%
+          abort_message = 'specific heat discrepancy'
+          CALL abort_gcm (modname,abort_message,1)
+        endif
+      endif
 
-      print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ'
-      print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...'
+! Additional initializations for aquaplanets
+!$OMP PARALLEL
+      if (iflag_phys>=100) then
+        call iniaqua(klon_omp,rlatd,rlond,iflag_phys)
+      endif
+!$OMP END PARALLEL
 
-      RETURN
-9999  CONTINUE
-      abort_message ='Cette version demande les fichier rnatur.dat
-     & et surf.def'
-      CALL abort_gcm (modname,abort_message,1)
+!      RETURN
+!9999  CONTINUE
+!      abort_message ='Cette version demande les fichier rnatur.dat
+!     & et surf.def'
+!      CALL abort_gcm (modname,abort_message,1)
 
       END