{5} Assigned, Active Tickets by Owner (Full Description) (9 matches)

The "Stop" instruction should not be used... Use "CALL abort_physic(...)" instead.

List of files to fix: (at rev 3115)

# grep -i -w "STOP" *
acama_gwd_rando_m.F90:       STOP 1
acama_gwd_rando_m.F90:       STOP 1
add_phys_tend_mod.F90:! Impression et stop en cas de probleme important
alpale_th.F90:    print *,'In order to run with iflag_coupl=2, you have to comment out the following stop'
alpale_th.F90:             STOP
clouds_gno.F90:          ! stop
convect2.F90:        ! stop
convect2.F90:    ! stop
grep: cosp: Is a directory
cv30_routines.F90:           stop
cv3_routines.F90:!jyg<  (loops stop at nl)
cv3_routines.F90:!jyg<  (loops stop at nl)
cv3_routines.F90:!jyg<  (loops stop at nl)
cv3_routines.F90:!           stop
cva_driver.F90:!          '. Might as well stop here.'
cva_driver.F90:!        STOP
cv_routines.F90:        ! stop
grep: Dust: Is a directory
grep: dyn1d: Is a directory
fisrtilp.F90:!        " but iflag_fisrtilp_qsat=",iflag_fisrtilp_qsat, ". Might as well stop here."
fisrtilp.F90:!         stop
flott_gwd_rando_m.F90:       STOP 1
flott_gwd_rando_m.F90:       STOP 1
inifis_mod.F90:        ! stop here if the relative difference is more than 1%
inifis_mod.F90:        ! stop here if the relative difference is more than 1%
inifis_mod.F90:        ! stop here if the relative difference is more than 1%
inifis_mod.F90:        ! stop here if the relative difference is more than 1%
inifis_mod.F90:        ! stop here if the relative difference is more than 1%
isccp_cloud_types.F90:  ! STOP
isccp_cloud_types.F90:  ! STOP
isccp_cloud_types.F90:  ! STOP
isccp_cloud_types.F90:  ! STOP
isccp_cloud_types.F90:  ! STOP
isccp_cloud_types.F90:  ! STOP
mod_synchro_omp.F90:      STOP 'synchro_omp'
phyaqua_mod.F90:      STOP 'probleme de dimensions dans iniaqua'
phyaqua_mod.F90:    IF (klon/=nlon) STOP 'probleme de dimensions dans iniaqua'
physiq_mod.F90:               '(H2Ov, H2Ol, H2Oi) but nqo=', nqo, '. Might as well stop here.'
physiq_mod.F90:          STOP
phytrac_mod.F90:          STOP
ran0_vec.F90:  ! stop
readaerosol_interp.F90:                 WRITE(lunout,*) 'stop for aerosol : ',name_aero(id_aero)
readaerosol_interp.F90:                 WRITE(lunout,*) 'stop for aerosol : ',name_aero(id_aero)
readaerosolstrato.F90:        STOP
readaerosolstrato.F90:       STOP
readaerosolstrato.F90:       STOP
readaerosolstrato.F90:       STOP
readaerosolstrato.F90:       STOP
readchlorophyll.F90:       STOP
readchlorophyll.F90:       STOP
readchlorophyll.F90:       STOP
grep: rrtm: Is a directory
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        write(*,*) 'som_dez = 0 STOP'
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
simu_airs.F90:        STOP
grep: sisvat: Is a directory
slab_heat_transp_mod.F90:      stop
slab_heat_transp_mod.F90:    STOP
grep: StratAer: Is a directory
stratosphere_mask.F90:  stop
thermcell_dry.F90:!               stop'On tombe sur le cas particulier de thermcell_dry'
thermcell.F90:    PRINT *, 'STOP dans convadj'
thermcell.F90:  ! if(wa_moy(1,4).gt.1.e-10) stop
thermcell_main.F90:         PRINT*,'STOP dans convadj'
thermcell_old.F90:    PRINT *, 'STOP dans convadj'
thermcell_old.F90:  ! if(wa_moy(1,4).gt.1.e-10) stop
thermcell_old.F90:    PRINT *, 'STOP dans convadj'
thermcell_old.F90:    PRINT *, 'STOP dans convadj'
thermcell_old.F90:    PRINT *, 'STOP dans convadj'
thermcell_old.F90:  ! if(wa_moy(1,4).gt.1.e-10) stop
thermcell_old.F90:    PRINT *, 'STOP dans convadj'
thermcell_plume.F90:!               stop'On tombe sur le cas particulier de thermcell_dry'
thermcell_plume.F90:!               stop'On tombe sur le cas particulier de thermcell_dry'
thermcellV0_main.F90:         PRINT*,'STOP dans convadj'
thermcellV0_main.F90:!     if(icount.eq.501) stop'au pas 301 dans thermcell_main'
vdif_kcay.F90:    ! stop
yamada4.F90:    STOP 'probleme de coherence dans appel a MY'
yamada4.F90:    STOP 'Cas nom prevu dans yamada4'
yamada_c.F90:!#     STOP'Ne pas utiliser iflag_pbl=29'
yamada_c.F90:!#     STOP'Ne pas utiliser iflag_pbl=25'

