wiki:DevelopmentActivities/MergeHydro/flag_couple_note

Version 1 (modified by nvuilsce, 10 years ago) (diff)

--

Fiche : flag couplé

Description

Plusieurs formulations (humrel, drysoil_frac, rsol, albedo) avaient été modifiées pour des applications en mode couplé afin de compenser des biais du modèle. Ces formulations proposées par Nathalie et Frédéric ne sont pas présentes dans la version LMD, version qui conserve donc les formulations originelles.

Choix

On décide de rajouter un flag spécifique pour le couplé (flag non actif par défaut) permettant d'activer les formulations proposées par Nathalie et Frédéric. Les modules impactées sont condveg.f90 pour alb_bare et hydrolc.f90 pour rsol, drysoilfrac, throughfall et humrel

Mise en oeuvre

La version à modifier est la version d'origine (1.9 ou 1.9.4)

  • dans condveg.f90
    • Traitement de alb_bare_model
      • On remplace
              alb_bare_model=.FALSE.
               CALL getin_p('ALB_BARE_MODEL', alb_bare_model)
        
        par
              ! Declaration d'une variable du module en tete de module
              LOGICAL, SAVE                         :: couple = .FALSE.
        
        +
        
              alb_bare_model = .TRUE.
              CALL getin_p('COUPLE', couple)	
              IF(couple) THEN
                     alb_bare_model = .FALSE.
               ENDIF
               CALL getin_p('ALB_BARE_MODEL', alb_bare_model)
        
  • dans hydrolc.f90
    • Traitement de rsol
  • On remplace dans hydrolc_var_init
        !
        ! Compute the resistance to bare soil evaporation
        !
        rsol(:) = -un
        DO ji = 1, kjpindex
           IF (veget(ji,1) .GE. min_sechiba) THEN
              !
              ! Correction Nathalie - le 28 mars 2006 - sur conseils Fred Hourdin
              ! on modifie le rsol pour que la resistance croisse subitement si on s'approche
              ! du fond. En gros, rsol=hdry*rsol_cste pour hdry < 1m70
              !Ancienne formulation
              rsol(ji) = dss(ji,1) * rsol_cste
              !Nouvelle formulation Nath          
              !rsol(ji) =  ( drysoil_frac(ji) + 1./(10.*(dpu_cste - drysoil_frac(ji))+1.e-10)**2 ) * rsol_cste
           ENDIF
        ENDDO
    

par

    ! Declaration d'une variable du module en tete de module
    LOGICAL, SAVE                                     ::  couple = .FALSE.

+

    ! Nouvelle Declaration locale
    REAL(r_std), DIMENSION (kjpindex) :: frsoil 

+  

    !
    ! Compute the resistance to bare soil evaporation
    !
    CALL getin_p('COUPLE', couple)	
    IF(couple) THEN
       frsoil(:)=( drysoil_frac(:) + 1./(10.*(dpu_cste – drysoil_frac(:))+1.e-10)**2 ) 
    ELSE
       frsoil(:)=dss(:,1) 
    ENDIF
    rsol(:) = -un
    DO ji = 1, kjpindex
       IF (veget(ji,1) .GE. min_sechiba) THEN
	 rsol(ji) = dss(ji,1) * rsol_cste
       ENDIF   
    ENDDO
  • On remplace dans hydrolc_soil
        ! The fraction of visibly dry soil (dry when dss(:,1) = 0.1 m)
        ! Correction Nathalie - le 28 Mars 2006 - re-ecriture drysoil_frac/hdry - Fred Hourdin
        !drysoil_frac(:) = MIN(MAX(dss(:,1),0.)*10._r_std, un)
        drysoil_frac(:) = a_subgrd(:,1)*dss(:,1) + (1.-a_subgrd(:,1))*dsp(:,1)
        !
        ! Compute the resistance to bare soil evaporation.
        !
        rsol(:) = -un
        DO ji = 1, kjpindex
           IF (veget(ji,1) .GE. min_sechiba) THEN
              !
              ! Correction Nathalie - le 28 mars 2006 - sur conseils Fred Hourdin
              ! on modifie le rsol pour que la resistance croisse subitement si on s'approche
              ! du fond. En gros, rsol=hdry*rsol_cste pour hdry < 1m70
              !rsol(ji) = dss(ji,1) * rsol_cste
              rsol(ji) =  ( drysoil_frac(ji) + 1./(10.*(dpu_cste - drysoil_frac(ji))+1.e-10)**2 ) * rsol_cste
           ENDIF
        ENDDO
    

