source: LMDZ5/branches/LMDZ5_SPLA/libf/phylmd/calcul_divers.h @ 2404

Last change on this file since 2404 was 1912, checked in by musat, 11 years ago

1) Modifications pour faire des simulations par an avec un calendrier realiste (365 jours ou autre).

Il faut mettre une frequence de sortie de -1 (variable phys_out_filetimesteps dans config.def) pour
que IOIPSL calcule les moyennes mensuels en prenant en compte des longuers variables de chaque
mois. Par exemple, pour le fichier histmth (1er fichier) et histmthNMC (7eme).
phys_out_filetimesteps= -1 1day 6hr 6hr 6hr 1d -1 1day 6hr

2) Corrections titres variables niveaux de pression des fichiers histmth, histday, etc

  • 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 974 bytes
Line 
1!
2! $Header$
3!
4! Initialisations diverses au tout debut
5      IF(itap.EQ.1) THEN
6         itapm1=0
7!        surface terre
8         DO i=1, klon
9            IF(pctsrf(i,is_ter).GT.0.) THEN
10               paire_ter(i)=airephy(i)*pctsrf(i,is_ter)
11            ENDIF
12         ENDDO
13      ENDIF
14
15! Initialisation debut de mois
16      IF(itap.EQ.itapm1+1) THEN
17        nday_rain(:)=0.
18!       print*,'initialisation mois suivants day_rain itap',itap
19      ENDIF
20
21! Calcul fin de journee : total_rain, nday_rain
22      IF(MOD(itap,NINT(un_jour/dtime)).EQ.0) THEN
23!        print*,'calcul nday_rain itap ',itap
24         DO i = 1, klon
25            total_rain(i)=rain_fall(i)+snow_fall(i) 
26            IF(total_rain(i).GT.0.) nday_rain(i)=nday_rain(i)+1.
27         ENDDO
28      ENDIF
29
30! Initialisation fin de mois
31      IF(MOD(itap-itapm1,NINT(mth_len*un_jour/dtime)).EQ.0) THEN
32        itapm1=itapm1+NINT(mth_len*un_jour/dtime)
33!       print*,'initialisation itapm1 ',itapm1
34      ENDIF
Note: See TracBrowser for help on using the repository browser.