Changeset 1707 for LMDZ5/branches/testing/libf/phydev/phyaqua.F

Timestamp:

Jan 11, 2013, 10:19:19 AM (12 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/phydev/phyaqua.F (modified) (1 diff)

LMDZ5/branches/testing/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