Ignore:
Timestamp:
2019-09-10T11:49:13+02:00 (15 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF: fix issue with ice sheet coupling and conservation + other minor changes (ticket #2142)

File:
1 edited

Legend:

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

    r11423 r11521  
    4848      IF( ln_timing )   CALL timing_start('tra_sbc') 
    4949      ! 
    50       ! cavity case 
    51       IF ( ln_isfcav_mlt ) CALL tra_isf_mlt(misfkt_cav, misfkb_cav, rhisf_tbl_cav, rfrac_tbl_cav, risf_cav_tsc, risf_cav_tsc_b, tsa) 
    52       ! 
    53       ! parametrisation case 
    54       IF ( ln_isfpar_mlt ) CALL tra_isf_mlt(misfkt_par, misfkb_par, rhisf_tbl_par, rfrac_tbl_par, risf_par_tsc, risf_par_tsc_b, tsa) 
     50      IF (ln_isf) THEN 
     51         ! 
     52         ! cavity case 
     53         IF ( ln_isfcav_mlt ) CALL tra_isf_mlt(misfkt_cav, misfkb_cav, rhisf_tbl_cav, rfrac_tbl_cav, risf_cav_tsc, risf_cav_tsc_b, tsa) 
     54         ! 
     55         ! parametrisation case 
     56         IF ( ln_isfpar_mlt ) CALL tra_isf_mlt(misfkt_par, misfkb_par, rhisf_tbl_par, rfrac_tbl_par, risf_par_tsc, risf_par_tsc_b, tsa) 
     57         ! 
     58      END IF 
    5559      ! 
    5660      ! ice sheet coupling case 
    57       IF ( ll_isfcpl .AND. kt == nit000 ) CALL tra_isf_cpl(risfcpl_tsc, tsa) 
    58       ! 
    59       ! ice sheet coupling case ( 
    60       IF ( ll_isfcpl_cons ) CALL tra_isf_cpl(risfcpl_cons_tsc, tsa) 
     61      IF ( ll_isfcpl ) THEN 
     62         ! 
     63         IF ( kt == nit000 ) CALL tra_isf_cpl(risfcpl_tsc, tsa) 
     64         ! 
     65         ! ensure 0 trend due to unconservation of the ice shelf coupling 
     66         IF ( ll_isfcpl_cons ) CALL tra_isf_cpl(risfcpl_cons_tsc, tsa) 
     67         ! 
     68      END IF 
    6169      ! 
    6270      IF( ln_timing )   CALL timing_stop('tra_isf') 
     
    95103            ! 
    96104            ! level fully include in the ice shelf boundary layer 
    97             ! sign - because fwf sign of evapo (rnf sign of precip) 
    98105            DO jk = ikt, ikb - 1 
    99106               pts(ji,jj,jk,jp_tem) = pts(ji,jj,jk,jp_tem) + ztc(ji,jj) 
     
    117124      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(inout) :: ptsa 
    118125      !!---------------------------------------------------------------------- 
    119       REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) :: ptsc  !( >0 out ) 
     126      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) :: ptsc 
    120127      !!---------------------------------------------------------------------- 
    121128      INTEGER :: jk 
Note: See TracChangeset for help on using the changeset viewer.