source: LMDZ5/branches/LMDZ5-DOFOCO/libf/phydev/phyaqua.F @ 3574

Last change on this file since 3574 was 1686, checked in by Ehouarn Millour, 12 years ago

Some additional stuff for "phydev": add what is required to be potentialy able to load a "startphy.nc" and also add illustrative examples of writting outputs in the physics using IOIPSL.
EM

File size: 1.0 KB
Line 
1!
2! $Id: $
3!
4
5      subroutine iniaqua(nlon,latfi,lonfi,iflag_phys)
6
7!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
8!  Create an initial state (startphy.nc) for the physics
9!  Usefull for idealised cases (e.g. aquaplanets or testcases)
10!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
11
12      use phys_state_var_mod, only : rlat,rlon,
13     &                               phys_state_var_init
14      use mod_phys_lmdz_para, only : klon_omp
15      use comgeomphy, only : rlond,rlatd
16      implicit none
17     
18      integer,intent(in) :: nlon,iflag_phys
19      real,intent(in) :: lonfi(nlon),latfi(nlon)
20
21! local variables
22      real :: pi
23
24! initializations:
25      pi=2.*asin(1.)
26
27      call phys_state_var_init()
28
29      rlat(1:klon_omp)=rlatd(1:klon_omp)*180./pi
30      rlon(1:klon_omp)=rlond(1:klon_omp)*180./pi
31
32
33! Here you could create an initial condition for the physics
34! ...
35! ... fill in the fields...
36! ...
37! ... and create a "startphy.nc" file
38      CALL phyredem ("startphy.nc")
39
40      end
41
Note: See TracBrowser for help on using the repository browser.