Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/obsolete
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (4 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/obsolete
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_conserv_m.F90
r2788 r5086 97 97 vt(it) = vt(it) + (xs(is + 1) - xs(is)) * vs(is) 98 98 is = is + 1 99 end do99 END DO 100 100 ! 1 <= is <= ns 101 101 vt(it) = (vt(it) + mean_lin(xs(is), xt(it + 1)) * (xt(it + 1) & … … 105 105 if (xs(is + 1) == xt(it + 1)) is = is + 1 106 106 ! 1 <= is <= ns .or. it == nt 107 end do107 END DO 108 108 109 109 contains … … 170 170 vt(it, :) = vt(it, :) + (xs(is + 1) - xs(is)) * vs(is, :) 171 171 is = is + 1 172 end do172 END DO 173 173 ! 1 <= is <= ns 174 174 vt(it, :) = (vt(it, :) + mean_lin(xs(is), xt(it + 1)) * (xt(it + 1) & … … 178 178 if (xs(is + 1) == xt(it + 1)) is = is + 1 179 179 ! 1 <= is <= ns .or. it == nt 180 end do180 END DO 181 181 182 182 contains … … 246 246 vt(it, :, :) = vt(it, :, :) + (xs(is + 1) - xs(is)) * vs(is, :, :) 247 247 is = is + 1 248 end do248 END DO 249 249 ! 1 <= is <= ns 250 250 vt(it, :, :) = (vt(it, :, :) + mean_lin(xs(is), xt(it + 1)) & … … 254 254 if (xs(is + 1) == xt(it + 1)) is = is + 1 255 255 ! 1 <= is <= ns .or. it == nt 256 end do256 END DO 257 257 258 258 contains … … 324 324 * vs(is, :, :, :) 325 325 is = is + 1 326 end do326 END DO 327 327 ! 1 <= is <= ns 328 328 vt(it, :, :, :) = (vt(it, :, :, :) + mean_lin(xs(is), xt(it + 1)) & … … 332 332 if (xs(is + 1) == xt(it + 1)) is = is + 1 333 333 ! 1 <= is <= ns .or. it == nt 334 end do334 END DO 335 335 336 336 contains -
LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_lint_m.F90
r2788 r5086 53 53 vt(it) = ((xs(is_b+1) - xt(it)) * vs(is_b) & 54 54 + (xt(it) - xs(is_b)) * vs(is_b+1)) / (xs(is_b+1) - xs(is_b)) 55 end do55 END DO 56 56 57 57 end function regr11_lint … … 92 92 vt(it, :) = ((xs(is_b+1) - xt(it)) * vs(is_b, :) & 93 93 + (xt(it) - xs(is_b)) * vs(is_b+1, :)) / (xs(is_b+1) - xs(is_b)) 94 end do94 END DO 95 95 96 96 end function regr12_lint -
LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_step_av_m.F90
r2440 r5086 73 73 is = is + 1 74 74 left_edge = xs(is) 75 end do75 END DO 76 76 ! 1 <= is <= ns 77 77 vt(it) = (vt(it) + (xt(it + 1) - left_edge) * vs(is)) & … … 79 79 if (xs(is + 1) == xt(it + 1)) is = is + 1 80 80 ! 1 <= is <= ns .or. it == nt 81 end do81 END DO 82 82 83 83 end function regr11_step_av … … 133 133 is = is + 1 134 134 left_edge = xs(is) 135 end do135 END DO 136 136 ! 1 <= is <= ns 137 137 vt(it, :) = (vt(it, :) + (xt(it + 1) - left_edge) * vs(is, :)) & … … 139 139 if (xs(is + 1) == xt(it + 1)) is = is + 1 140 140 ! 1 <= is <= ns .or. it == nt 141 end do141 END DO 142 142 143 143 end function regr12_step_av … … 194 194 is = is + 1 195 195 left_edge = xs(is) 196 end do196 END DO 197 197 ! 1 <= is <= ns 198 198 vt(it, :, :) = (vt(it, :, :) & … … 200 200 if (xs(is + 1) == xt(it + 1)) is = is + 1 201 201 ! 1 <= is <= ns .or. it == nt 202 end do202 END DO 203 203 204 204 end function regr13_step_av … … 256 256 is = is + 1 257 257 left_edge = xs(is) 258 end do258 END DO 259 259 ! 1 <= is <= ns 260 260 vt(it, :, :, :) = (vt(it, :, :, :) + (xt(it + 1) - left_edge) & … … 262 262 if (xs(is + 1) == xt(it + 1)) is = is + 1 263 263 ! 1 <= is <= ns .or. it == nt 264 end do264 END DO 265 265 266 266 end function regr14_step_av -
LMDZ6/branches/Amaury_dev/libf/obsolete/regr3_lint_m.F90
r2788 r5086 53 53 vt(:, :, it) = ((xs(is_b+1) - xt(it)) * vs(:, :, is_b) & 54 54 + (xt(it) - xs(is_b)) * vs(:, :, is_b+1)) / (xs(is_b+1) - xs(is_b)) 55 end do55 END DO 56 56 57 57 end function regr33_lint … … 94 94 + (xt(it) - xs(is_b)) * vs(:, :, is_b+1, :)) & 95 95 / (xs(is_b+1) - xs(is_b)) 96 end do96 END DO 97 97 98 98 end function regr34_lint -
LMDZ6/branches/Amaury_dev/libf/obsolete/regr_lat_time_climoz_m.F90
r2788 r5086 249 249 do while (o3_in(j, 1, l, m) == missing_value) 250 250 j = j + 1 251 end do251 END DO 252 252 if (j > 1) o3_in(:j-1, :, l, m) = & 253 253 spread(o3_in(j, :, l, m), dim=1, ncopies=j-1) … … 257 257 do while (o3_in(j, 1, l, m) == missing_value) 258 258 j = j - 1 259 end do259 END DO 260 260 if (j < n_lat) o3_in(j+1:, :, l, m) = & 261 261 spread(o3_in(j, :, l, m), dim=1, ncopies=n_lat-j) … … 270 270 do while (o3_in(j, k, l, m) /= missing_value .and. k < n_plev) 271 271 k = k + 1 272 end do272 END DO 273 273 ! Replace missing values with the valid value at the 274 274 ! lowest level above missing values: 275 275 if (o3_in(j, k, l, m) == missing_value) & 276 276 o3_in(j, k:n_plev, l, m) = o3_in(j, k-1, l, m) 277 end do278 end do277 END DO 278 END DO 279 279 else 280 280 print *, "regr_lat_time_climoz: field ", m, & 281 281 ", no missing value attribute" 282 282 end if 283 end do283 END DO 284 284 285 285 call nf95_close(ncid_in) … … 322 322 call nf95_put_var(ncid_out, varid_out(m), o3_out(nbp_lat:1:-1, :, :, m)) 323 323 ! (The order of "rlatu" is inverted in the output file) 324 end do324 END DO 325 325 326 326 call nf95_close(ncid_out) -
LMDZ6/branches/Amaury_dev/libf/obsolete/regr_pr_av_m.F90
r2788 r5086 99 99 call handle_err("regr_pr_av nf90_get_var " // trim(name(i)), ncerr, & 100 100 ncid) 101 end do101 END DO 102 102 103 103 ! Latitudes are in ascending order in the input file while … … 118 118 slopes(v2(i, :, :), press_in_edg)) 119 119 ! (invert order of indices because "paprs" is in descending order) 120 end do120 END DO 121 121 122 122 end subroutine regr_pr_av
Note: See TracChangeset
for help on using the changeset viewer.