Simp^{lification et vérification}

Il faut mieux gérer le rename des executables issus de la compilation par makelmdz_fcm dans install_lmdz.sh

Sauf erreur, les schémas d'advection qui utilisent ppm3d ne peuvent pas marcher avec les fichiers sources tels qu'ils sont actuellement (révision 2449). En effet, ligne 463 de ppm3d.F :

     do 55 i=1,IMR*JMR
	55    CRY(i,2) = DTDY*V(i,1,k)

On passe forcément par cette ligne. Le tableau cry est déclaré ainsi :

real CRY(IMR,JNP)

et :

jmr = jnp - 1

Il y a donc un dépassement des limites du tableau. (Pour que ça marche, il faudrait qu'il y ait 2 * imr * jmr éléments dans le tableau.)

Le fichier est comme ça depuis au moins 12 ans (révision 524).

Olivier writes:

j'aimerais passer la décomposition direct/diffus pour le SW à la surface à destination d'orchidee pour des histoires de NPP. Pour cela, je propose de déclarer une nouvelle variable dans phys_state_var_mod.F90, la faire transiter de radlwsw.F90 vers pbl_surface_mod.F90 via physiq_mod.F90. Puis la rendre visible à ORCHIDEE via le tableau fields_out dans surf_land_orchidee_mod.F90.

Frédérique replied: C'est la methode que j'ai adoptée apres discussion avec josefine pour le freinage par les bosquets.

Puis Le 09/07/2020 à 10:06, Olivier Boucher a écrit :

Ce faisant, je butte sur la ligne suivante dans pbl_surface.F90

solsw(i,nsrf) = solsw_m(i) * (1.-alb(i,nsrf)) / (1.-alb_m(i))

On reconstruit le SW net pour la sous-surface nsrf en faisant l'hypothèse que le SW down est le même sur toutes les sous-surfaces. Mais alb correspond à l'albédo pour le rayonnement direct, alors qu'il me semble qu'il faudrait prendre l'albédo total. Je soupçonne que c'est quelque chose que Sunghye n'a pas mis à jour quand elle a introduit les albédos diffus (?).

Oui, très probablement. Tant que tu y es, ça vaut le coupe de vérifier que l'albédo est bien l'albédo moyen sur le spectre solaire, et non l'albédo d'une bande SW particulière. Ca n'a pas toujours été très carré...

oui c'est le cas, mais avec un spectre typique pas le spectre de chaque point de grille.

Je me dis que tout ça ne conserve pas l'énergie puisque le rayonnement voit un albédo et les surfaces en voit un autre. Il faudrait avoir toute l'info de ce qui se passe dans le rayonnement pour mieux faire cette répartition.

Normalement ça conserve par construction. En effet alb_m(i) = sum_nsrf w(nsrf) alb(i,nsrf) avec w la fraction de la sous surface et sum_nsrf w(nsrf)=1

donc, par construction sum_nsrf w(nsrf) solsw(i,nsrf) = solsw_m(i)

donc on conserve l'énergie

JL