Ignore:
Timestamp:
2019-08-08T16:02:49+02:00 (15 months ago)
Author:
mathiot
Message:

ENHANCE-02_ISF_nemo : add UKESM ice sheet coupling method (ticket #2142)

File:
1 edited

Legend:

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

    r11403 r11423  
    1616   USE isf            ! Ice shelf variable 
    1717   USE isfutils       ! 
    18    USE iscplini       ! Ice sheet coupling 
    1918   ! 
    2019   USE in_out_manager ! I/O manager 
     
    5655      ! 
    5756      ! ice sheet coupling case 
    58       !CALL tra_isf_cpl() 
     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) 
    5961      ! 
    6062      IF( ln_timing )   CALL timing_stop('tra_isf') 
     
    7678      REAL(wp), DIMENSION(jpi,jpj)     , INTENT(in   ) :: phtbl, pfrac 
    7779      REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in   ) :: ptsc , ptsc_b 
    78       REAL(wp), DIMENSION(jpi,jpj) :: ztc 
    7980      !!---------------------------------------------------------------------- 
    8081      INTEGER                      :: ji,jj,jk  ! loop index    
     
    107108   END SUBROUTINE tra_isf_mlt 
    108109   ! 
    109    !SUBROUTINE tra_isf_cpl 
    110        ! 
    111 !      !---------------------------------------- 
    112 !      !        Ice Sheet coupling imbalance correction to have conservation 
    113 !      !---------------------------------------- 
    114 !      ! 
    115 !      IF( ln_iscpl .AND. ln_iscpl_hsb) THEN         ! input of heat and salt due to river runoff  
    116 !         DO jk = 1,jpk 
    117 !            DO jj = 2, jpj  
    118 !               DO ji = fs_2, fs_jpim1 
    119 !                  zdep = 1._wp / e3t_n(ji,jj,jk)  
    120 !                  tsa(ji,jj,jk,jp_tem) = tsa(ji,jj,jk,jp_tem) - htsc_iscpl(ji,jj,jk,jp_tem) * zdep 
    121 !                  tsa(ji,jj,jk,jp_sal) = tsa(ji,jj,jk,jp_sal) - htsc_iscpl(ji,jj,jk,jp_sal) * zdep   
    122 !               END DO   
    123 !            END DO   
    124 !         END DO 
    125 !      ENDIF 
    126 !      ! 
    127 !   END SUBROUTINE tra_isf_cpl 
     110   SUBROUTINE tra_isf_cpl( ptsc, ptsa ) 
     111      !!---------------------------------------------------------------------- 
     112      !!                  ***  ROUTINE tra_isf_cpl  *** 
     113      !! 
     114      !! *** Action :: Update tsa with the ice shelf coupling trend  
     115      !! 
     116      !!---------------------------------------------------------------------- 
     117      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(inout) :: ptsa 
     118      !!---------------------------------------------------------------------- 
     119      REAL(wp), DIMENSION(jpi,jpj,jpk,jpts), INTENT(in   ) :: ptsc  !( >0 out ) 
     120      !!---------------------------------------------------------------------- 
     121      INTEGER :: jk 
     122      !!---------------------------------------------------------------------- 
     123      ! 
     124      DO jk = 1,jpk 
     125         ptsa(:,:,jk,jp_tem) = ptsa(:,:,jk,jp_tem) - ptsc(:,:,jk,jp_tem) * r1_e1e2t(:,:) / e3t_n(:,:,jk) 
     126         ptsa(:,:,jk,jp_sal) = ptsa(:,:,jk,jp_sal) - ptsc(:,:,jk,jp_sal) * r1_e1e2t(:,:) / e3t_n(:,:,jk) 
     127      END DO 
     128      ! 
     129   END SUBROUTINE tra_isf_cpl 
    128130   ! 
    129131END MODULE traisf 
Note: See TracChangeset for help on using the changeset viewer.