- Timestamp:
- 2017-09-13T18:46:56+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerdgrft.F90
r8517 r8518 34 34 PUBLIC ice_rdgrft_strength ! called by icerhg_evp 35 35 PUBLIC ice_rdgrft_init ! called by ice_stp 36 PUBLIC ice_rdgrft_alloc ! called by ice_init37 36 38 37 ! Variables shared among ridging subroutines 39 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: asum ! sum of total ice and open water area40 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: aksum ! ratio of area removed to area ridged41 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: athorn ! participation function; fraction of ridging/closing associated w/ category n42 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: hrmin ! minimum ridge thickness43 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: hrmax ! maximum ridge thickness44 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: hraft ! thickness of rafted ice45 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: krdg ! thickness of ridging ice / mean ridge thickness46 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: aridge ! participating ice ridging47 REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) :: araft ! participating ice rafting38 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: asum ! sum of total ice and open water area 39 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: aksum ! ratio of area removed to area ridged 40 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: athorn ! participation function; fraction of ridging/closing associated w/ category n 41 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: hrmin ! minimum ridge thickness 42 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: hrmax ! maximum ridge thickness 43 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: hraft ! thickness of rafted ice 44 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: krdg ! thickness of ridging ice / mean ridge thickness 45 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: aridge ! participating ice ridging 46 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: araft ! participating ice rafting 48 47 ! 49 48 REAL(wp), PARAMETER :: krdgmin = 1.1_wp ! min ridge thickness multiplier … … 52 51 ! 53 52 ! ** namelist (namice_rdgrft) ** 53 LOGICAL :: ln_str_H79 ! ice strength parameterization (Hibler79) 54 REAL(wp) :: rn_pstar ! determines ice strength, Hibler JPO79 55 REAL(wp) :: rn_crhg ! determines changes in ice strength 56 LOGICAL :: ln_str_R75 ! ice strength parameterization (Rothrock75) 57 REAL(wp) :: rn_perdg ! ridging work divided by pot. energy change in ridging 54 58 REAL(wp) :: rn_csrdg ! fraction of shearing energy contributing to ridging 55 59 LOGICAL :: ln_partf_lin ! participation function linear (Thorndike et al. (1975)) … … 82 86 & hrmin(jpi,jpj,jpl) , hraft (jpi,jpj,jpl) , aridge(jpi,jpj,jpl) , & 83 87 & hrmax(jpi,jpj,jpl) , krdg (jpi,jpj,jpl) , araft (jpi,jpj,jpl) , STAT=ice_rdgrft_alloc ) 84 ! 88 89 IF( lk_mpp ) CALL mpp_sum ( ice_rdgrft_alloc ) 85 90 IF( ice_rdgrft_alloc /= 0 ) CALL ctl_warn( 'ice_rdgrft_alloc: failed to allocate arrays' ) 86 91 ! … … 909 914 CALL ctl_stop( 'ice_rdgrft_init: choose one and only one participation function (ln_partf_lin or ln_partf_exp)' ) 910 915 ENDIF 916 ! ! allocate tke arrays 917 IF( ice_rdgrft_alloc() /= 0 ) CALL ctl_stop( 'STOP', 'ice_rdgrft_init : unable to allocate arrays' ) 911 918 ! 912 919 END SUBROUTINE ice_rdgrft_init
Note: See TracChangeset
for help on using the changeset viewer.