Changeset 2221 for LMDZ5/trunk
- Timestamp:
- Mar 9, 2015, 7:38:03 AM (10 years ago)
- Location:
- LMDZ5/trunk/libf
- Files:
-
- 1 deleted
- 21 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/trunk/libf/dyn3d/calfis.F
r2037 r2221 21 21 $ pdq, 22 22 $ flxw, 23 $ clesphy0,24 23 $ pdufi, 25 24 $ pdvfi, … … 131 130 REAL,INTENT(OUT) :: pdpsfi(iip1,jjp1) ! tendency on surface pressure (Pa/s) 132 131 133 INTEGER,PARAMETER :: longcles = 20134 REAL,INTENT(IN) :: clesphy0( longcles ) ! unused135 136 132 137 133 c Local variables : … … 465 461 . zphis, 466 462 . presnivs, 467 . clesphy0,468 463 . zufi, 469 464 . zvfi, -
LMDZ5/trunk/libf/dyn3d/ce0l.F90
r1984 r2221 41 41 #include "temps.h" 42 42 #include "logic.h" 43 INTEGER, PARAMETER :: longcles=2044 REAL, DIMENSION(longcles) :: clesphy045 43 REAL, DIMENSION(iip1,jjp1) :: masque 46 44 CHARACTER(LEN=15) :: calnd 47 45 REAL, DIMENSION(iip1,jjp1) :: phis ! geopotentiel au sol 48 46 !------------------------------------------------------------------------------- 49 CALL conf_gcm( 99, .TRUE. , clesphy0)47 CALL conf_gcm( 99, .TRUE. ) 50 48 51 49 CALL Init_Phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/)) -
LMDZ5/trunk/libf/dyn3d/conf_gcm.F90
r2218 r2221 2 2 ! $Id$ 3 3 4 SUBROUTINE conf_gcm( tapedef, etatinit , clesphy0)4 SUBROUTINE conf_gcm( tapedef, etatinit ) 5 5 6 6 USE control_mod … … 23 23 ! etatinit : = TRUE , on ne compare pas les valeurs des para- 24 24 ! -metres du zoom avec celles lues sur le fichier start . 25 ! clesphy0 : sortie . 26 27 LOGICAL etatinit 28 INTEGER tapedef 29 30 INTEGER longcles 31 PARAMETER( longcles = 20 ) 32 REAL clesphy0( longcles ) 25 26 LOGICAL,INTENT(IN) :: etatinit 27 INTEGER,INTENT(IN) :: tapedef 33 28 34 29 ! Declarations : … … 41 36 include "temps.h" 42 37 include "comconst.h" 43 44 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique45 ! include "clesphys.h"46 38 include "iniprint.h" 47 39 … … 381 373 ip_ebil_dyn = 0 382 374 CALL getin('ip_ebil_dyn',ip_ebil_dyn) 383 384 DO i = 1, longcles385 clesphy0(i) = 0.386 ENDDO387 375 388 376 !cc .... P. Le Van , ajout le 7/03/95 .pour le zoom ... -
LMDZ5/trunk/libf/dyn3d/gcm.F
r2151 r2221 90 90 !#include "indicesol.h" 91 91 #endif 92 INTEGER longcles93 PARAMETER ( longcles = 20 )94 REAL clesphy0( longcles )95 SAVE clesphy096 97 98 92 99 93 REAL zdtvr … … 175 169 c --------------------------------------- 176 170 c 177 ! Ehouarn: dump possibility of using defrun 178 !#ifdef CPP_IOIPSL 179 CALL conf_gcm( 99, .TRUE. , clesphy0 ) 171 CALL conf_gcm( 99, .TRUE.) 180 172 if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm", 181 173 s "iphysiq must be a multiple of iperiod", 1) 182 !#else183 ! CALL defrun( 99, .TRUE. , clesphy0 )184 !#endif185 174 186 175 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 530 519 531 520 532 CALL leapfrog(ucov,vcov,teta,ps,masse,phis,q,clesphy0, 533 . time_0) 521 CALL leapfrog(ucov,vcov,teta,ps,masse,phis,q,time_0) 534 522 535 523 END -
LMDZ5/trunk/libf/dyn3d/leapfrog.F
r2039 r2221 4 4 c 5 5 c 6 SUBROUTINE leapfrog(ucov,vcov,teta,ps,masse,phis,q,clesphy0, 7 & time_0) 6 SUBROUTINE leapfrog(ucov,vcov,teta,ps,masse,phis,q,time_0) 8 7 9 8 … … 70 69 #include "academic.h" 71 70 72 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique73 ! #include "clesphys.h"74 75 INTEGER,PARAMETER :: longcles = 2076 REAL,INTENT(IN) :: clesphy0( longcles ) ! not used77 71 REAL,INTENT(IN) :: time_0 ! not used 78 72 … … 446 440 $ ucov,vcov,teta,q,masse,ps,p,pk,phis,phi , 447 441 $ du,dv,dteta,dq, 448 $ flxw, 449 $ clesphy0, dufi,dvfi,dtetafi,dqfi,dpfi ) 442 $ flxw,dufi,dvfi,dtetafi,dqfi,dpfi ) 450 443 451 444 c ajout des tendances physiques: -
LMDZ5/trunk/libf/dyn3d_common/grilles_gcm_netcdf_sub.F90
r1984 r2221 31 31 INTEGER out_latudim,out_latvdim,out_dim(3) 32 32 INTEGER out_levdim 33 34 INTEGER, PARAMETER :: longcles = 2035 REAL clesphy0(longcles)36 33 37 34 INTEGER start(4),COUNT(4) … … 60 57 pa= 50000. 61 58 62 CALL conf_gcm( 99, .TRUE. , clesphy0)59 CALL conf_gcm( 99, .TRUE. ) 63 60 CALL iniconst 64 61 CALL inigeom -
LMDZ5/trunk/libf/dyn3dmem/calfis_loc.F
r2037 r2221 21 21 $ pdq, 22 22 $ flxw, 23 $ clesphy0,24 23 $ pdufi, 25 24 $ pdvfi, … … 141 140 REAL,INTENT(OUT) :: pdqfi(iip1,jjb_u:jje_u,llm,nqtot) ! tendency on tracers 142 141 REAL,INTENT(OUT) :: pdpsfi(iip1,jjb_u:jje_u) ! tendency on surface pressure (Pa/s) 143 144 INTEGER,PARAMETER :: longcles = 20145 REAL,INTENT(IN) :: clesphy0( longcles ) ! unused146 147 142 148 143 #ifdef CPP_PHYS … … 674 669 . zphis_omp, 675 670 . presnivs_omp, 676 . clesphy0,677 671 . zufi_omp, 678 672 . zvfi_omp, 679 673 . ztfi_omp, 680 674 . zqfi_omp, 681 c#ifdef INCA682 675 . flxwfi_omp, 683 c#endif684 676 . zdufi_omp, 685 677 . zdvfi_omp, -
LMDZ5/trunk/libf/dyn3dmem/call_calfis_mod.F90
r2021 r2221 68 68 69 69 70 SUBROUTINE call_calfis(itau,lafin, clesphy0,ucov_dyn,vcov_dyn,teta_dyn,masse_dyn,ps_dyn, &70 SUBROUTINE call_calfis(itau,lafin,ucov_dyn,vcov_dyn,teta_dyn,masse_dyn,ps_dyn, & 71 71 phis_dyn,q_dyn,flxw_dyn) 72 72 USE dimensions_mod … … 91 91 INTEGER,INTENT(IN) :: itau ! (time) iteration step number 92 92 LOGICAL,INTENT(IN) :: lafin ! .true. if final time step 93 REAL,INTENT(IN) :: clesphy0( : ) ! not used94 93 REAL,INTENT(INOUT) :: ucov_dyn(ijb_u:ije_u,llm) ! covariant zonal wind 95 94 REAL,INTENT(INOUT) :: vcov_dyn(ijb_v:ije_v,llm) ! covariant meridional wind … … 231 230 ucov,vcov,teta,q,masse,ps,p,pk,phis,phi , & 232 231 du,dv,dteta,dq, & 233 flxw, & 234 clesphy0, dufi,dvfi,dtetafi,dqfi,dpfi ) 232 flxw, dufi,dvfi,dtetafi,dqfi,dpfi ) 235 233 236 234 ijb=ij_begin -
LMDZ5/trunk/libf/dyn3dmem/ce0l.F90
r1984 r2221 48 48 #endif 49 49 50 INTEGER, PARAMETER :: longcles=2051 50 INTEGER :: ierr 52 REAL, DIMENSION(longcles) :: clesphy053 51 REAL, DIMENSION(iip1,jjp1) :: masque 54 52 CHARACTER(LEN=15) :: calnd 55 53 REAL, DIMENSION(iip1,jjp1) :: phis ! geopotentiel au sol 56 54 !------------------------------------------------------------------------------- 57 CALL conf_gcm( 99, .TRUE. , clesphy0)55 CALL conf_gcm( 99, .TRUE. ) 58 56 59 57 #ifdef CPP_MPI -
LMDZ5/trunk/libf/dyn3dmem/conf_gcm.F90
r2218 r2221 2 2 ! $Id$ 3 3 4 SUBROUTINE conf_gcm( tapedef, etatinit , clesphy0)4 SUBROUTINE conf_gcm( tapedef, etatinit ) 5 5 6 6 USE control_mod … … 27 27 ! etatinit : = TRUE , on ne compare pas les valeurs des para- 28 28 ! -metres du zoom avec celles lues sur le fichier start . 29 ! clesphy0 : sortie . 30 31 LOGICAL etatinit 32 INTEGER tapedef 33 34 INTEGER longcles 35 PARAMETER( longcles = 20 ) 36 REAL clesphy0( longcles ) 29 30 LOGICAL,INTENT(IN) :: etatinit 31 INTEGER,INTENT(IN) :: tapedef 37 32 38 33 ! Declarations : … … 45 40 include "temps.h" 46 41 include "comconst.h" 47 48 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique49 ! include "clesphys.h"50 42 include "iniprint.h" 51 43 -
LMDZ5/trunk/libf/dyn3dmem/gcm.F
r2180 r2221 83 83 #endif 84 84 85 INTEGER longcles86 PARAMETER ( longcles = 20 )87 REAL clesphy0( longcles )88 SAVE clesphy089 90 91 92 85 REAL zdtvr 93 86 … … 172 165 c 173 166 ! Ehouarn: dump possibility of using defrun 174 !#ifdef CPP_IOIPSL 175 CALL conf_gcm( 99, .TRUE. , clesphy0 ) 167 CALL conf_gcm( 99, .TRUE. ) 176 168 if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm", 177 169 s "iphysiq must be a multiple of iperiod", 1) 178 !#else179 ! CALL defrun( 99, .TRUE. , clesphy0 )180 !#endif181 170 c 182 171 c … … 547 536 548 537 c$OMP PARALLEL DEFAULT(SHARED) COPYIN(/temps/,/logici/,/logicl/) 549 CALL leapfrog_loc(ucov,vcov,teta,ps,masse,phis,q,clesphy0, 550 . time_0) 538 CALL leapfrog_loc(ucov,vcov,teta,ps,masse,phis,q,time_0) 551 539 c$OMP END PARALLEL 552 540 -
LMDZ5/trunk/libf/dyn3dmem/leapfrog_loc.F
r2185 r2221 9 9 10 10 SUBROUTINE leapfrog_loc(ucov0,vcov0,teta0,ps0, 11 & masse0,phis0,q0,clesphy0, 12 & time_0) 11 & masse0,phis0,q0,time_0) 13 12 14 13 USE misc_mod … … 82 81 ! include "mpif.h" 83 82 84 INTEGER,PARAMETER :: longcles = 2085 REAL,INTENT(IN) :: clesphy0( longcles ) ! not used86 83 REAL,INTENT(IN) :: time_0 ! not used 87 84 … … 757 754 IF( apphys ) THEN 758 755 759 CALL call_calfis(itau,lafin, clesphy0,ucov,vcov,teta,masse,ps,756 CALL call_calfis(itau,lafin,ucov,vcov,teta,masse,ps, 760 757 & phis,q,flxw) 761 758 ! #ifdef DEBUG_IO … … 882 879 ! $ du,dv,dteta,dq, 883 880 ! $ flxw, 884 ! $ clesphy0,dufi,dvfi,dtetafi,dqfi,dpfi )881 ! $ dufi,dvfi,dtetafi,dqfi,dpfi ) 885 882 ! ! CALL FTRACE_REGION_END("calfis") 886 883 ! ! ijb=ij_begin -
LMDZ5/trunk/libf/dyn3dpar/calfis_p.F
r2037 r2221 21 21 $ pdq, 22 22 $ flxw, 23 $ clesphy0,24 23 $ pdufi, 25 24 $ pdvfi, … … 140 139 REAL,INTENT(OUT) :: pdqfi(iip1,jjp1,llm,nqtot) ! tendency on tracers 141 140 REAL,INTENT(OUT) :: pdpsfi(iip1,jjp1) ! tendency on surface pressure (Pa/s) 142 143 INTEGER,PARAMETER :: longcles = 20144 REAL,INTENT(IN) :: clesphy0( longcles ) ! unused145 141 146 142 #ifdef CPP_PHYS … … 636 632 . zphis_omp, 637 633 . presnivs_omp, 638 . clesphy0,639 634 . zufi_omp, 640 635 . zvfi_omp, 641 636 . ztfi_omp, 642 637 . zqfi_omp, 643 c#ifdef INCA644 638 . flxwfi_omp, 645 c#endif646 639 . zdufi_omp, 647 640 . zdvfi_omp, -
LMDZ5/trunk/libf/dyn3dpar/ce0l.F90
r1984 r2221 48 48 #endif 49 49 50 INTEGER, PARAMETER :: longcles=2051 50 INTEGER :: ierr 52 REAL, DIMENSION(longcles) :: clesphy053 51 REAL, DIMENSION(iip1,jjp1) :: masque 54 52 CHARACTER(LEN=15) :: calnd 55 53 REAL, DIMENSION(iip1,jjp1) :: phis ! geopotentiel au sol 56 54 !------------------------------------------------------------------------------- 57 CALL conf_gcm( 99, .TRUE. , clesphy0)55 CALL conf_gcm( 99, .TRUE. ) 58 56 59 57 #ifdef CPP_MPI -
LMDZ5/trunk/libf/dyn3dpar/conf_gcm.F90
r2151 r2221 2 2 ! $Id$ 3 3 4 SUBROUTINE conf_gcm( tapedef, etatinit , clesphy0)4 SUBROUTINE conf_gcm( tapedef, etatinit ) 5 5 6 6 USE control_mod … … 26 26 ! etatinit : = TRUE , on ne compare pas les valeurs des para- 27 27 ! -metres du zoom avec celles lues sur le fichier start . 28 ! clesphy0 : sortie . 29 30 LOGICAL etatinit 31 INTEGER tapedef 32 33 INTEGER longcles 34 PARAMETER( longcles = 20 ) 35 REAL clesphy0( longcles ) 28 29 LOGICAL,INTENT(IN) :: etatinit 30 INTEGER,INTENT(IN) :: tapedef 36 31 37 32 ! Declarations : … … 44 39 include "temps.h" 45 40 include "comconst.h" 46 47 ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique48 ! include "clesphys.h"49 41 include "iniprint.h" 50 42 -
LMDZ5/trunk/libf/dyn3dpar/gcm.F
r2180 r2221 84 84 #endif 85 85 86 INTEGER longcles87 PARAMETER ( longcles = 20 )88 REAL clesphy0( longcles )89 SAVE clesphy090 91 92 93 86 REAL zdtvr 94 87 … … 172 165 c --------------------------------------- 173 166 c 174 ! Ehouarn: dump possibility of using defrun 175 !#ifdef CPP_IOIPSL 176 CALL conf_gcm( 99, .TRUE. , clesphy0 ) 167 CALL conf_gcm( 99, .TRUE. ) 177 168 if (mod(iphysiq, iperiod) /= 0) call abort_gcm("conf_gcm", 178 169 s "iphysiq must be a multiple of iperiod", 1) 179 !#else180 ! CALL defrun( 99, .TRUE. , clesphy0 )181 !#endif182 170 c 183 171 c … … 550 538 551 539 c$OMP PARALLEL DEFAULT(SHARED) COPYIN(/temps/,/logici/,/logicl/) 552 CALL leapfrog_p(ucov,vcov,teta,ps,masse,phis,q,clesphy0, 553 . time_0) 540 CALL leapfrog_p(ucov,vcov,teta,ps,masse,phis,q,time_0) 554 541 c$OMP END PARALLEL 555 542 -
LMDZ5/trunk/libf/dyn3dpar/leapfrog_p.F
r2180 r2221 5 5 c 6 6 7 SUBROUTINE leapfrog_p(ucov,vcov,teta,ps,masse,phis,q,clesphy0, 8 & time_0) 7 SUBROUTINE leapfrog_p(ucov,vcov,teta,ps,masse,phis,q,time_0) 9 8 10 9 use exner_hyb_m, only: exner_hyb … … 77 76 #include "academic.h" 78 77 79 INTEGER,PARAMETER :: longcles = 2080 REAL,INTENT(IN) :: clesphy0( longcles ) ! not used81 78 REAL,INTENT(IN) :: time_0 ! not used 82 79 … … 831 828 $ ucov,vcov,teta,q,masse,ps,p,pk,phis,phi , 832 829 $ du,dv,dteta,dq, 833 $ flxw, 834 $ clesphy0, dufi,dvfi,dtetafi,dqfi,dpfi ) 830 $ flxw, dufi,dvfi,dtetafi,dqfi,dpfi ) 835 831 ! CALL FTRACE_REGION_END("calfis") 836 832 ijb=ij_begin -
LMDZ5/trunk/libf/phydev/physiq.F90
r2097 r2221 4 4 SUBROUTINE physiq (nlon,nlev, & 5 5 & debut,lafin,jD_cur, jH_cur,pdtphys, & 6 & paprs,pplay,pphi,pphis,presnivs, clesphy0,&6 & paprs,pplay,pphi,pphis,presnivs, & 7 7 & u,v,t,qx, & 8 8 & flxmass_w, & … … 45 45 real,intent(in) :: pphis(klon) ! surface geopotential 46 46 real,intent(in) :: presnivs(klev) ! pseudo-pressure (Pa) of mid-layers 47 integer,parameter :: longcles=2048 real,intent(in) :: clesphy0(longcles) ! Not used49 47 real,intent(in) :: u(klon,klev) ! eastward zonal wind (m/s) 50 48 real,intent(in) :: v(klon,klev) ! northward meridional wind (m/s) -
LMDZ5/trunk/libf/phylmd/lmdz1d.F90
r2191 r2221 203 203 ! Call to physiq 204 204 !--------------------------------------------------------------------- 205 integer, parameter :: longcles=20206 205 logical :: firstcall=.true. 207 206 logical :: lastcall=.false. 208 207 real :: phis = 0.0 209 real :: clesphy0(longcles) = 0.0210 208 real :: dpsrf 211 209 … … 365 363 !--------------------------------------------------------------------- 366 364 367 call conf_gcm( 99, .TRUE. , clesphy0)365 call conf_gcm( 99, .TRUE. ) 368 366 !----------------------------------------------------------------------- 369 367 ! Choix du calendrier … … 618 616 ! Ecriture du startphy avant le premier appel a la physique. 619 617 ! On le met juste avant pour avoir acces a tous les champs 620 ! NB: les clesphy0 seront remplies dans phyredem d'apres les flags lus dans gcm.def621 618 622 619 if (ok_writedem) then … … 859 856 & firstcall,lastcall, & 860 857 & day,time,timestep, & 861 & plev,play,phi,phis,presnivs, clesphy0,&858 & plev,play,phi,phis,presnivs, & 862 859 & u,v,temp,q,omega2, & 863 860 & du_phys,dv_phys,dt_phys,dq,dpsrf, & -
LMDZ5/trunk/libf/phylmd/physiq.F90
r2213 r2221 4 4 SUBROUTINE physiq (nlon,nlev, & 5 5 debut,lafin,jD_cur, jH_cur,pdtphys, & 6 paprs,pplay,pphi,pphis,presnivs, clesphy0,&6 paprs,pplay,pphi,pphis,presnivs, & 7 7 u,v,t,qx, & 8 8 flxmass_w, & … … 283 283 !$OMP THREADPRIVATE(ok_hf) 284 284 285 INTEGER longcles286 PARAMETER ( longcles = 20)287 REAL clesphy0( longcles)285 INTEGER,PARAMETER :: longcles=20 286 REAL,SAVE :: clesphy0(longcles) 287 !$OMP THREADPRIVATE(clesphy0) 288 288 ! 289 289 ! Variables propres a la physique -
LMDZ5/trunk/libf/phymar/physiq.F90
r2104 r2221 10 10 SUBROUTINE physiq (nlon,nlev, & 11 11 & debut,lafin,jD_cur, jH_cur,pdtphys, & 12 & paprs,pplay,pphi,pphis,ppresnivs, clesphy0,&12 & paprs,pplay,pphi,pphis,ppresnivs, & 13 13 & u,v,t,qx, & 14 14 & flxmass_w, & … … 118 118 real,intent(in) :: pphis(klon) ! surface geopotential 119 119 real,intent(in) :: ppresnivs(klev) ! pseudo-pressure (Pa) of mid-layers 120 integer,parameter :: longcles=20121 real,intent(in) :: clesphy0(longcles) ! Not used122 120 real,intent(in) :: u(klon,klev) ! eastward zonal wind (m/s) 123 121 real,intent(in) :: v(klon,klev) ! northward meridional wind (m/s)
Note: See TracChangeset
for help on using the changeset viewer.