Ignore:
Timestamp:
2019-11-19T18:30:57+01:00 (11 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo: add comments, improve memory usage of ln_isfcpl_cons option, fix issue in ISOMIP+ configuration

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/ISF/isfpar.F90

    r11876 r11931  
    1414   !!   isfpar       : compute ice shelf melt using a prametrisation of ice shelf cavities 
    1515   !!---------------------------------------------------------------------- 
    16    !USE oce            ! ocean dynamics and tracers 
    1716   USE isf            ! ice shelf 
    1817   ! 
     
    2322   USE isfdiags , ONLY: isf_diags_flx  ! ice shelf diags subroutine 
    2423   ! 
    25    USE dom_oce  , ONLY: jpi,jpj, bathy ! ocean space and time domain 
     24   USE dom_oce  , ONLY: bathy          ! ocean space and time domain 
     25   USE par_oce  , ONLY: jpi,jpj        ! ocean space and time domain 
    2626   USE phycst   , ONLY: r1_rau0_rcp    ! physical constants 
    2727   ! 
     
    6262      CALL isfpar_mlt( kt, zqhc, zqoce, pqfwf  ) 
    6363      ! 
    64       ! compute heat and water flux  (change signe directly in the melt subroutine) 
     64      ! compute heat and water flux ( > 0 out ) 
    6565      pqfwf(:,:) = pqfwf(:,:) * mskisf_par(:,:) 
    6666      zqoce(:,:) = zqoce(:,:) * mskisf_par(:,:) 
    6767      zqhc (:,:) = zqhc(:,:)  * mskisf_par(:,:) 
    6868      ! 
    69       ! compute heat content flux 
     69      ! compute heat content flux ( > 0 out ) 
    7070      zqlat(:,:) = pqfwf(:,:) * rLfusisf    ! 2d latent heat flux (W/m2) 
    7171      ! 
    72       ! total heat flux 
     72      ! total heat flux ( > 0 out ) 
    7373      zqh(:,:) = ( zqhc (:,:) + zqoce(:,:) ) 
    7474      ! 
     
    107107      CALL isf_alloc_par() 
    108108      ! 
    109       ! par 
     109      ! initialisation 
    110110      misfkt_par(:,:)     = 1         ; misfkb_par(:,:)       = 1          
    111111      rhisf_tbl_par(:,:)  = 1e-20     ; rfrac_tbl_par(:,:)    = 0.0_wp 
Note: See TracChangeset for help on using the changeset viewer.