Changeset 313 for LMDZ.3.3/branches


Ignore:
Timestamp:
Jan 25, 2002, 11:59:34 AM (22 years ago)
Author:
lmdzadmin
Message:

Prise en compte de la fonte des glaciers antarctiques pour le couplé
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90

    r295 r313  
    4444!!$PB
    4545  REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: tmp_rriv, tmp_rcoa
     46!! pour simuler la fonte des glaciers antarctiques
     47  REAL, ALLOCATABLE, DIMENSION(:,:), SAVE :: coeff_iceberg
     48  real, save                              :: surf_maille
     49  real, save                              :: cte_flux_iceberg = 6.3e7
     50  integer, save                           :: num_antarctic = 1
    4651!!$
    4752  CONTAINS
     
    11231128!   itime        numero du pas de temps
    11241129!   iim, jjm     nbres de pts de grille
    1125 !   dtime        pas de tempsde la physique
     1130!   dtime        pas de temps de la physique
    11261131!   klon         nombre total de points de grille
    11271132!   nisurf       index de la surface a traiter (1 = sol continental)
     
    13431348    CALL histsync(nidcs)
    13441349
     1350! pour simuler la fonte des glaciers antarctiques
     1351!
     1352    surf_maille = (4. * rpi * ra**2) / (iim * (jjm +1))
     1353    ALLOCATE(coeff_iceberg(iim,jjm+1), stat=error)
     1354    if (error /= 0) then
     1355      abort_message='Pb allocation variable coeff_iceberg'
     1356      call abort_gcm(modname,abort_message,1)
     1357    endif
     1358    open (12,file='flux_iceberg',form='formatted',status='old')
     1359    read (12,*) coeff_iceberg
     1360    close (12)
     1361    num_antarctic = max(1, count(coeff_iceberg > 0))
     1362   
    13451363    first_appel = .false.
    13461364  endif ! fin if (first_appel)
     
    15441562      endwhere
    15451563!
     1564! pour simuler la fonte des glaciers antarctiques
     1565!
     1566        wri_rain = wri_rain      &
     1567      &     + coeff_iceberg * cte_flux_iceberg / (num_antarctic * surf_maille)
     1568
     1569!
    15461570! on passe les coordonnées de la grille
    15471571!
Note: See TracChangeset for help on using the changeset viewer.