Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/cva_driver.F90
- Timestamp:
- Jul 22, 2024, 9:29:09 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cva_driver.F90
r5081 r5099 367 367 ! 14 No moist convection; atmosphere is very 368 368 ! stable (=> no computation) 369 !370 369 371 370 ! ft: Array of temperature tendency (K/s) of dimension ND, defined at same … … 857 856 ! If the fraction of convective points is larger than comp_threshold, then compression 858 857 ! is assumed useless. 859 ! 858 860 859 compress = ncum < len*comp_threshold 861 ! 860 862 861 IF (.not. compress) THEN 863 862 DO i = 1,len … … 865 864 ENDDO 866 865 ENDIF 867 ! 866 868 867 !>jyg 869 868 if (prt_level >= 9) & … … 1096 1095 !! wdtrainA, wdtrainM) ! RomP 1097 1096 wdtrainA, wdtrainS, wdtrainM) !!jygprl 1098 ! 1097 1099 1098 IF (prt_level >= 10) THEN 1100 1099 Print *, 'cva_driver after cv3_unsat:mp , water, ice, evap, fondue ' … … 1110 1109 ENDDO 1111 1110 ENDIF 1112 ! 1111 1113 1112 END IF !(iflag_con==3) 1114 1113 … … 1152 1151 !! ftd, fqd, qnk, qtc, sigt, tau_cld_cv, coefw_cld_cv) 1153 1152 ftd, fqd, qta, qtc, sigt, detrain, tau_cld_cv, coefw_cld_cv) !!jygprl 1154 ! 1153 1155 1154 ! Test conseravtion de l'eau 1156 ! 1155 1157 1156 IF (debut) THEN 1158 1157 PRINT *, ' cv3_yield -> fqd(1) = ', fqd(igout, 1) 1159 1158 END IF !(debut) THEN 1160 ! 1159 1161 1160 IF (prt_level >= 10) THEN 1162 1161 Print *, 'cva_driver after cv3_yield:ft(1) , ftd(1) ', & … … 1165 1164 fq(igout,1), fqd(igout,1) 1166 1165 ENDIF 1167 ! 1166 1168 1167 END IF 1169 1168 … … 1235 1234 wdtrainA1, wdtrainS1, wdtrainM1, & ! RomP 1236 1235 qtc1, sigt1, detrain1, epmax_diag1) ! epmax_cape 1237 ! 1236 1238 1237 IF (prt_level >= 10) THEN 1239 1238 Print *, 'cva_driver after cv3_uncompress:ft1(1) , ftd1(1) ', & … … 1242 1241 fq1(igout,1), fqd1(igout,1) 1243 1242 ENDIF 1244 ! 1243 1245 1244 END IF 1246 1245 … … 1260 1259 1261 1260 END IF ! ncum>0 1262 ! 1263 ! 1261 1262 1264 1263 DO i = 1,len 1265 1264 IF (iflag1(i) == 14) THEN … … 1269 1268 ENDDO 1270 1269 1271 !1272 1270 ! In order take into account the possibility of changing the compression, 1273 1271 ! reset m, sig and w0 to zero for non-convective points.
Note: See TracChangeset
for help on using the changeset viewer.