Changeset 1671 for LMDZ5/trunk/libf/phydev

Timestamp:

Oct 24, 2012, 9:10:10 AM (12 years ago)

Author:

Ehouarn Millour

Message:

• fixed "aquaplanet case" so that initializations (creation of files startphy.nc and limit.nc) now also works in parallel (mpi,omp,mixed).
• call to "iniaqua" is now done from within "iniphysiq" ; also added some tests to check consistency between essential variables shared by dynamics and physics (planetary radius, gravity, Cp, ...)
• simillarily adapted "phydev" routines, and added necessary routines to also be able to generate startphy.nc/restartphy.nc files there. Also removed common file "comcstphy.h" and replaced it with a module "comcstphy.F90"

EM

Location:

LMDZ5/trunk/libf/phydev

Files:

• comcstphy.F90 (added)
• comcstphy.h (deleted)
• iniphysiq.F (modified) (6 diffs)
• iostart.F90 (added)
• phyaqua.F (modified) (1 diff)
• phyredem.F (added)
• phys_state_var_mod.F90 (added)
• physiq.F90 (modified) (3 diffs)

LMDZ5/trunk/libf/phydev/iniphysiq.F

@@ -2 +2 @@
 ! $Id: iniphysiq.F 1403 2010-07-01 09:02:53Z fairhead $
 !
-c
-c
       SUBROUTINE iniphysiq(ngrid,nlayer,
      $           punjours,
      $           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
+      USE comcstphy, only : rradius,rg,rr,rcpp
 
       IMPLICIT NONE
@@ -18 +18 @@
 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 +26 @@
 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"
 cym#include "dimphy.h"
 cym#include "comgeomphy.h"
-#include "comcstphy.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
+#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
- 
-      REAL ptimestep
       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 +65 @@
 
       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 +84 @@
       rlatd(1:klon_omp)=plat(offset+klon_omp_begin:offset+klon_omp_end)
 
-!     call suphel
-!     prad,pg,pr,pcpp
+! copy some fundamental parameters to physics
       rradius=prad
       rg=pg
@@ -99 +90 @@
       rcpp=pcpp
 
-      
+!$OMP END PARALLEL
 
-c$OMP END PARALLEL
+!      print*,'ATTENTION !!! TRAVAILLER SUR INIPHYSIQ'
+!      print*,'CONTROLE DES LATITUDES, LONGITUDES, PARAMETRES ...'
 
-      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

LMDZ5/trunk/libf/phydev/phyaqua.F

@@ -1 @@
-! Routines complementaires pour la physique planetaire.
-
+!
+! $Id: $
+!
 
       subroutine iniaqua(nlon,latfi,lonfi,iflag_phys)
 
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-!  Creation d'un etat initial et de conditions aux limites 
-!  (resp startphy.nc et limit.nc) pour des configurations idealisees 
-! du modele LMDZ dans sa version terrestre.
-!  iflag_phys est un parametre qui controle
-!  iflag_phys = N  
-!    de 100 a 199 : aqua planetes avec SST forcees
-!                 N-100 determine le type de SSTs
-!    de 200 a 299 : terra planetes avec Ts calcule
-!        
+!  Create an initial state (startphy.nc) for the physics
+!  Usefull for idealised cases (e.g. aquaplanets or testcases)
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
+      use phys_state_var_mod, only : rlat,rlon,
+     &                               phys_state_var_init
+      use mod_phys_lmdz_para, only : klon_omp
+      use comgeomphy, only : rlond,rlatd
+      implicit none
+      
+      integer,intent(in) :: nlon,iflag_phys
+      real,intent(in) :: lonfi(nlon),latfi(nlon)
 
-      integer nlon,iflag_phys
-cIM ajout latfi, lonfi
-      REAL, DIMENSION (nlon) :: lonfi, latfi
+! local variables
+      real :: pi
+
+! initializations:
+      pi=2.*asin(1.)
+
+      call phys_state_var_init()
+
+      rlat(1:klon_omp)=rlatd(1:klon_omp)*180./pi
+      rlon(1:klon_omp)=rlond(1:klon_omp)*180./pi
 
 
-      return
+! Here you could create an initial condition for the physics
+! ...
+! ... fill in the fields...
+! ...
+! ... and create a "startphy.nc" file
+!      CALL phyredem ("startphy.nc")
+
       end
 

LMDZ5/trunk/libf/phydev/physiq.F90

@@ -11 +11 @@
      &            , PVteta)
 
-      USE dimphy
-      USE infotrac
-      USE comgeomphy
+      USE dimphy, only : klon,klev
+      USE infotrac, only : nqtot
+      USE comgeomphy, only : rlatd
+      USE comcstphy, only : rg
 
       IMPLICIT none
@@ -50 +51 @@
 !======================================================================
 #include "dimensions.h"
-#include "comcstphy.h"
+!#include "comcstphy.h"
 
       integer jjmp1
@@ -101 +102 @@
     PARAMETER    ( longcles = 20 )
 
-real temp_newton(klon,klev)
+real :: temp_newton(klon,klev)
 integer k
 logical, save :: first=.true.
+!$OMP THREADPRIVATE(first)
 
       REAL clesphy0( longcles      )
 
-d_u=0.
-d_v=0.
-d_t=0.
-d_qx=0.
-d_ps=0.
+! initializations
+if (first) then
+! ...
 
-     d_u(:,1)=-u(:,1)/86400.
-     do k=1,klev
-        temp_newton(:,k)=280.+cos(rlatd(:))*40.-pphi(:,k)/rg*6.e-3
-        d_t(:,k)=(temp_newton(:,k)-t(:,k))/1.e5
-     enddo
+  first=.false.
+endif
+
+! set all tendencies to zero
+d_u(:,:)=0.
+d_v(:,:)=0.
+d_t(:,:)=0.
+d_qx(:,:,:)=0.
+d_ps(:)=0.
+
+! compute tendencies to return to the dynamics:
+! "friction" on the first layer
+d_u(:,1)=-u(:,1)/86400.
+! newtonian rlaxation towards temp_newton()
+do k=1,klev
+  temp_newton(:,k)=280.+cos(rlatd(:))*40.-pphi(:,k)/rg*6.e-3
+  d_t(:,k)=(temp_newton(:,k)-t(:,k))/1.e5
+enddo
 
 
-      print*,'COUCOU PHYDEV'
-      return
-      end
+print*,'COUCOU PHYDEV'
+
+! if lastcall, then it is time to write "restartphy.nc" file
+if (lafin) then
+  call phyredem("restartphy.nc")
+endif
+
+end