Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/calcul_divers.h
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/calcul_divers.h
r4556 r5099 1 ! 1 2 2 ! $Id$ 3 ! 3 4 4 ! itap : nombre de pas de temps de la physique 5 5 ! itapm1 : somme du nombre de pas de temps du/des mois precedent/s. 6 6 ! (itap - itapm1) : nombre de pas de temps du mois courant 7 ! 7 8 8 ! Ne pas modifier les IFs ci-dessous impliquant itapm1, 9 9 ! autrement les resultats seront faux !! 10 ! 10 11 11 ! Ici on utilise MOD(itap - itapm1,NINT(mth_len*un_jour/phys_tstep)).EQ.1) 12 12 ! pour detecter le debut de chaque mois lorsque l on tourne par an. 13 ! 13 14 14 ! IM, 26.05.2023 15 ! 16 ! 15 16 17 17 ! Initialisations itapm1 du premier mois 18 18 IF(itap.EQ.1) THEN … … 21 21 ENDIF 22 22 23 !24 23 ! Initialisation debut de mois 25 24 IF(itap.EQ.itapm1+1) THEN … … 27 26 ! print*,'Initialisation ndayrain_mth ',itap 28 27 ENDIF 29 ! 28 30 29 ! Initialisation debut de chaque jour 31 30 IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.1) THEN … … 33 32 ! print*,'initialisation mois suivants day_rain itap',itap 34 33 ENDIF 35 ! 34 36 35 ! Calcul a chaque pas de temps de la physique 37 36 DO i = 1, klon … … 39 38 IF(total_rain(i).GT.0.) nday_rain(i)=1. 40 39 ENDDO 41 ! 40 42 41 ! Cumul en fin de journee 43 42 IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN … … 46 45 ENDDO 47 46 ENDIF 48 ! 47 49 48 ! Initialisation fin de mois 50 49 ! Ne pas changer le IF ci-dessous, car le compteur itapm1 est augmente … … 55 54 ! print*,'fmois i i1 mjt',itap,itapm1,mth_len,un_jour,phys_tstep 56 55 ENDIF 57 ! 56 58 57 ! calcul temperatures minimale et maximale moyennees sur le mois 59 ! 58 60 59 !initialisation debut de mois pour les fichiers mensuels annuels 61 60 IF(itap.EQ.itapm1+1) THEN … … 63 62 t2m_max_mon=0. 64 63 ENDIF 65 ! 64 66 65 !initialisation debut de journee pour les fichiers mensuels annuels 67 66 IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.1) THEN … … 69 68 zt2m_max_mon=zt2m 70 69 ENDIF 71 ! 70 72 71 !calcul sur tous les pas de temps pour les fichiers mensuels annuels 73 72 DO i = 1, klon … … 75 74 zt2m_max_mon(i)=MAX(zt2m(i),zt2m_max_mon(i)) 76 75 ENDDO 77 ! 76 78 77 !fin journee 79 78 IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN … … 81 80 t2m_max_mon=t2m_max_mon+zt2m_max_mon 82 81 ENDIF 83 ! 82 84 83 !fin mois 85 84 IF(itap==itapm1) THEN … … 87 86 t2m_max_mon=t2m_max_mon/mth_len 88 87 ENDIF 89 !
Note: See TracChangeset
for help on using the changeset viewer.