source: LMDZ5/trunk/libf/phydev/phyaqua.F @ 1907

Last change on this file since 1907 was 1907, checked in by lguez, 10 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.

  • 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.