Changeset 655


Ignore:
Timestamp:
Jun 22, 2005, 2:51:34 PM (19 years ago)
Author:
lmdzadmin
Message:

Correction dans l'appel de chemhook_end (INCA) AC
MAF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/physiq.F

    r654 r655  
    28942894#include "calcul_divers.h"
    28952895c
    2896 c=============================================================
    2897 c
    2898 c Convertir les incrementations en tendances
    2899 c
    2900       DO k = 1, klev
    2901       DO i = 1, klon
    2902          d_u(i,k) = ( u_seri(i,k) - u(i,k) ) / dtime
    2903          d_v(i,k) = ( v_seri(i,k) - v(i,k) ) / dtime
    2904          d_t(i,k) = ( t_seri(i,k)-t(i,k) ) / dtime
    2905          d_qx(i,k,ivap) = ( q_seri(i,k) - qx(i,k,ivap) ) / dtime
    2906          d_qx(i,k,iliq) = ( ql_seri(i,k) - qx(i,k,iliq) ) / dtime
    2907       ENDDO
    2908       ENDDO
    2909 c
    2910       IF (nqmax.GE.3) THEN
    2911       DO iq = 3, nqmax
    2912       DO  k = 1, klev
    2913       DO  i = 1, klon
    2914          d_qx(i,k,iq) = ( tr_seri(i,k,iq-2) - qx(i,k,iq) ) / dtime
    2915       ENDDO
    2916       ENDDO
    2917       ENDDO
    2918       ENDIF
    2919 c
    2920 cIM rajout diagnostiques bilan KP pour analyse MJO par Jun-Ichi Yano
    2921 c#include "write_bilKP_ins.h"
    2922 c#include "write_bilKP_ave.h"
    2923 c
    2924 c Sauvegarder les valeurs de t et q a la fin de la physique:
    2925 c
    2926       DO k = 1, klev
    2927       DO i = 1, klon
    2928          t_ancien(i,k) = t_seri(i,k)
    2929          q_ancien(i,k) = q_seri(i,k)
    2930       ENDDO
    2931       ENDDO
    2932 c
    2933 c 22.03.04 BEG
    2934 c=============================================================
    2935 c   Ecriture des sorties
    2936 c=============================================================
    2937 #ifdef CPP_IOIPSL
    2938 
    2939 #ifdef histhf
    2940 #include "write_histhf.h"
    2941 #endif
    2942 
    2943 #ifdef histday
    2944 #include "write_histday.h"
    2945 #include "write_histday_seri.h"
    2946 #endif
    2947 
    2948 #ifdef histmth
    2949 #include "write_histmth.h"
    2950 #endif
    2951 
    2952 #ifdef histins
    2953 #include "write_histins.h"
    2954 #endif
    2955 
    2956 #ifdef histREGDYN
    2957 #include "write_histREGDYN.h"
    2958 #endif
    2959 
    2960 #ifdef histISCCP
    2961 #include "write_histISCCP.h"
    2962 #endif
    2963 
    2964 
    2965 #ifdef histmthNMC
    2966 #include "write_histmthNMC.h"
    2967 #endif
    2968 
    2969 #endif
    2970 
    29712896#ifdef INCA
    29722897#ifdef INCAINFO
     
    29972922           WRITE(lunout,*)'OK.'
    29982923#endif
     2924#endif
     2925
     2926c=============================================================
     2927c
     2928c Convertir les incrementations en tendances
     2929c
     2930      DO k = 1, klev
     2931      DO i = 1, klon
     2932         d_u(i,k) = ( u_seri(i,k) - u(i,k) ) / dtime
     2933         d_v(i,k) = ( v_seri(i,k) - v(i,k) ) / dtime
     2934         d_t(i,k) = ( t_seri(i,k)-t(i,k) ) / dtime
     2935         d_qx(i,k,ivap) = ( q_seri(i,k) - qx(i,k,ivap) ) / dtime
     2936         d_qx(i,k,iliq) = ( ql_seri(i,k) - qx(i,k,iliq) ) / dtime
     2937      ENDDO
     2938      ENDDO
     2939c
     2940      IF (nqmax.GE.3) THEN
     2941      DO iq = 3, nqmax
     2942      DO  k = 1, klev
     2943      DO  i = 1, klon
     2944         d_qx(i,k,iq) = ( tr_seri(i,k,iq-2) - qx(i,k,iq) ) / dtime
     2945      ENDDO
     2946      ENDDO
     2947      ENDDO
     2948      ENDIF
     2949c
     2950cIM rajout diagnostiques bilan KP pour analyse MJO par Jun-Ichi Yano
     2951c#include "write_bilKP_ins.h"
     2952c#include "write_bilKP_ave.h"
     2953c
     2954c Sauvegarder les valeurs de t et q a la fin de la physique:
     2955c
     2956      DO k = 1, klev
     2957      DO i = 1, klon
     2958         t_ancien(i,k) = t_seri(i,k)
     2959         q_ancien(i,k) = q_seri(i,k)
     2960      ENDDO
     2961      ENDDO
     2962c
     2963c 22.03.04 BEG
     2964c=============================================================
     2965c   Ecriture des sorties
     2966c=============================================================
     2967#ifdef CPP_IOIPSL
     2968
     2969#ifdef histhf
     2970#include "write_histhf.h"
     2971#endif
     2972
     2973#ifdef histday
     2974#include "write_histday.h"
     2975#include "write_histday_seri.h"
     2976#endif
     2977
     2978#ifdef histmth
     2979#include "write_histmth.h"
     2980#endif
     2981
     2982#ifdef histins
     2983#include "write_histins.h"
     2984#endif
     2985
     2986#ifdef histREGDYN
     2987#include "write_histREGDYN.h"
     2988#endif
     2989
     2990#ifdef histISCCP
     2991#include "write_histISCCP.h"
     2992#endif
     2993
     2994
     2995#ifdef histmthNMC
     2996#include "write_histmthNMC.h"
     2997#endif
     2998
    29992999#endif
    30003000
Note: See TracChangeset for help on using the changeset viewer.