- Timestamp:
- 2019-08-08T16:02:49+02:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/ENHANCE-02_ISF_nemo/src/OCE/TRA/traisf.F90
r11403 r11423 16 16 USE isf ! Ice shelf variable 17 17 USE isfutils ! 18 USE iscplini ! Ice sheet coupling19 18 ! 20 19 USE in_out_manager ! I/O manager … … 56 55 ! 57 56 ! 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) 59 61 ! 60 62 IF( ln_timing ) CALL timing_stop('tra_isf') … … 76 78 REAL(wp), DIMENSION(jpi,jpj) , INTENT(in ) :: phtbl, pfrac 77 79 REAL(wp), DIMENSION(jpi,jpj,jpts), INTENT(in ) :: ptsc , ptsc_b 78 REAL(wp), DIMENSION(jpi,jpj) :: ztc79 80 !!---------------------------------------------------------------------- 80 81 INTEGER :: ji,jj,jk ! loop index … … 107 108 END SUBROUTINE tra_isf_mlt 108 109 ! 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 128 130 ! 129 131 END MODULE traisf
Note: See TracChangeset
for help on using the changeset viewer.