source: trunk/LMDZ.VENUS/libf/phyvenus/phyredem.F90 @ 1198

Last change on this file since 1198 was 973, checked in by slebonnois, 12 years ago

EM+SL: bug corrections in Venus physics

File size: 3.1 KB
Line 
1!
2! $Id: $
3!
4      SUBROUTINE phyredem (fichnom)
5
6      USE dimphy
7      USE mod_grid_phy_lmdz
8      USE mod_phys_lmdz_para
9      USE iophy
10      USE phys_state_var_mod
11      USE iostart, only : open_restartphy,close_restartphy, &
12                          put_var,put_field
13      USE infotrac
14      USE comgeomphy,  only: rlatd,rlond
15      USE control_mod, only: raz_date
16
17      implicit none
18!======================================================================
19! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
20! Objet: Lecture de l'etat initial pour la physique
21!======================================================================
22#include "dimensions.h"
23#include "netcdf.inc"
24#include "dimsoil.h"
25#include "clesphys.h"
26#include "tabcontrol.h"
27#include "temps.h"
28!======================================================================
29
30character(len=*),intent(in) :: fichnom
31REAL    :: tab_cntrl(length)
32integer :: isoil
33CHARACTER(len=2) :: str2
34
35
36! open file
37
38      CALL open_restartphy(fichnom)
39
40! tab_cntrl() contains run parameters
41
42      tab_cntrl(:)=0.0
43 
44      tab_cntrl(1) = dtime
45      tab_cntrl(2) = radpas
46      tab_cntrl(3) = 0.0
47      tab_cntrl(4) = solaire
48      tab_cntrl(5) = 0
49      tab_cntrl(6) = nbapp_rad
50
51      IF( cycle_diurne ) tab_cntrl( 7 ) = 1.
52      IF(   soil_model ) tab_cntrl( 8 ) = 1.
53      IF(     ok_orodr ) tab_cntrl(10 ) = 1.
54      IF(     ok_orolf ) tab_cntrl(11 ) = 1.
55      IF( ok_gw_nonoro ) tab_cntrl(12 ) = 1.
56
57      tab_cntrl(13) = day_end
58      tab_cntrl(14) = annee_ref
59      tab_cntrl(15) = itau_phy
60
61      CALL put_var("controle","Parametres de controle",tab_cntrl)
62
63! coordinates
64
65      CALL put_field("longitude", &
66                     "Longitudes de la grille physique",rlond)
67     
68      CALL put_field("latitude","Latitudes de la grille physique",rlatd)
69
70! variables
71
72      CALL put_field("TS","Temperature de surface",ftsol)
73
74      DO isoil=1, nsoilmx
75        IF (isoil.LE.99) THEN
76        WRITE(str2,'(i2.2)') isoil
77        CALL put_field("Tsoil"//str2, &
78                       "Temperature du sol No."//str2,ftsoil(:,isoil))
79        ELSE
80        PRINT*, "Trop de couches"
81        CALL abort
82        ENDIF
83      ENDDO
84
85      CALL put_field("ALBE","albedo de surface",falbe)
86      CALL put_field("solsw","Rayonnement solaire a la surface",solsw)
87      CALL put_field("sollw","Rayonnement IR a la surface",sollw)
88      CALL put_field("fder","Derive de flux",fder)
89      CALL put_field("dlw","Derivee flux IR",dlw)
90      CALL put_field("sollwdown","Flux IR vers le bas a la surface",sollwdown)
91      CALL put_field("RADS","Rayonnement net a la surface",radsol)
92      CALL put_field("ZMEA","zmea Orographie sous-maille",zmea)
93      CALL put_field("ZSTD","zstd Orographie sous-maille",zstd)
94      CALL put_field("ZSIG","zsig Orographie sous-maille",zsig)
95      CALL put_field("ZGAM","zgam Orographie sous-maille",zgam)
96      CALL put_field("ZTHE","zthe Orographie sous-maille",zthe)
97      CALL put_field("ZPIC","zpic Orographie sous-maille",zpic)
98      CALL put_field("ZVAL","zval Orographie sous-maille",zval)
99
100      CALL put_field("TANCIEN","T Previous iteration",t_ancien)
101
102! close file
103
104      CALL close_restartphy
105!$OMP BARRIER
106
107      END SUBROUTINE phyredem
Note: See TracBrowser for help on using the repository browser.