Ignore:
Timestamp:
Jul 22, 2024, 9:29:09 PM (4 months ago)
Author:
abarral
Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/phylmd/calcul_divers.h

    r4556 r5099  
    1 !
     1
    22! $Id$
    3 !
     3
    44! itap : nombre de pas de temps de la physique
    55! itapm1 : somme du nombre de pas de temps du/des mois precedent/s.
    66! (itap - itapm1) : nombre de pas de temps du mois courant
    7 !
     7
    88! Ne pas modifier les IFs ci-dessous impliquant itapm1,
    99! autrement les resultats seront faux !!
    10 !
     10
    1111! Ici on utilise MOD(itap - itapm1,NINT(mth_len*un_jour/phys_tstep)).EQ.1)
    1212! pour detecter le debut de chaque mois lorsque  l on tourne par an.
    13 !
     13
    1414! IM, 26.05.2023
    15 !
    16 !
     15
     16
    1717! Initialisations itapm1 du premier mois
    1818      IF(itap.EQ.1) THEN
     
    2121      ENDIF
    2222
    23 !
    2423! Initialisation debut de mois
    2524      IF(itap.EQ.itapm1+1) THEN
     
    2726!        print*,'Initialisation ndayrain_mth ',itap
    2827      ENDIF
    29 !
     28
    3029! Initialisation debut de chaque jour
    3130      IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.1) THEN
     
    3332!       print*,'initialisation mois suivants day_rain itap',itap
    3433      ENDIF
    35 !
     34
    3635! Calcul a chaque pas de temps de la physique
    3736         DO i = 1, klon
     
    3938            IF(total_rain(i).GT.0.) nday_rain(i)=1.
    4039         ENDDO
    41 !
     40
    4241! Cumul en fin de journee
    4342      IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN
     
    4645         ENDDO
    4746      ENDIF
    48 !
     47
    4948! Initialisation fin de mois
    5049! Ne pas changer le IF ci-dessous, car le compteur itapm1 est augmente
     
    5554!       print*,'fmois i i1 mjt',itap,itapm1,mth_len,un_jour,phys_tstep
    5655      ENDIF
    57 !
     56
    5857! calcul temperatures minimale et maximale moyennees sur le mois
    59 !
     58
    6059!initialisation debut de mois pour les fichiers mensuels annuels
    6160  IF(itap.EQ.itapm1+1) THEN
     
    6362     t2m_max_mon=0.
    6463  ENDIF
    65 !
     64
    6665!initialisation debut de journee pour les fichiers mensuels annuels
    6766  IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.1) THEN
     
    6968     zt2m_max_mon=zt2m
    7069  ENDIF
    71 !
     70
    7271!calcul sur tous les pas de temps pour les fichiers mensuels annuels
    7372     DO i = 1, klon
     
    7574        zt2m_max_mon(i)=MAX(zt2m(i),zt2m_max_mon(i))
    7675     ENDDO
    77 !
     76
    7877!fin journee
    7978  IF(MOD(itap,NINT(un_jour/phys_tstep)).EQ.0) THEN
     
    8180   t2m_max_mon=t2m_max_mon+zt2m_max_mon
    8281  ENDIF
    83 !
     82
    8483!fin mois
    8584  IF(itap==itapm1) THEN
     
    8786   t2m_max_mon=t2m_max_mon/mth_len
    8887  ENDIF
    89 !
Note: See TracChangeset for help on using the changeset viewer.