Changeset 1412 for LMDZ4/trunk


Ignore:
Timestamp:
Jul 9, 2010, 3:06:40 PM (14 years ago)
Author:
idelkadi
Message:

Nettoyage dans physiq.F des anciennes versions des options iflag_cldcon=5 ou 6.

File:
1 edited

Legend:

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

    r1411 r1412  
    11! $Id$
    2 !
    32c#define IO_DEBUG
    43
     
    262261      CHARACTER*4 bb2
    263262      CHARACTER*2 bb3
    264 c
    265263
    266264      real twriteSTD(klon,nlevSTD,nfiles)
     
    865863      REAL rflag(klon)          ! flag fonctionnement de convect
    866864      INTEGER iflagctrl(klon)          ! flag fonctionnement de convect
     865
    867866c -- convect43:
    868867      INTEGER ntra              ! nb traceurs pour convect4.3
     
    12331232      call phys_state_var_init(read_climoz)
    12341233      call phys_output_var_init
    1235 
    12361234      print*, '================================================='
    12371235cIM for NMC files
     
    12531251c         pmflxr=0.
    12541252c         pmflxs=0.
     1253
    12551254        itau_con=0
    12561255        first=.false.
     
    14241423           ema_pcb(i)  = 0.
    14251424           ema_pct(i)  = 0.
    1426            ema_workcbmf(i) = 0.
     1425c          ema_workcbmf(i) = 0.
    14271426          ENDDO
    14281427cIM15/11/02 rajout initialisation ibas_con,itop_con cf. SB =>BEG
     
    14491448c================================================================================
    14501449
    1451          ENDIF
     1450         ENDIF !debut
    14521451
    14531452           DO i=1,klon
     
    21712170     .        d_t_con,d_q_con,d_u_con,d_v_con,d_tr,
    21722171     .        rain_con, snow_con, ibas_con, itop_con, sigd,
    2173      .        upwd,dnwd,dnwd0,
     2172     .        ema_cbmf,upwd,dnwd,dnwd0,
    21742173     .        Ma,mip,Vprecip,cape,cin,tvp,Tconv,iflagctrl,
    21752174     .        pbase,bbase,dtvpdt1,dtvpdq1,dplcldt,dplcldr,qcondc,wd,
     
    21902189
    21912190          ELSE ! ok_cvl
     2191
    21922192c MAF conema3 ne contient pas les traceurs
    21932193          CALL conema3 (dtime,
     
    22292229
    22302230          DO i = 1, klon
    2231             ema_pcb(i)  = pbase(i)
     2231            ema_pcb(i)  = paprs(i,ibas_con(i))
    22322232          ENDDO
    22332233          DO i = 1, klon
    2234 
    22352234! L'idicage de itop_con peut cacher un pb potentiel
    22362235! FH sous la dictee de JYG, CR
     
    22432242              endif
    22442243            endif
    2245           ENDDO
    2246           DO i = 1, klon
    2247             ema_cbmf(i) = ema_workcbmf(i)
    2248           ENDDO     
     2244          ENDDO     
    22492245      ELSE IF (iflag_con.eq.0) THEN
    22502246          write(lunout,*) 'On n appelle pas la convection'
     
    24632459c  ==============
    24642460
    2465 ! Dans le cas où on active les thermiques, on fait partir l'ajustement
     2461! Dans le cas où on active les thermiques, on fait partir l'ajustement
    24662462! a partir du sommet des thermiques.
    24672463! Dans le cas contraire, on demarre au niveau 1.
     
    25472543     s   (klon,klev,q_seri,zqsat,clwcon0th,ptconvth,ratqsc,rnebcon0th)
    25482544         if(prt_level.ge.9) print*,' CLOUDS_GNO OK'
    2549 
    2550       endif
     2545       
     2546       endif
    25512547
    25522548c   ratqs stables
     
    26112607
    26122608      if (iflag_cldcon.eq.1 .or.iflag_cldcon.eq.2
    2613      s    .or.iflag_cldcon.eq.4) then
     2609     s    .or.iflag_cldcon.ge.4) then
    26142610
    26152611! On ajoute une constante au ratqsc*2 pour tenir compte de
     
    26462642      endif
    26472643
    2648       print*,'PHSYIQ NUAGES4'
    26492644
    26502645c
     
    28422837      else
    28432838! Ancienne version
    2844          cldfra(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.)
    2845          cldliq(:,:)=cldliq(:,:)+rnebcon(:,:)*clwcon(:,:)
     2839      cldfra(:,:)=min(max(cldfra(:,:),rnebcon(:,:)),1.)
     2840      cldliq(:,:)=cldliq(:,:)+rnebcon(:,:)*clwcon(:,:)
    28462841      endif
    28472842
     
    33123307!     s        ref_liq,ref_ice
    33133308          call phys_cosp(itap,dtime,freq_cosp,
    3314      $                 ecrit_mth,ecrit_day,ecrit_hf,overlap,
    3315      $                   klon,klev,rlon,rlat,presnivs,
     3309     $                   ok_mensuelCOSP,ok_journeCOSP,ok_hfCOSP,
     3310     $                   ecrit_mth,ecrit_day,ecrit_hf,
     3311     $                   klon,klev,rlon,rlat,presnivs,overlap,
    33163312     $                   ref_liq,ref_ice,
    33173313     $                   pctsrf(:,is_ter)+pctsrf(:,is_lic),
     
    33223318     $                   pmflxr(:,1:klev),pmflxs(:,1:klev),
    33233319     $                   mr_ozone,cldtau, cldemi)
     3320
    33243321!     L          calipso2D,calipso3D,cfadlidar,parasolrefl,atb,betamol,
    33253322!     L          cfaddbze,clcalipso2,dbze,cltlidarradar,
     
    34693466      wwriteSTD(:,:,4)=wlevSTD(:,:)
    34703467c
     3468cIM initialisation 5eme fichier de sortie
    34713469cIM ajoute 5eme niveau 170310 BEG
    34723470      twriteSTD(:,:,5)=tlevSTD(:,:)
     
    35813579cIM global posePB#include "write_bilKP_ave.h"
    35823580c
     3581
    35833582c Sauvegarder les valeurs de t et q a la fin de la physique:
    35843583c
     
    36393638      DO k = 1, klev
    36403639      DO i = 1, klon
    3641         theta(i,k)=t(i,k)*(100000./pplay(i,k))**(RD/RCPD)
     3640cJYG/IM theta en debut du pas de temps
     3641cJYG/IM       theta(i,k)=t(i,k)*(100000./pplay(i,k))**(RD/RCPD)
     3642cJYG/IM theta en fin de pas de temps de physique
     3643        theta(i,k)=t_seri(i,k)*(100000./pplay(i,k))**(RD/RCPD)
    36423644      ENDDO
    36433645      ENDDO
Note: See TracChangeset for help on using the changeset viewer.