|
Last change
on this file since 2101 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.
