Index: /LMDZ5/trunk/libf/phy1d/1DUTILS.h =================================================================== --- /LMDZ5/trunk/libf/phy1d/1DUTILS.h (revision 1959) +++ /LMDZ5/trunk/libf/phy1d/1DUTILS.h (revision 1960) @@ -4,5 +4,5 @@ c c - SUBROUTINE conf_unicol( tapedef ) + SUBROUTINE conf_unicol c #ifdef CPP_IOIPSL @@ -15,10 +15,4 @@ c----------------------------------------------------------------------- c Auteurs : A. Lahellec . -c -c Arguments : -c -c tapedef : - - INTEGER tapedef c c Declarations : @@ -367,7 +361,4 @@ c Variables locales pour NetCDF: c ------------------------------ - INTEGER nid, nvarid - INTEGER idim_s - INTEGER ierr, ierr_file INTEGER iq INTEGER length @@ -378,6 +369,4 @@ character*80 abort_message LOGICAL found -c - INTEGER nb modname = 'dyn1deta0 : ' @@ -508,5 +497,4 @@ c ---------- CHARACTER*(*) fichnom - REAL time cAl1 plev tronque pour .nc mais plev(klev+1):=0 real :: plev(klon,klev),play (klon,klev),phi(klon,klev) @@ -520,7 +508,6 @@ c Variables locales pour NetCDF: c ------------------------------ - INTEGER nid, nvarid - INTEGER idim_s - INTEGER ierr, ierr_file + INTEGER nid + INTEGER ierr INTEGER iq,l INTEGER length @@ -535,9 +522,4 @@ DATA nb / 0 / - REAL zan0,zjulian,hours - INTEGER yyears0,jjour0, mmois0 - character*30 unites - -cDbg CALL open_restartphy(fichnom) print*,'redm1 ',fichnom,klon,klev,nqtot @@ -550,6 +532,6 @@ ierr = NF_OPEN(fichnom, NF_WRITE, nid) IF (ierr .NE. NF_NOERR) THEN - PRINT*, "Pb. d ouverture "//fichnom - CALL abort + abort_message="Pb. d ouverture "//fichnom + CALL abort_gcm('Modele 1D',abort_message,1) ENDIF @@ -661,7 +643,7 @@ ! traitement des point normaux DO j=2,jm-1 - ig=2+(j-2)*(im-1) + ig=2+(j-2)*(im-1) CALL SCOPY(im-1,pfi(ig,ifield),1,pdyn(1,j,ifield),1) - pdyn(im,j,ifield)=pdyn(1,j,ifield) + pdyn(im,j,ifield)=pdyn(1,j,ifield) ENDDO ENDDO @@ -992,5 +974,5 @@ DO ifield=1,nfield DO j=2,jm-1 - ig=2+(j-2)*(im-1) + ig=2+(j-2)*(im-1) CALL SCOPY(im-1,pdyn(1,j,ifield),1,pfi(ig,ifield),1) ENDDO @@ -1151,5 +1133,5 @@ !====================================================================== - SUBROUTINE read_tsurf1d(knon,knindex,sst_out) + SUBROUTINE read_tsurf1d(knon,sst_out) ! This subroutine specifies the surface temperature to be used in 1D simulations @@ -1158,5 +1140,4 @@ INTEGER, INTENT(IN) :: knon ! nomber of points on compressed grid - INTEGER, DIMENSION(klon), INTENT(IN) :: knindex ! grid point number for compressed grid REAL, DIMENSION(klon), INTENT(OUT) :: sst_out ! tsurf used to force the single-column model @@ -1220,6 +1201,5 @@ SUBROUTINE advect_va(llm,omega,d_t_va,d_q_va,d_u_va,d_v_va, - ! q,temp,u,v, - ! play,plev) + s q,temp,u,v,play) !itlmd !---------------------------------------------------------------------- @@ -1237,5 +1217,5 @@ real q(llm,3),temp(llm) real u(llm),v(llm) - real play(llm),plev(llm+1) + real play(llm) ! interne integer l @@ -1323,6 +1303,6 @@ real dph(llm),dqdp(llm),dtdp(llm) ! interne - integer l,k - real alpha,omdn,omup + integer k + real omdn,omup ! dudp=0. @@ -1403,6 +1383,6 @@ character*80 fich_toga - integer no,l,k,ip - real riy,rim,rid,rih,bid + integer k,ip + real bid integer iy,im,id,ih @@ -1422,5 +1402,5 @@ do k = 1, nlev_toga - read(21,230) plev_toga(k,ip), t_toga(k,ip), q_toga(k,ip) + read(21,230) plev_toga(k,ip), t_toga(k,ip), q_toga(k,ip) : ,u_toga(k,ip), v_toga(k,ip), w_toga(k,ip) : ,ht_toga(k,ip), vt_toga(k,ip), hq_toga(k,ip), vq_toga(k,ip) @@ -1443,6 +1423,4 @@ 223 format(4i3,6f8.2) - 226 format(f7.1,1x,10f8.2) - 227 format(f7.1,1x,1p,4e11.3) 230 format(6f9.3,4e11.3) @@ -1462,12 +1440,11 @@ character*80 fich_sandu - integer no,l,k,ip - real riy,rim,rid,rih,bid - + integer ip integer iy,im,id,ih - real plev_min - - plev_min = 55000. ! pas de tendance de vap. d eau au-dessus de 55 hPa + real plev_min + + print*,'nlev_sandu',nlev_sandu + plev_min = 55000. ! pas de tendance de vap. d eau au-dessus de 55 hPa open(21,file=trim(fich_sandu),form='formatted') @@ -1482,7 +1459,4 @@ 223 format(4i3,f8.2) - 226 format(f7.1,1x,10f8.2) - 227 format(f7.1,1x,1p,4e11.3) - 230 format(6f9.3,4e11.3) return @@ -1504,11 +1478,10 @@ character*80 fich_astex - integer no,l,k,ip - real riy,rim,rid,rih,bid - + integer ip integer iy,im,id,ih real plev_min + print*,'nlev_astex',nlev_astex plev_min = 55000. ! pas de tendance de vap. d eau au-dessus de 55 hPa @@ -1528,7 +1501,4 @@ 223 format(4i3,e13.2,f7.2,f7.3,f7.2,f7.3,f7.2) - 226 format(f7.1,1x,10f8.2) - 227 format(f7.1,1x,1p,4e11.3) - 230 format(6f9.3,4e11.3) return @@ -1551,5 +1521,5 @@ character*80 :: fich_twpice real*8 time(ntime) - real*8 lat, lon, alt, phis + real*8 lat, lon, alt, phis real*8 lev(nlevel) real*8 plev(nlevel,ntime) @@ -1561,5 +1531,5 @@ real*8 T_adv_h(nlevel,ntime) real*8 T_adv_v(nlevel,ntime), q_adv_h(nlevel,ntime) - real*8 q_adv_v(nlevel,ntime) + real*8 q_adv_v(nlevel,ntime) real*8 s(nlevel,ntime), s_adv_h(nlevel,ntime) real*8 s_adv_v(nlevel,ntime) @@ -1976,5 +1946,4 @@ integer ierr - integer i integer timevar,levvar integer timelen,levlen @@ -2050,6 +2019,6 @@ real omega_mod(llm),o3mmr_mod(llm) - integer l,k,k1,k2,kp - real aa,frac,frac1,frac2,fact + integer l,k,k1,k2 + real frac,frac1,frac2,fact do l = 1, llm @@ -2168,6 +2137,6 @@ real o3mmr_mod(llm),ql_mod(llm),qt_mod(llm) - integer l,k,k1,k2,kp - real aa,frac,frac1,frac2,fact + integer l,k,k1,k2 + real frac,frac1,frac2,fact do l = 1, llm @@ -2444,5 +2413,5 @@ real ts_prof ! local: - integer it_sandu1, it_sandu2,k + integer it_sandu1, it_sandu2 real timeit,time_sandu1,time_sandu2,frac ! Check that initial day of the simulation consistent with SANDU period: @@ -2511,8 +2480,9 @@ character*80 fich_armcu - integer no,l,k,ip - real riy,rim,rid,rih,bid + integer ip integer iy,im,id,ih,in + + print*,'nlev_armcu',nlev_armcu open(21,file=trim(fich_armcu),form='formatted') @@ -2529,7 +2499,4 @@ 223 format(5i3,5f8.3) - 226 format(f7.1,1x,10f8.2) - 227 format(f7.1,1x,1p,4e11.3) - 230 format(6f9.3,4e11.3) return @@ -2571,6 +2538,6 @@ real hq_mod(llm),vq_mod(llm) - integer l,k,k1,k2,kp - real aa,frac,frac1,frac2,fact + integer l,k,k1,k2 + real frac,frac1,frac2,fact do l = 1, llm @@ -2684,5 +2651,5 @@ real div_prof,ts_prof,ug_prof,vg_prof,ufa_prof,vfa_prof ! local: - integer it_astex1, it_astex2,k + integer it_astex1, it_astex2 real timeit,time_astex1,time_astex2,frac @@ -2985,9 +2952,5 @@ integer, parameter :: ilesfile=1 - integer :: ierr,irad,imax,jtot,k,itrac,nt1,nt2 - logical :: lmoist,lcoriol,ltimedep - real :: xsize,ysize - real :: ustin,wsvsurf,timerad - character(80) :: chmess + integer :: ierr,k,itrac,nt1,nt2 if(.not.(llesread)) return @@ -3043,5 +3006,5 @@ implicit none - integer nlev_max,kmax,kmax2 + integer nlev_max,kmax logical :: llesread = .true. @@ -3052,9 +3015,5 @@ integer, parameter :: ilesfile=1 - integer :: ierr,irad,imax,jtot,k - logical :: lmoist,lcoriol,ltimedep - real :: xsize,ysize - real :: ustin,wsvsurf,timerad - character(80) :: chmess + integer :: k,ierr if(.not.(llesread)) return @@ -3079,5 +3038,5 @@ implicit none - integer nlev_max,kmax,kmax2 + integer nlev_max,kmax logical :: llesread = .true. @@ -3088,9 +3047,5 @@ integer, parameter :: ilesfile=1 - integer :: ierr,irad,imax,jtot,k - logical :: lmoist,lcoriol,ltimedep - real :: xsize,ysize - real :: ustin,wsvsurf,timerad - character(80) :: chmess + integer :: ierr,k if(.not.(llesread)) return @@ -3117,5 +3072,5 @@ implicit none - integer nlev_max,kmax,kmax2 + integer nlev_max,kmax logical :: llesread = .true. @@ -3127,9 +3082,5 @@ integer, parameter :: ilesfile=1 integer, parameter :: ifile=2 - integer :: ierr,irad,imax,jtot,k - logical :: lmoist,lcoriol,ltimedep - real :: xsize,ysize - real :: ustin,wsvsurf,timerad - character(80) :: chmess + integer :: ierr,jtot,k if(.not.(llesread)) return @@ -3182,8 +3133,6 @@ integer ntime,nlevel - integer l,k character*80 :: fich_amma - real*8 time(ntime) - real*8 zz(nlevel) + real*8 zz(nlevel) real*8 temp(nlevel),pp(nlevel) @@ -3192,5 +3141,5 @@ real*8 dw(nlevel,ntime) real*8 dt(nlevel,ntime) - real*8 dq(nlevel,ntime) + real*8 dq(nlevel,ntime) real*8 flat(ntime),sens(ntime) @@ -3522,8 +3471,6 @@ integer ntime,nlevel - integer l,k character*80 :: fich_fire - real*8 time(ntime) - real*8 zz(nlevel) + real*8 zz(nlevel) real*8 thl(nlevel) @@ -3532,5 +3479,5 @@ real*8 ug(nlevel,ntime),vg(nlevel,ntime),wls(nlevel,ntime) real*8 dqtdx(nlevel,ntime),dqtdy(nlevel,ntime) - real*8 dqtdt(nlevel,ntime),thl_rad(nlevel,ntime) + real*8 dqtdt(nlevel,ntime),thl_rad(nlevel,ntime) integer nid, ierr Index: /LMDZ5/trunk/libf/phy1d/1D_decl_cases.h =================================================================== --- /LMDZ5/trunk/libf/phy1d/1D_decl_cases.h (revision 1959) +++ /LMDZ5/trunk/libf/phy1d/1D_decl_cases.h (revision 1960) @@ -16,8 +16,4 @@ real sec_print !! - integer nn - integer it_toga1, it_toga2 - real time_toga1,time_toga2 - real ts_toga(nt_toga) real plev_toga(nlev_toga,nt_toga),w_toga(nlev_toga,nt_toga) @@ -34,5 +30,5 @@ real hq_prof(nlev_toga),vq_prof(nlev_toga) - real plev_mod(llm),w_mod(llm), t_mod(llm),q_mod(llm) + real w_mod(llm), t_mod(llm),q_mod(llm) real u_mod(llm),v_mod(llm), ht_mod(llm),vt_mod(llm) real hq_mod(llm),vq_mod(llm),qv_mod(llm),ql_mod(llm),qt_mod(llm) @@ -87,5 +83,4 @@ character*80 :: fich_amma ! Option du cas AMMA ou on impose la discretisation verticale (Ap,Bp) - logical :: fixe_disvert=.true. integer nlev_amma, nt_amma ! parameter (nlev_amma=29, nt_amma=48) ! Fleur, juillet 2012 @@ -104,6 +99,4 @@ !profils initiaux: real plev_amma(nlev_amma) - real tv_amma(nlev_amma),rho_amma(nlev_amma) - real thv_amma(nlev_amma) real z_amma(nlev_amma) @@ -111,6 +104,4 @@ real u_amma(nlev_amma) real v_amma(nlev_amma) - - real thvsurf_amma,tvsurf_amma,rhosurf_amma,thsurf real th_ammai(nlev_amma),q_ammai(nlev_amma) @@ -130,5 +121,5 @@ !champs interpoles - real plev_profamma(nlev_amma),vitw_profamma(nlev_amma) + real vitw_profamma(nlev_amma) real ht_profamma(nlev_amma) real hq_profamma(nlev_amma) @@ -148,5 +139,4 @@ integer year_ini_fire, day_ini_fire, mth_ini_fire real heure_ini_fire - real day_ju_ini_fire ! Julian day of fire first day parameter (year_ini_fire=1987) parameter (mth_ini_fire=7) @@ -154,19 +144,4 @@ parameter (heure_ini_fire=0.) !0h en secondes -!profils initiaux: - real z_fire(nlev_fire) - real thl_fire(nlev_fire),qt_fire(nlev_fire) - real u_fire(nlev_fire), v_fire(nlev_fire) - real tke_fire(nlev_fire) - -!forcings - real ugeo_fire(nlev_fire),vgeo_fire(nlev_fire) - real wls_fire(nlev_fire),dqtdx_fire(nlev_fire) - real dqtdy_fire(nlev_fire) - real dqtdt_fire(nlev_fire),thl_rad_fire(nlev_fire) - - real ugeo_mod(llm),vgeo_mod(llm),wls_mod(llm) - real dqtdx_mod(llm),dqtdy_mod(llm),dqtdt_mod(llm) - real thl_rad_mod(llm) !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! Declarations specifiques au cas GCSSold @@ -180,10 +155,8 @@ real Ts_gcssold real dtime_frcg -cAl1 logical :: imp_fcg_gcssold -c logical :: ts_fcg_gcssold -c logical :: Tp_fcg_gcssold logical :: Turb_fcg_gcssold - common /turb_forcing/ dtime_frcg, - $ Turb_fcg_gcssold, hthturb_gcssold, hqturb_gcssold + + common /turb_forcing/ + s dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! Declarations specifiques au cas Arm_cu @@ -206,8 +179,5 @@ real adv_qt_armcu(nt_armcu) real theta_mod(llm),rv_mod(llm),play_mod(llm) - real d_t_dyn_ls(llm),d_q_dyn_ls(llm) ! profc comme "profil armcu" - real h_profc,play_profc,t_profc,th_profc,plev_profc - real u_profc,v_profc,qv_profc,rv_profc ! forcages interpoles dans le temps @@ -226,12 +196,9 @@ parameter (mth_ini_sandu=7) parameter (day_ini_sandu=15) ! 196 = 15 juillet 2006 - real dt_sandu, tau_sandu + real dt_sandu logical :: trouve_700=.true. parameter (dt_sandu=6.*3600.) ! forcages donnes ttes les 6 heures par ifa_sandu.txt ! parameter (tau_sandu=3600.) ! temps de relaxation u,v,thetal,qt vers profil init et au dessus 700hPa !! - integer it_sandu1, it_sandu2 - real time_sandu1,time_sandu2 - real ts_sandu(nt_sandu) ! profs comme "profil sandu" @@ -243,6 +210,4 @@ real thl_mod(llm),omega_mod(llm),o3mmr_mod(llm),tke_mod(llm) -! pour relaxer u,v,thl et qt vers les profils initiaux au dessus de 700hPa - real relax_u(llm),relax_v(llm),relax_thl(llm),relax_q(llm,2) !vertical advection computation real d_t_z(llm), d_q_z(llm) @@ -260,8 +225,6 @@ parameter (mth_ini_astex=6) parameter (day_ini_astex=13) ! 165 = 13 juin 1992 - real dt_astex, tau_astex + real dt_astex parameter (dt_astex=3600.) ! forcages donnes ttes les heures par ifa_astex.txt - integer it_astex1, it_astex2 - real time_astex1,time_astex2 real ts_astex(nt_astex),div_astex(nt_astex),ug_astex(nt_astex) real vg_astex(nt_astex),ufa_astex(nt_astex),vfa_astex(nt_astex) Index: /LMDZ5/trunk/libf/phy1d/1Dconv.h =================================================================== --- /LMDZ5/trunk/libf/phy1d/1Dconv.h (revision 1959) +++ /LMDZ5/trunk/libf/phy1d/1Dconv.h (revision 1960) @@ -22,11 +22,10 @@ REAL hplaym(100) !pression en hPa milieux des couches Meso-NH - integer i,j,k,ii,ll,in - REAL tsol,qsol + integer i,j,k,ll,in CHARACTER*80 file_forctl,file_fordat - COMMON/com1_phys_gcss/klev,play,JM,coef1,coef2 - COMMON/com2_phys_gcss/nblvlm,playm,hplaym + COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev + COMMON/com2_phys_gcss/playm,hplaym,nblvlm c====================================================================== @@ -170,5 +169,5 @@ hqTurbbef(i)=hqTurbaft(i) enddo - tsbef = tsaft + tsbef = tsaft timebef=pasprev*dt timeaft=timebef+dt @@ -213,5 +212,5 @@ print*,'hqTurb_mes ',(hqTurb_mes(i),i=1,nblvlm) endif - IF (ts_fcg) print*,'ts_subr', ts_subr + IF (ts_fcg) print*,'ts_subr', ts_subr c*** on interpole les champs meso_NH sur les niveaux de pression*** c*** gcm . on obtient le nouveau champ after *** @@ -263,5 +262,5 @@ hqTurb(ll)=hqTurbaft(ll) enddo - ts_subr = tsaft + ts_subr = tsaft else ! temps.ge.pasmax c*** on interpole sur les pas de temps de 10mn du gcm a partir *** @@ -282,5 +281,5 @@ endif ! Turb_fcg enddo - ts_subr = ((timeaft-time)*tsbef + (time-timebef)*tsaft)/dt + ts_subr = ((timeaft-time)*tsbef + (time-timebef)*tsaft)/dt endif ! temps.ge.pasmax c @@ -440,5 +439,5 @@ endif ! Turb_fcg enddo - ts_subr = tsaft + ts_subr = tsaft close(99) close(98) @@ -505,8 +504,8 @@ REAL hplaym(100)!pression en hecto-Pa des milieux de couche Meso-NH - COMMON/com1_phys_gcss/klev,play,JM,coef1,coef2 - COMMON/com2_phys_gcss/nblvlm,playm,hplaym - - integer i,k,klevgcm + COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev + COMMON/com2_phys_gcss/playm,hplaym,nblvlm + + integer k,klevgcm real playgcm(klevgcm) ! pression en milieu de couche du gcm real psolgcm @@ -577,7 +576,7 @@ REAL playm(100) !pression en Pa milieu de chaque couche Meso-NH REAL hplaym(100) !pression en hPa des milieux de couche Meso-NH - COMMON/com2_phys_gcss/nblvlm,playm,hplaym - - INTEGER i,lu,k,mlz,mlzh,j + COMMON/com2_phys_gcss/playm,hplaym,nblvlm + + INTEGER i,lu,mlz,mlzh character*80 file_forctl @@ -644,5 +643,5 @@ real ts c - INTEGER i, k + INTEGER k c LOGICAL imp_fcg,ts_fcg,Turb_fcg @@ -725,10 +724,10 @@ REAL hplaym(100)!pression en hPa des milieux de couche Meso-NH - COMMON/com1_phys_gcss/klev,play,JM,coef1,coef2 - COMMON/com2_phys_gcss/nblvlm,playm,hplaym + COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev + COMMON/com2_phys_gcss/playm,hplaym,nblvlm REAL psol REAL val - INTEGER k, mlz, mlzh + INTEGER k, mlz Index: /LMDZ5/trunk/libf/phy1d/add_pbl_tend.F =================================================================== --- /LMDZ5/trunk/libf/phy1d/add_pbl_tend.F (revision 1959) +++ /LMDZ5/trunk/libf/phy1d/add_pbl_tend.F (revision 1960) @@ -19,6 +19,7 @@ real dtime_frcg logical Turb_fcg_gcssold - common /turb_forcing/ dtime_frcg, - $ Turb_fcg_gcssold, hthturb_gcssold, hqturb_gcssold + common /turb_forcing/ + s dtime_frcg,hthturb_gcssold, hqturb_gcssold,Turb_fcg_gcssold + ! Arguments : !------------ Index: /LMDZ5/trunk/libf/phy1d/compar1d.h =================================================================== --- /LMDZ5/trunk/libf/phy1d/compar1d.h (revision 1959) +++ /LMDZ5/trunk/libf/phy1d/compar1d.h (revision 1960) @@ -27,7 +27,9 @@ logical :: ok_old_disvert - common/com_par1d/forcing_type,nat_surf,tsurf,rugos, & + common/com_par1d/ + & nat_surf,tsurf,rugos, & & qsol,qsurf,psurf,zsurf,albedo,time,time_ini,xlat,xlon,airefi, & & wtsurf,wqsurf,restart_runoff,xagesno,qsolinp,zpicinp, & + & forcing_type, & restart,ok_old_disvert Index: /LMDZ5/trunk/libf/phy1d/lmdz1d.F =================================================================== --- /LMDZ5/trunk/libf/phy1d/lmdz1d.F (revision 1959) +++ /LMDZ5/trunk/libf/phy1d/lmdz1d.F (revision 1960) @@ -81,8 +81,4 @@ integer :: an -! - real :: paire = 1. ! aire de la maille -!** common /flux_arp/fsens,flat,ok_flux_surf - !--------------------------------------------------------------------- ! Declarations related to forcing and initial profiles @@ -92,5 +88,5 @@ integer llm700,nq1,nq2 INTEGER, PARAMETER :: nlev_max=1000, nqmx=1000 - real timestep, frac, timeit + real timestep, frac real height(nlev_max),tttprof(nlev_max),qtprof(nlev_max), . uprof(nlev_max),vprof(nlev_max),e12prof(nlev_max), @@ -100,5 +96,4 @@ . qprof(nlev_max,nqmx) - real :: fff c integer :: forcing_type logical :: forcing_les = .false. @@ -138,5 +133,4 @@ real :: pzero=1.e5 real :: play (llm),zlay (llm),sig_s(llm),plev(llm+1) - real :: playd(llm),zlayd(llm),ap_amma(llm+1),bp_amma(llm+1),poub !--------------------------------------------------------------------- @@ -144,5 +138,4 @@ !--------------------------------------------------------------------- - integer :: iq real :: phi(llm) real :: teta(llm),tetal(llm),temp(llm),u(llm),v(llm),w(llm) @@ -152,9 +145,8 @@ real :: sfdt, cfdt real :: du_phys(llm),dv_phys(llm),dt_phys(llm) - real :: du_dyn(llm),dv_dyn(llm),dt_dyn(llm) - real :: dt_cooling(llm),d_t_cool(llm),d_th_adv(llm) - real :: dq_cooling(llm),d_q_cool(llm) - real :: tmpvar(llm) + real :: dt_dyn(llm) + real :: dt_cooling(llm),d_th_adv(llm) real :: alpha + real :: ttt REAL, ALLOCATABLE, DIMENSION(:,:):: q @@ -205,6 +197,5 @@ ! Fichiers et d'autres variables !--------------------------------------------------------------------- - real ttt,bow,q1 - integer :: ierr,k,l,i,it=1,mxcalc + integer :: k,l,i,it=1,mxcalc integer jjmp1 parameter (jjmp1=jjm+1-1/jjm) @@ -231,5 +222,5 @@ !--------------------------------------------------------------------- cAl1 - call conf_unicol(99) + call conf_unicol cAl1 moves this gcssold var from common fcg_gcssold to Turb_fcg_gcssold = xTurb_fcg_gcssold @@ -375,5 +366,5 @@ itau_phy = 0 call ymds2ju(annee_ref,mois,day_ref,heure,day) - day_ini = day + day_ini = int(day) day_end = day_ini + fnday @@ -859,11 +850,11 @@ ! du_age(1:mxcalc)= -2.*sfdt/timestep* - : (sfdt*(u(1:mxcalc)-ug(1:mxcalc)) - - : cfdt*(v(1:mxcalc)-vg(1:mxcalc)) ) + s (sfdt*(u(1:mxcalc)-ug(1:mxcalc)) - + s cfdt*(v(1:mxcalc)-vg(1:mxcalc)) ) !! : fcoriolis*(v(1:mxcalc)-vg(1:mxcalc)) ! dv_age(1:mxcalc)= -2.*sfdt/timestep* - : (cfdt*(u(1:mxcalc)-ug(1:mxcalc)) + - : sfdt*(v(1:mxcalc)-vg(1:mxcalc)) ) + s (cfdt*(u(1:mxcalc)-ug(1:mxcalc)) + + s sfdt*(v(1:mxcalc)-vg(1:mxcalc)) ) !! : -fcoriolis*(u(1:mxcalc)-ug(1:mxcalc)) ! @@ -879,12 +870,12 @@ !!!!!!!!!!!!!!!!!!!!!!!!!!!!! u(1:mxcalc)=u(1:mxcalc) + timestep*( - : du_phys(1:mxcalc) - : +du_age(1:mxcalc) ) + s du_phys(1:mxcalc) + s +du_age(1:mxcalc) ) v(1:mxcalc)=v(1:mxcalc) + timestep*( - : dv_phys(1:mxcalc) - : +dv_age(1:mxcalc) ) + s dv_phys(1:mxcalc) + s +dv_age(1:mxcalc) ) q(1:mxcalc,:)=q(1:mxcalc,:)+timestep*( - : dq(1:mxcalc,:) - : +d_q_adv(1:mxcalc,:) ) + s dq(1:mxcalc,:) + s +d_q_adv(1:mxcalc,:) ) if (prt_level.ge.1) then