Changeset 91 for branches


Ignore:
Timestamp:
10/28/16 09:11:13 (8 years ago)
Author:
aquiquet
Message:

GRISLI coupled using basal melting rates coming from CLIO. Might need further testing for LGM conditions.

Location:
branches/iLoveclim/SOURCES
Files:
1 added
6 edited

Legend:

Unmodified
Added
Removed

  • branches/iLoveclim/SOURCES/3D-physique-gen_mod.f90

    r90 r91  
    253253  real (kind = 8),dimension(nx,ny) :: calvingGRIS 
    254254  real (kind = 8),dimension(nx,ny),target :: calvin_GRIS ! calving flux[m3/s] given to CLIO 
     255  real (kind = 8),dimension(nx,ny,20),target :: bmshelfCLIO  ! basal melting given by CLIO   
    255256  real, dimension(nx,ny) :: coef_drag   !< coefficient de la loi de friction non lineaire : depend de la valeur de alpha_drag 
    256257                                        !< si alpha_drag = 1, coef_drag = drag_centre   

  • branches/iLoveclim/SOURCES/Hemin40_files/module_choix-hemin40-0.4.f90

    r88 r91  
    117117!use bmelt_ant_regions ! pour l'Antarctique avec régions 
    118118!use bmelt_nor_regions ! pour le nord avec régions 
    119 use bmelt_nor_depth   ! pour le nord avec profondeur d'eau 
    120  
     119!use bmelt_nor_depth   ! pour le nord avec profondeur d'eau 
     120use bmelt_clio_coupl_mod ! coupled version: CLIO bm rates (needs SHELFMELT==1 iLOVECLIM flag) 
    121121 
    122122! fake -------fausses routines pour compatibilités 

  • branches/iLoveclim/SOURCES/Makefile.grisli-gfortran.inc

    r78 r91  
    8282        conserv-mass-adv-diff_sept2009_mod.o  \ 
    8383        eaubasale-0.5_mod.o spinup_mod.o no_spinup_mod.o \ 
    84         tridiag_mod-0.3.o out_horiz_mod.o bmelt-seuil-profondeur_mod.o  
     84        tridiag_mod-0.3.o out_horiz_mod.o \ 
     85        bmelt_clio_coupl_mod.o bmelt-seuil-profondeur_mod.o  
    8586 
    8687mod_communs_RGL =   util_recovery.o \ 

  • branches/iLoveclim/SOURCES/Makefile.grisli.inc

    r78 r91  
    8181        conserv-mass-adv-diff_sept2009_mod.o  \ 
    8282        eaubasale-0.5_mod.o spinup_mod.o no_spinup_mod.o \ 
    83         tridiag_mod-0.3.o out_horiz_mod.o bmelt-seuil-profondeur_mod.o  
     83        tridiag_mod-0.3.o out_horiz_mod.o \ 
     84        bmelt_clio_coupl_mod.o bmelt-seuil-profondeur_mod.o  
    8485 
    8586mod_communs_RGL =   util_recovery.o \ 

  • branches/iLoveclim/SOURCES/main3D-0.4-40km.f90

    r40 r91  
    137137  end do time_loop 
    138138 
     139  bmshelfCLIO(:,:,:) = 0d0 
     140 
    139141  if (itracebug.eq.1)  call tracebug('dans main avant call out_recovery ') 
    140142  call out_recovery(iout) 

  • branches/iLoveclim/SOURCES/steps_time_loop.f90

    r77 r91  
    246246     call ablation 
    247247 
    248  
     248     ! oceanic forcing 
     249     !===================== 
     250     call bmeltshelf  ! afq -- icetemp is not called before the 4th dtt, so a first call to bmelt here 
    249251 
    250252     ! SIA dynamic velocities  (deformation + loi de glissement ) 
    … …  
    289291        call icetemp 
    290292        ! subroutines pour le calcul de la fusion basale 
     293        ! afq -- 28-10-2016: there is a reset of bmelt in temp_col in the 
     294        ! current version of the code, so shelf and grounded melt have to  
     295        ! be called here. TODO: remove reset in icetemp routines?? 
    291296        call bmeltshelf 
    292297        call bmelt_grounded 
Note: See TracChangeset for help on using the changeset viewer.