source: LMDZ5/branches/testing/libf/phydev/phyaqua.F @ 1910

Last change on this file since 1910 was 1910, checked in by Laurent Fairhead, 11 years ago

Merged trunk changes r1860:1909 into testing branch

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
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.