[11403] | 1 | MODULE isfload |
---|
| 2 | !!====================================================================== |
---|
| 3 | !! *** MODULE isfload *** |
---|
| 4 | !! isfload module : compute ice shelf load (needed for the hpg) |
---|
| 5 | !!====================================================================== |
---|
| 6 | !! History : 4.1 ! 2019-09 (P. Mathiot) original code |
---|
| 7 | !!---------------------------------------------------------------------- |
---|
| 8 | |
---|
| 9 | !!---------------------------------------------------------------------- |
---|
| 10 | !! isfload : compute ice shelf load |
---|
| 11 | !!---------------------------------------------------------------------- |
---|
| 12 | |
---|
[12077] | 13 | USE isf_oce, ONLY: cn_isfload, rn_isfload_T, rn_isfload_S ! ice shelf variables |
---|
[11852] | 14 | |
---|
[12068] | 15 | USE dom_oce, ONLY: e3w, gdept, risfdep, mikt ! vertical scale factor |
---|
[11852] | 16 | USE eosbn2 , ONLY: eos ! eos routine |
---|
| 17 | |
---|
[12077] | 18 | USE lib_mpp, ONLY: ctl_stop ! ctl_stop routine |
---|
| 19 | USE in_out_manager ! |
---|
[11403] | 20 | |
---|
| 21 | IMPLICIT NONE |
---|
| 22 | |
---|
| 23 | PRIVATE |
---|
| 24 | |
---|
| 25 | PUBLIC isf_load |
---|
[12340] | 26 | !! * Substitutions |
---|
| 27 | # include "do_loop_substitute.h90" |
---|
[11403] | 28 | |
---|
| 29 | CONTAINS |
---|
| 30 | |
---|
[12068] | 31 | SUBROUTINE isf_load ( Kmm, pisfload ) |
---|
[11403] | 32 | !!-------------------------------------------------------------------- |
---|
| 33 | !! *** SUBROUTINE isf_load *** |
---|
| 34 | !! |
---|
| 35 | !! ** Purpose : compute the ice shelf load |
---|
| 36 | !! |
---|
| 37 | !!-------------------------------------------------------------------- |
---|
| 38 | !!-------------------------- OUT ------------------------------------- |
---|
| 39 | REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: pisfload |
---|
| 40 | !!-------------------------- IN ------------------------------------- |
---|
[12068] | 41 | INTEGER, INTENT(in) :: Kmm ! ocean time level index |
---|
[11403] | 42 | !!---------------------------------------------------------------------- |
---|
| 43 | ! |
---|
| 44 | ! quality test: ice shelf in a stratify/uniform ocean should not drive any flow. |
---|
| 45 | ! the smaller the residual flow is, the better it is. |
---|
| 46 | ! |
---|
| 47 | ! ice shelf cavity |
---|
| 48 | SELECT CASE ( cn_isfload ) |
---|
[12077] | 49 | CASE ( 'uniform' ) |
---|
| 50 | CALL isf_load_uniform ( Kmm, pisfload ) |
---|
[11403] | 51 | CASE DEFAULT |
---|
| 52 | CALL ctl_stop('STOP','method cn_isfload to compute ice shelf load does not exist (isomip), check your namelist') |
---|
| 53 | END SELECT |
---|
| 54 | ! |
---|
| 55 | END SUBROUTINE isf_load |
---|
| 56 | |
---|
[12077] | 57 | SUBROUTINE isf_load_uniform( Kmm, pisfload ) |
---|
[11403] | 58 | !!-------------------------------------------------------------------- |
---|
| 59 | !! *** SUBROUTINE isf_load *** |
---|
| 60 | !! |
---|
| 61 | !! ** Purpose : compute the ice shelf load |
---|
| 62 | !! |
---|
| 63 | !! ** Method : The ice shelf is assumed to be in hydro static equilibrium |
---|
| 64 | !! in water at -1.9 C and 34.4 PSU. Weight of the ice shelf is |
---|
| 65 | !! integrated from top to bottom. |
---|
| 66 | !! |
---|
| 67 | !!-------------------------------------------------------------------- |
---|
| 68 | !!-------------------------- OUT ------------------------------------- |
---|
| 69 | REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: pisfload |
---|
| 70 | !!-------------------------- IN ------------------------------------- |
---|
[12068] | 71 | INTEGER, INTENT(in) :: Kmm ! ocean time level index |
---|
[11403] | 72 | !!-------------------------------------------------------------------- |
---|
| 73 | INTEGER :: ji, jj, jk |
---|
| 74 | INTEGER :: ikt |
---|
| 75 | REAL(wp) :: znad ! |
---|
| 76 | REAL(wp), DIMENSION(jpi,jpj) :: zrhdtop_isf ! water density displaced by the ice shelf (at the interface) |
---|
| 77 | REAL(wp), DIMENSION(jpi,jpj,jpts) :: zts_top ! water properties displaced by the ice shelf |
---|
| 78 | REAL(wp), DIMENSION(jpi,jpj,jpk) :: zrhd ! water density displaced by the ice shelf |
---|
| 79 | !!---------------------------------------------------------------------- |
---|
| 80 | ! |
---|
| 81 | znad = 1._wp !- To use density and not density anomaly |
---|
| 82 | ! |
---|
[12077] | 83 | ! !- assume water displaced by the ice shelf is at T=rn_isfload_T and S=rn_isfload_S (rude) |
---|
| 84 | zts_top(:,:,jp_tem) = rn_isfload_T ; zts_top(:,:,jp_sal) = rn_isfload_S |
---|
[11403] | 85 | ! |
---|
| 86 | DO jk = 1, jpk !- compute density of the water displaced by the ice shelf |
---|
[12068] | 87 | CALL eos( zts_top(:,:,:), gdept(:,:,jk,Kmm), zrhd(:,:,jk) ) |
---|
[11403] | 88 | END DO |
---|
| 89 | ! |
---|
| 90 | ! !- compute rhd at the ice/oce interface (ice shelf side) |
---|
| 91 | CALL eos( zts_top , risfdep, zrhdtop_isf ) |
---|
| 92 | ! |
---|
| 93 | ! !- Surface value + ice shelf gradient |
---|
[11852] | 94 | pisfload(:,:) = 0._wp ! compute pressure due to ice shelf load |
---|
[12340] | 95 | DO_2D_11_11 |
---|
| 96 | ikt = mikt(ji,jj) |
---|
| 97 | ! |
---|
| 98 | IF ( ikt > 1 ) THEN |
---|
[11403] | 99 | ! |
---|
[12340] | 100 | ! top layer of the ice shelf |
---|
| 101 | pisfload(ji,jj) = pisfload(ji,jj) + (znad + zrhd(ji,jj,1) ) * e3w(ji,jj,1,Kmm) |
---|
| 102 | ! |
---|
| 103 | ! core layers of the ice shelf |
---|
| 104 | DO jk = 2, ikt-1 |
---|
| 105 | pisfload(ji,jj) = pisfload(ji,jj) + (2._wp * znad + zrhd(ji,jj,jk-1) + zrhd(ji,jj,jk)) * e3w(ji,jj,jk,Kmm) |
---|
| 106 | END DO |
---|
| 107 | ! |
---|
| 108 | ! deepest part of the ice shelf (between deepest T point and ice/ocean interface |
---|
| 109 | pisfload(ji,jj) = pisfload(ji,jj) + (2._wp * znad + zrhdtop_isf(ji,jj) + zrhd(ji,jj,ikt-1)) & |
---|
| 110 | & * ( risfdep(ji,jj) - gdept(ji,jj,ikt-1,Kmm) ) |
---|
| 111 | ! |
---|
| 112 | END IF |
---|
| 113 | END_2D |
---|
[11403] | 114 | ! |
---|
[12077] | 115 | END SUBROUTINE isf_load_uniform |
---|
[11403] | 116 | |
---|
| 117 | END MODULE isfload |
---|