- Timestamp:
- Aug 2, 2024, 2:12:03 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/prec_scops.F90
r5099 r5158 66 66 if (cv_col .eq. 0) cv_col=1 67 67 68 doilev=1,nlev69 doibox=1,ncol70 do j=1,npoints68 DO ilev=1,nlev 69 DO ibox=1,ncol 70 DO j=1,npoints 71 71 prec_frac(j,ibox,ilev) = 0 72 72 enddo … … 74 74 enddo 75 75 76 doj=1,npoints77 doibox=1,ncol76 DO j=1,npoints 77 DO ibox=1,ncol 78 78 frac_out_ls(j,ibox)=0 79 79 frac_out_cv(j,ibox)=0 80 80 flag_ls=0 81 81 flag_cv=0 82 doilev=1,nlev82 DO ilev=1,nlev 83 83 if (frac_out(j,ibox,ilev) .eq. 1) then 84 84 flag_ls=1 … … 98 98 99 99 ! initialize the top layer 100 doj=1,npoints100 DO j=1,npoints 101 101 flag_ls=0 102 102 flag_cv=0 103 103 104 104 if (ls_p_rate(j,1) .gt. 0.) then 105 doibox=1,ncol ! possibility ONE105 DO ibox=1,ncol ! possibility ONE 106 106 if (frac_out(j,ibox,1) .eq. 1) then 107 107 prec_frac(j,ibox,1) = 1 … … 110 110 enddo ! loop over ncol 111 111 if (flag_ls .eq. 0) then ! possibility THREE 112 doibox=1,ncol112 DO ibox=1,ncol 113 113 if (frac_out(j,ibox,2) .eq. 1) then 114 114 prec_frac(j,ibox,1) = 1 … … 118 118 endif 119 119 if (flag_ls .eq. 0) then ! possibility Four 120 doibox=1,ncol120 DO ibox=1,ncol 121 121 if (frac_out_ls(j,ibox) .eq. 1) then 122 122 prec_frac(j,ibox,1) = 1 … … 126 126 endif 127 127 if (flag_ls .eq. 0) then ! possibility Five 128 doibox=1,ncol128 DO ibox=1,ncol 129 129 ! prec_frac(j,1:ncol,1) = 1 130 130 prec_frac(j,ibox,1) = 1 … … 135 135 136 136 if (cv_p_rate(j,1) .gt. 0.) then 137 doibox=1,ncol ! possibility ONE137 DO ibox=1,ncol ! possibility ONE 138 138 if (frac_out(j,ibox,1) .eq. 2) then 139 139 if (prec_frac(j,ibox,1) .eq. 0) then … … 146 146 enddo ! loop over ncol 147 147 if (flag_cv .eq. 0) then ! possibility THREE 148 doibox=1,ncol148 DO ibox=1,ncol 149 149 if (frac_out(j,ibox,2) .eq. 2) then 150 150 if (prec_frac(j,ibox,1) .eq. 0) then … … 158 158 endif 159 159 if (flag_cv .eq. 0) then ! possibility Four 160 doibox=1,ncol160 DO ibox=1,ncol 161 161 if (frac_out_cv(j,ibox) .eq. 1) then 162 162 if (prec_frac(j,ibox,1) .eq. 0) then … … 170 170 endif 171 171 if (flag_cv .eq. 0) then ! possibility Five 172 doibox=1,cv_col172 DO ibox=1,cv_col 173 173 if (prec_frac(j,ibox,1) .eq. 0) then 174 174 prec_frac(j,ibox,1) = 2 … … 187 187 188 188 ! working on the levels from top to surface 189 doilev=2,nlev190 doj=1,npoints189 DO ilev=2,nlev 190 DO j=1,npoints 191 191 flag_ls=0 192 192 flag_cv=0 193 193 194 194 if (ls_p_rate(j,ilev) .gt. 0.) then 195 doibox=1,ncol ! possibility ONE&TWO195 DO ibox=1,ncol ! possibility ONE&TWO 196 196 if ((frac_out(j,ibox,ilev) .eq. 1) .or. ((prec_frac(j,ibox,ilev-1) .eq. 1) & 197 197 .or. (prec_frac(j,ibox,ilev-1) .eq. 3))) then … … 201 201 enddo ! loop over ncol 202 202 if ((flag_ls .eq. 0) .and. (ilev .lt. nlev)) then ! possibility THREE 203 doibox=1,ncol203 DO ibox=1,ncol 204 204 if (frac_out(j,ibox,ilev+1) .eq. 1) then 205 205 prec_frac(j,ibox,ilev) = 1 … … 209 209 endif 210 210 if (flag_ls .eq. 0) then ! possibility Four 211 doibox=1,ncol211 DO ibox=1,ncol 212 212 if (frac_out_ls(j,ibox) .eq. 1) then 213 213 prec_frac(j,ibox,ilev) = 1 … … 217 217 endif 218 218 if (flag_ls .eq. 0) then ! possibility Five 219 doibox=1,ncol219 DO ibox=1,ncol 220 220 ! prec_frac(j,1:ncol,ilev) = 1 221 221 prec_frac(j,ibox,ilev) = 1 … … 225 225 226 226 if (cv_p_rate(j,ilev) .gt. 0.) then 227 doibox=1,ncol ! possibility ONE&TWO227 DO ibox=1,ncol ! possibility ONE&TWO 228 228 if ((frac_out(j,ibox,ilev) .eq. 2) .or. ((prec_frac(j,ibox,ilev-1) .eq. 2) & 229 229 .or. (prec_frac(j,ibox,ilev-1) .eq. 3))) then … … 237 237 enddo ! loop over ncol 238 238 if ((flag_cv .eq. 0) .and. (ilev .lt. nlev)) then ! possibility THREE 239 doibox=1,ncol239 DO ibox=1,ncol 240 240 if (frac_out(j,ibox,ilev+1) .eq. 2) then 241 241 if (prec_frac(j,ibox,ilev) .eq. 0) then … … 249 249 endif 250 250 if (flag_cv .eq. 0) then ! possibility Four 251 doibox=1,ncol251 DO ibox=1,ncol 252 252 if (frac_out_cv(j,ibox) .eq. 1) then 253 253 if (prec_frac(j,ibox,ilev) .eq. 0) then … … 261 261 endif 262 262 if (flag_cv .eq. 0) then ! possibility Five 263 doibox=1,cv_col263 DO ibox=1,cv_col 264 264 if (prec_frac(j,ibox,ilev) .eq. 0) then 265 265 prec_frac(j,ibox,ilev) = 2
Note: See TracChangeset
for help on using the changeset viewer.