par

    ! Nouvelle Declaration locale
    REAL(r_std), DIMENSION (kjpindex) :: frsoil 

    ! The fraction of visibly dry soil (dry when dss(:,1) = 0.1 m)
    IF (couple) THEN
       ! Correction Nathalie - le 28 Mars 2006 - re-ecriture drysoil_frac/hdry - Fred Hourdin	 
       ! on modifie le rsol pour que la resistance croisse subitement si on s'approche
       ! du fond. En gros, rsol=hdry*rsol_cste pour hdry < 1m70
       drysoil_frac(:) = a_subgrd(:,1)*dss(:,1) + (1.-a_subgrd(:,1))*dsp(:,1)
       frsoil(:)=( drysoil_frac(:) + 1./(10.*(dpu_cste – drysoil_frac(:))+1.e-10)**2 ) 
    ELSE
       drysoil_frac(:) = MIN(MAX(dss(:,1),0.)*10._r_std, un)
       frsoil(:)=dss(:,1)
    ENDIF
    !
    ! Compute the resistance to bare soil evaporation.
    !
    rsol(:) = -un
    DO ji = 1, kjpindex
       IF (veget(ji,1) .GE. min_sechiba) THEN
          rsol(ji) = frsoil(ji) * rsol_cste
       ENDIF
    ENDDO
  • Traitement de throughfall (pour calcul de precisol et qsintveg)
  • On remplace dans hydrolc_canop
        IF ( firstcall ) THEN
           !Config  Key  = PERCENT_THROUGHFALL_PFT
           !Config  Desc = Percent by PFT of precip that is not intercepted by the canopy
           !Config  Def  = 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30.
           !Config  Help = During one rainfall event, PERCENT_THROUGHFALL_PFT% of the incident rainfall
           !Config         will get directly to the ground without being intercepted, for each PFT.
           
           throughfall_by_pft = (/ 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30. /)
           CALL getin_p('PERCENT_THROUGHFALL_PFT',throughfall_by_pft)
           throughfall_by_pft = throughfall_by_pft / 100.
    
           firstcall=.FALSE.
        ENDIF
    
    par
        IF ( firstcall ) THEN
           !Config  Key  = PERCENT_THROUGHFALL_PFT
           !Config  Desc = Percent by PFT of precip that is not intercepted by the canopy
           !Config  Def  = 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30.
           !Config  Help = During one rainfall event, PERCENT_THROUGHFALL_PFT% of the incident rainfall
           !Config         will get directly to the ground without being intercepted, for each PFT.
           
           throughfall_by_pft = (/ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0. /)
    
           IF (couple) THEN
    	throughfall_by_pft = (/ 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30., 30. /)
           ENDIF
           CALL getin_p('PERCENT_THROUGHFALL_PFT',throughfall_by_pft)
     
          throughfall_by_pft = throughfall_by_pft / 100.
    
           firstcall=.FALSE.
        ENDIF
    
  • Traitement de humrel
  • On remplace dans hydrolc_soil
                ELSE
                   ! Corrections Nathalie - le 28 Mars 2006 - sur conseils Fred Hourdin
                   !zhumrel_lo(ji) = EXP( - humcste(jv) * dpu_cste * (dsp(ji,jv)/dpu_cste) )
                   !zhumrel_up(ji) = EXP( - humcste(jv) * dpu_cste * (dss(ji,jv)/dsg(ji,jv)) )
                   !humrel(ji,jv) = MAX(zhumrel_lo(ji),zhumrel_up(ji))
                   !
                   ! As we need a slower variable for vegetation growth the stress is computed
                   ! differently than in humrel.
                   !
                   zhumrel_lo(ji) = EXP( - humcste(jv) * dsp(ji,jv))
                   zhumrel_up(ji) = EXP( - humcste(jv) * dss(ji,jv))
                   ! Ajouts Nathalie - Fred - le 28 Mars 2006
                   a_subgrd(ji,jv)=MIN(MAX(dsg(ji,jv)-dss(ji,jv),0.)/dsg_min,1.)
                   humrel(ji,jv)=a_subgrd(ji,jv)*zhumrel_up(ji)+(1.-a_subgrd(ji,jv))*zhumrel_lo(ji)
                   !
                   vegstress(ji,jv) = zhumrel_lo(ji) + zhumrel_up(ji) - EXP( - humcste(jv) * dsg(ji,jv) ) 
                   !
                ENDIF
    

par

           ELSE
	   IF(couple)
	       ! Ajouts Nathalie - Fred - le 28 Mars 2006
                     zhumrel_lo(ji) = EXP( - humcste(jv) * dsp(ji,jv))
                     zhumrel_up(ji) = EXP( - humcste(jv) * dss(ji,jv))
                    a_subgrd(ji,jv)=MIN(MAX(dsg(ji,jv)-dss(ji,jv),0.)/dsg_min,1.)
                    humrel(ji,jv)=a_subgrd(ji,jv)*zhumrel_up(ji)+(1.-a_subgrd(ji,jv))*zhumrel_lo(ji)
	   ELSE
                    zhumrel_lo(ji) = EXP( - humcste(jv) * dpu_cste * (dsp(ji,jv)/dpu_cste) )
                    zhumrel_up(ji) = EXP( - humcste(jv) * dpu_cste * (dss(ji,jv)/dsg(ji,jv)) )
                    humrel(ji,jv) = MAX(zhumrel_lo(ji),zhumrel_up(ji))
	   ENDIF
             ENDIF

Tests