Changeset 113


Ignore:
Timestamp:
06/13/17 16:58:07 (7 years ago)
Author:
dumas
Message:

Update water cycle : ablation_bord and calving_frange corrected : results are now very close to diff_H | In main3D, lineartemp, lect-anteis_mod and lect-greeneem_mod replace H=1 by H=0

Location:
trunk/SOURCES
Files:
7 edited

Legend:

Unmodified
Added
Removed

  • trunk/SOURCES/Ant40_files/lect-anteis_mod.f90

    r86 r113  
    182182!    endwhere 
    183183    S(:,:)=max(S(:,:),0.)   ! pour etre au niveau des mers : ATTENTION si SEALEV <0 
    184     H(:,:)=max(H(:,:),1.)   ! pour avoir au moins 1 m 
     184    H(:,:)=max(H(:,:),0.)   ! pas d epaisseurs negatives 
    185185 
    186186 

  • trunk/SOURCES/Ant40_files/output_anta40_mod-0.4.f90

    r99 r113  
    1414 
    1515USE module3D_phy 
     16use bilan_eau_mod 
    1617 
    1718 

  • trunk/SOURCES/Greeneem_files/lect-greeneem_mod.f90

    r85 r113  
    7979    call lect_input(1,'H',1,H,topo_dep,file_ncdf)          ! epaisseur 
    8080    S(:,:)=max(S(:,:),0.)   ! pour etre au niveau des mers : ATTENTION si SEALEV <0 
    81     H(:,:)=max(H(:,:),1.)   ! pour avoir au moins 1 m 
     81    H(:,:)=max(H(:,:),0.)   ! pas d epaisseurs negatives 
    8282 
    8383           

  • trunk/SOURCES/ablation_bord.f90

    r102 r113  
    2525integer :: somm 
    2626 
    27  
    28 ! Dans les zones sans glace ou le calcul icethick a donne un H negatif 
    29 ! ablbord contient pour l'instant Hneg/dt  
    30 ! ailleurs ablbord=0 
    31  
    32 ! mise a jour du masque "ice" 
    33 where (flot(:,:))        ! points flottants, sera éventuellement réévalué dans flottab 
    34 !cdc 1m  where(H(:,:).gt.max(Hmin,Hmin+BM(:,:)-Bmelt(:,:))) 
    35    where(H(:,:).gt.0.) 
    36       ice(:,:)=1 
    37    elsewhere 
    38       ice(:,:)=0 
    39 !cdc 1m      H(:,:)=max(1.,min(0.,(sealevel - Bsoc(:,:))*row/ro-0.01)) 
    40                         H(:,:)=0. 
    41    endwhere 
    42 elsewhere                ! points posés 
    43    where(H(:,:).gt.0.) 
    44       ice(:,:)=1 
    45    elsewhere 
    46       ice(:,:)=0 
    47    endwhere 
    48 endwhere 
    49  
    50  
    51  
    52  
    53  
    54 where (ice(:,:).eq.1)    !  quand glace : abl=bm-acc  
    55    ablbord(:,:)=0. 
    56 endwhere 
    57 !~ where (calv(:,:).LT.0.) ! il ne peut pas y avoir calving et ablbord sur un meme point 
    58 !~       ablbord(:,:)=0. 
    59 !~ endwhere 
     27! Dans les zones ou le calcul icethick a donne un H negatif : ablbord=H 
    6028 
    6129! pour être au bord, il suffit que la somme des ice voisins soit > 0 

  • trunk/SOURCES/calving_frange.f90

    r112 r113  
    410410       
    411411! on met en calving les points detectes iceberg : 
    412        where (iceberg(:,:))    
     412       where (iceberg(:,:)) 
     413                                                 calv(:,:)=-h(:,:) 
    413414             ice(:,:)=0 
    414415             H(:,:)=0. 

  • trunk/SOURCES/lineartemp-0.2.f90

    r29 r113  
    4747        !     T(I,J,NZ)=TPMP(I,J,NZ)-10. 
    4848 
    49         if (H(i,j).gt.1.)  T(I,J,NZ)=TPMP(I,J,NZ)-20. 
     49        if (H(i,j).gt.0.)  T(I,J,NZ)=TPMP(I,J,NZ)-20. 
    5050 
    5151 

  • trunk/SOURCES/main3D-0.4-40km.f90

    r112 r113  
    232232  do j=1,ny 
    233233     do i=1,nx 
    234         H(i,j)=max(1.,H(i,j)) 
     234        H(i,j)=max(0.,H(i,j)) 
    235235     enddo 
    236236  enddo 
Note: See TracChangeset for help on using the changeset viewer.