New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 11467 for NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icedyn_rhg_evp.F90 – NEMO

Ignore:
Timestamp:
2019-08-22T11:49:08+02:00 (5 years ago)
Author:
andmirek
Message:

Ticket #2197 allocate arrays at the beggining of the run

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/branches/UKMO/dev_r10037_GPU/src/ICE/icedyn_rhg_evp.F90

    r10843 r11467  
    107107      !!              Kimmritz et al., Ocean Modelling 2016 & 2017 
    108108      !!------------------------------------------------------------------- 
     109      USE scoce, ONLY : z1_e1t0 => scr2D1, z1_e2t0 => scr2D2, &                                       ! scale factors 
     110                        zp_delt => scr2D3, &                                                          ! P/delta at T points 
     111                        zbeta => scr2D4, &                                                            ! beta coef from Kimmritz 2017 
     112      ! 
     113                        zdt_m => scr2D5, &                                                            ! (dt / ice-snow_mass) on T points 
     114                        zaU => scr2D6  , zaV => scr2D7, &                                             ! ice fraction on U/V points 
     115                        zmU_t => scr2D8, zmV_t => scr2D9, &                                           ! (ice-snow_mass / dt) on U/V points 
     116                        zmf => scr2D10, &                                                             ! coriolis parameter at T points 
     117                        zTauU_ia => scr2D11, ztauV_ia => scr2D12, &                                   ! ice-atm. stress at U-V points 
     118                        zspgU => scr2D13, zspgV => scr2D14, &                                         ! surface pressure gradient at U/V points 
     119                        v_oceU => scr2D15, u_oceV => scr2D16, v_iceU => scr2D17, u_iceV => scr2D18, & ! ocean/ice u/v component on V/U points                            
     120                        zfU => scr2D19  , zfV => scr2D20, &                                           ! internal stresses 
     121      ! 
     122                        zds => scr2D21, &                                                             ! shear 
     123                        zs1 => scr2D22, zs2 => scr2D23, zs12 => scr2D38, &                            ! stress tensor components 
     124                        zu_ice => scr2D24, zv_ice => scr2D25, zresr => scr2D26, &                     ! check convergence 
     125                        zpice  => scr2D27, &                                                          ! array used for the calculation of ice surface slope: 
     126      !                                                                                               !    ocean surface (ssh_m) if ice is not embedded 
     127      !                                                                                               !    ice top surface if ice is embedded    
     128                        zCorx => scr2D28, zCory => scr2D29, &                                         ! Coriolis stress array 
     129                        ztaux_oi => scr2D30, ztauy_oi => scr2D31, &                                   ! Ocean-to-ice stress array 
     130      ! 
     131                        zswitchU => scr2D32, zswitchV => scr2D33, &                                   ! dummy arrays 
     132                        zmaskU => scr2D34, zmaskV => scr2D35, &                                       ! mask for ice presence 
     133                        zfmask => scr2D36, zwf => scr2D37                                             ! mask at F points for the ice 
     134 
     135                         
     136 
    109137      INTEGER                 , INTENT(in   ) ::   kt                                    ! time step 
    110138      REAL(wp), DIMENSION(:,:), INTENT(inout) ::   pstress1_i, pstress2_i, pstress12_i   ! 
     
    127155      REAL(wp) ::   zshear, zdum1, zdum2 
    128156      ! 
    129       REAL(wp), DIMENSION(jpi,jpj) ::   z1_e1t0, z1_e2t0                ! scale factors 
    130       REAL(wp), DIMENSION(jpi,jpj) ::   zp_delt                         ! P/delta at T points 
    131       REAL(wp), DIMENSION(jpi,jpj) ::   zbeta                           ! beta coef from Kimmritz 2017 
    132       ! 
    133       REAL(wp), DIMENSION(jpi,jpj) ::   zdt_m                           ! (dt / ice-snow_mass) on T points 
    134       REAL(wp), DIMENSION(jpi,jpj) ::   zaU   , zaV                     ! ice fraction on U/V points 
    135       REAL(wp), DIMENSION(jpi,jpj) ::   zmU_t, zmV_t                    ! (ice-snow_mass / dt) on U/V points 
    136       REAL(wp), DIMENSION(jpi,jpj) ::   zmf                             ! coriolis parameter at T points 
    137       REAL(wp), DIMENSION(jpi,jpj) ::   zTauU_ia , ztauV_ia             ! ice-atm. stress at U-V points 
    138       REAL(wp), DIMENSION(jpi,jpj) ::   zspgU , zspgV                   ! surface pressure gradient at U/V points 
    139       REAL(wp), DIMENSION(jpi,jpj) ::   v_oceU, u_oceV, v_iceU, u_iceV  ! ocean/ice u/v component on V/U points                            
    140       REAL(wp), DIMENSION(jpi,jpj) ::   zfU   , zfV                     ! internal stresses 
    141       ! 
    142       REAL(wp), DIMENSION(jpi,jpj) ::   zds                             ! shear 
    143       REAL(wp), DIMENSION(jpi,jpj) ::   zs1, zs2, zs12                  ! stress tensor components 
    144       REAL(wp), DIMENSION(jpi,jpj) ::   zu_ice, zv_ice, zresr           ! check convergence 
    145       REAL(wp), DIMENSION(jpi,jpj) ::   zpice                           ! array used for the calculation of ice surface slope: 
    146       !                                                                 !    ocean surface (ssh_m) if ice is not embedded 
    147       !                                                                 !    ice top surface if ice is embedded    
    148       REAL(wp), DIMENSION(jpi,jpj) ::   zCorx, zCory                    ! Coriolis stress array 
    149       REAL(wp), DIMENSION(jpi,jpj) ::   ztaux_oi, ztauy_oi              ! Ocean-to-ice stress array 
    150       ! 
    151       REAL(wp), DIMENSION(jpi,jpj) ::   zswitchU, zswitchV              ! dummy arrays 
    152       REAL(wp), DIMENSION(jpi,jpj) ::   zmaskU, zmaskV                  ! mask for ice presence 
    153       REAL(wp), DIMENSION(jpi,jpj) ::   zfmask, zwf                     ! mask at F points for the ice 
    154  
    155157      REAL(wp), PARAMETER          ::   zepsi  = 1.0e-20_wp             ! tolerance parameter 
    156158      REAL(wp), PARAMETER          ::   zmmin  = 1._wp                  ! ice mass (kg/m2) below which ice velocity equals ocean velocity 
Note: See TracChangeset for help on using the changeset viewer.