! $Id: ini_wake.f90 5268 2024-10-23 17:02:39Z abarral $ SUBROUTINE ini_wake(wape, fip, it_wape_prescr, wape_prescr, fip_prescr, & alp_bl_prescr, ale_bl_prescr) USE print_control_mod, ONLY: lunout IMPLICIT NONE ! ************************************************************** ! * ! INI_WAKE : variables d'initialisation de la poche * ! froide, necessaires au declenchement * ! de la convection. * ! * ! * ! ************************************************************** ! Arguments ! ========= ! Input ! ----- ! wape : valeur de l'energie potentielle de la poche (WAPE) ! dans l'etat initial ! fip : valeur de la puissance incidente sur le front (FIP) ! dans l'etat initial ! Output ! ------ ! it_wape_prescr : nombre de pas de temps pendant lesquels la WAPE ! doit etre imposee. ! wape_prescr : valeur prescrite de la WAPE. ! fip_prescr : valeur prescrite de la FIP. ! ale_bl_prescr : valeur prescrite de la Ale de PBL. ! alp_bl_prescr : valeur prescrite de la Alp de PBL. ! Variables internes ! ================== ! it = nbre de pas de temps lu ! w = WAPE lue ! f = FIP lue ! alebl = Ale de PBL lue ! alpbl = Alp de PBL lue ! declarations REAL wape, fip, wape_prescr, fip_prescr INTEGER it_wape_prescr REAL ale_bl_prescr REAL alp_bl_prescr REAL it REAL w,f,alebl,alpbl ! FH A mettre si besoin dans physiq.def ! FH : voir avec JYG it = 0. w = 4. f = 0.1 alebl = 4. alpbl = 0.1 ! CR: on rajoute ale et alp de la PBL precrits OPEN (99, FILE='ini_wake_param.data', FORM='formatted', STATUS='old', & ERR=902) READ (99, *) it READ (99, *) w READ (99, *) f READ (99, *, END=901) alebl READ (99, *, END=901) alpbl 901 CLOSE (99) 902 CONTINUE WRITE (lunout, *) ' it,wape ', it, wape it_wape_prescr = it IF (w<0) THEN wape_prescr = wape fip_prescr = fip ELSE wape_prescr = w fip_prescr = f END IF WRITE (lunout, *) ' alebl, alpbl ', alebl, alpbl ale_bl_prescr = alebl alp_bl_prescr = alpbl PRINT *, 'Initialisation de la poche : WAPE, FIP imposees =', wape_prescr, & fip_prescr PRINT *, ' pendant ', it_wape_prescr, ' steps' PRINT *, 'Initialisation de la BL: ALP, ALE imposees =', alp_bl_prescr, & ale_bl_prescr RETURN END SUBROUTINE ini_wake