Changeset 136 in lmdz_wrf
- Timestamp:
- Jul 29, 2014, 11:32:34 AM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/LMDZ_WRFmeas/WRFV3/lmdz/thermcell_dq.F90
r77 r136 38 38 CHARACTER (LEN=80) :: abort_message 39 39 40 41 ! Lluis42 INTEGER :: llp43 CHARACTER(LEN=50) :: lvarname, lfname44 REAL :: largest45 46 llp = 73447 lfname = 'physiq'48 largest = 10.e549 50 40 ! Old explicite scheme 51 lfname='before old scheme'52 lvarname = 'dq'53 CALL check_var3D(lfname, lvarname, dq, ngrid, nlay, largest, .FALSE.)54 55 41 if (impl==-1) then 56 42 call thermcell_dq_o(ngrid,nlay,ptimestep,fm,entr, & … … 58 44 return 59 45 endif 60 lfname='after old scheme'61 lvarname = 'dq'62 CALL check_var3D(lfname, lvarname, dq, ngrid, nlay, largest, .FALSE.)63 46 64 47 ! Calcul du critere CFL pour l'advection dans la subsidence … … 131 114 fqa(:,1)=0. ; fqa(:,nlay)=0. 132 115 133 134 lfname='after Plume vertical flux'135 lvarname = 'qa'136 CALL check_var3D(lfname, lvarname, qa, ngrid, nlay, largest, .FALSE.)137 lvarname = 'fm'138 CALL check_var3D(lfname, lvarname, fm, ngrid, nlay, largest, .FALSE.)139 lvarname = 'fqa'140 CALL check_var3D(lfname, lvarname, fqa, ngrid, nlay, largest, .FALSE.)141 lvarname = 'q'142 CALL check_var3D(lfname, lvarname, q, ngrid, nlay, largest, .FALSE.)143 lvarname = 'masse'144 CALL check_var3D(lfname, lvarname, masse, ngrid, nlay, largest, .FALSE.)145 146 116 ! Trace species evolution 147 117 if (impl==0) then … … 156 126 enddo 157 127 endif 158 lfname='after tendencies _dq 000'159 lvarname = 'q'160 CALL check_var3D(lfname, lvarname, q, ngrid, nlay, largest, .FALSE.)161 lvarname = 'dq'162 CALL check_var3D(lfname, lvarname, dq, ngrid, nlay, largest, .FALSE.)163 lvarname = 'qold'164 CALL check_var3D(lfname, lvarname, qold, ngrid, nlay, largest, .FALSE.)165 128 166 129 ! Tendencies … … 171 134 enddo 172 135 enddo 173 lfname='after tendencies _dq'174 lvarname = 'q'175 CALL check_var3D(lfname, lvarname, q, ngrid, nlay, largest, .FALSE.)176 lvarname = 'dq'177 CALL check_var3D(lfname, lvarname, dq, ngrid, nlay, largest, .FALSE.)178 lvarname = 'qold'179 CALL check_var3D(lfname, lvarname, qold, ngrid, nlay, largest, .FALSE.)180 136 181 137 return
Note: See TracChangeset
for help on using the changeset viewer.