- Timestamp:
- Jul 20, 2024, 12:00:23 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/surf_land_bucket_mod.F90
r5082 r5087 130 130 IF (iso_eau > 0) THEN 131 131 CALL iso_verif_egalite_choix(precip_snow(i), & 132 &xtprecip_snow(iso_eau,i),'surf_land_bucket 131', &133 &errmax,errmaxrel)132 xtprecip_snow(iso_eau,i),'surf_land_bucket 131', & 133 errmax,errmaxrel) 134 134 CALL iso_verif_egalite_choix(qsol(i), & 135 &xtsol(iso_eau,i),'surf_land_bucket 134', &136 &errmax,errmaxrel)135 xtsol(iso_eau,i),'surf_land_bucket 134', & 136 errmax,errmaxrel) 137 137 ENDIF 138 138 ENDDO … … 171 171 IF (soil_model) THEN 172 172 CALL soil(dtime, is_ter, knon, snow, tsurf, qsol, & 173 &longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)173 longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux) 174 174 175 175 DO i=1, knon … … 204 204 IF (iso_eau > 0) THEN 205 205 CALL iso_verif_egalite_choix(xtsnow(iso_eau,i), & 206 &snow(i),'surf_land_bucket 522', &207 & errmax,errmaxrel)206 snow(i),'surf_land_bucket 522', & 207 errmax,errmaxrel) 208 208 ENDIF !IF (iso_eau > 0) then 209 209 ENDDO !DO i=1,knon … … 229 229 IF (iso_eau > 0) THEN 230 230 CALL iso_verif_egalite(qsol_prec(i),xtsol_prec(iso_eau,i), & 231 &'surf_land_bucket 141')231 'surf_land_bucket 141') 232 232 ENDIF 233 233 ENDDO !DO i=1,knon … … 241 241 snow, qsol, tsurf_new, evap & 242 242 #ifdef ISO 243 &,fq_fonte_diag,fqfonte_diag,snow_evap_diag,fqcalving_diag &244 &,max_eau_sol_diag,runoff_diag,run_off_lic_diag,coeff_rel_diag &245 #endif 246 &)243 ,fq_fonte_diag,fqfonte_diag,snow_evap_diag,fqcalving_diag & 244 ,max_eau_sol_diag,runoff_diag,run_off_lic_diag,coeff_rel_diag & 245 #endif 246 ) 247 247 248 248 #ifdef ISO … … 259 259 IF (iso_eau > 0) THEN 260 260 CALL iso_verif_egalite_choix(qsol_prec(i), & 261 &xtsol_prec(iso_eau,i),'surf_land_bucket 628', &262 &errmax,errmaxrel)261 xtsol_prec(iso_eau,i),'surf_land_bucket 628', & 262 errmax,errmaxrel) 263 263 CALL iso_verif_egalite_choix(precip_snow(i), & 264 &xtprecip_snow(iso_eau,i),'surf_land_bucket 227', &265 &errmax,errmaxrel)264 xtprecip_snow(iso_eau,i),'surf_land_bucket 227', & 265 errmax,errmaxrel) 266 266 ! attention, dans fonte_neige, on modifie snow sans modifier 267 267 ! xtsnow … … 270 270 ! write(*,*) 'snow(i)=',snow(i) 271 271 CALL iso_verif_egalite_choix(xtsnow(iso_eau,i), & 272 &snow_prec(i),'surf_land_bucket 245', &273 & errmax,errmaxrel)272 snow_prec(i),'surf_land_bucket 245', & 273 errmax,errmaxrel) 274 274 ENDIF 275 275 IF ((iso_O17 > 0).AND.(iso_O18 > 0)) THEN 276 276 IF (qsol_prec(i) > ridicule_qsol) THEN 277 277 CALL iso_verif_aberrant_o17(xtsol_prec(iso_O17,i)/qsol_prec(i) & 278 &,xtsol_prec(iso_O18,i)/qsol_prec(i) &279 &,'surf_land_bucket 642')278 ,xtsol_prec(iso_O18,i)/qsol_prec(i) & 279 ,'surf_land_bucket 642') 280 280 ENDIF !IF ((qsol_prec(i) > ridicule_qsol) & 281 281 ENDIF !IF ((iso_O17 > 0).AND.(iso_O18 > 0)) THEN … … 285 285 #endif 286 286 CALL calcul_iso_surf_ter_vectall(klon,knon, & 287 &evap,snow_evap_diag,snow, &288 &fq_fonte_diag,fqfonte_diag,dtime,precip_rain,xtprecip_rain, &289 &precip_snow,xtprecip_snow, snow_prec,xtsnow_prec, &290 &tsurf_new,xtspechum,pref,spechum,t_coup,u1_lay,v1_lay,p1lay, &291 &qsol,xtsol,qsol_prec,xtsol_prec, &292 &max_eau_sol_diag, &293 &xtevap,xtsnow,h1,runoff_diag,xtrunoff_diag,fqcalving_diag, &294 &knindex,is_ter,run_off_lic_diag,coeff_rel_diag,Rland_ice &295 &)287 evap,snow_evap_diag,snow, & 288 fq_fonte_diag,fqfonte_diag,dtime,precip_rain,xtprecip_rain, & 289 precip_snow,xtprecip_snow, snow_prec,xtsnow_prec, & 290 tsurf_new,xtspechum,pref,spechum,t_coup,u1_lay,v1_lay,p1lay, & 291 qsol,xtsol,qsol_prec,xtsol_prec, & 292 max_eau_sol_diag, & 293 xtevap,xtsnow,h1,runoff_diag,xtrunoff_diag,fqcalving_diag, & 294 knindex,is_ter,run_off_lic_diag,coeff_rel_diag,Rland_ice & 295 ) 296 296 !#ifdef ISOVERIF 297 297 ! write(*,*) 'surf_land_bucket 303'
Note: See TracChangeset
for help on using the changeset viewer.