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/OCE/DYN/dynspg.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/OCE/DYN/dynspg.F90

    r10843 r11467  
    7171      !!             period is used to prevent the divergence of odd and even time step. 
    7272      !!---------------------------------------------------------------------- 
     73      USE scoce, ONLY : zpice => scr2D1, ztrdu => scr1, ztrdv => scr2 
    7374      INTEGER, INTENT(in) ::   kt   ! ocean time-step index 
    7475      ! 
    7576      INTEGER  ::   ji, jj, jk                   ! dummy loop indices 
    7677      REAL(wp) ::   z2dt, zg_2, zintp, zgrau0r, zld   ! local scalars 
    77       REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zpice 
    78       REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   ztrdu, ztrdv 
    7978      !!---------------------------------------------------------------------- 
    8079      ! 
     
    8281      ! 
    8382      IF( l_trddyn )   THEN                      ! temporary save of ta and sa trends 
    84          ALLOCATE( ztrdu(jpi,jpj,jpk) , ztrdv(jpi,jpj,jpk) )  
    8583         ztrdu(:,:,:) = ua(:,:,:) 
    8684         ztrdv(:,:,:) = va(:,:,:) 
     
    134132         ! 
    135133         IF( ln_ice_embd ) THEN              !== embedded sea ice: Pressure gradient due to snow-ice mass ==! 
    136             ALLOCATE( zpice(jpi,jpj) ) 
    137134            zintp = REAL( MOD( kt-1, nn_fsbc ) ) / REAL( nn_fsbc ) 
    138135            zgrau0r     = - grav * r1_rau0 
     
    144141               END DO 
    145142            END DO 
    146             DEALLOCATE( zpice )          
    147143         ENDIF 
    148144         ! 
     
    169165         ztrdv(:,:,:) = va(:,:,:) - ztrdv(:,:,:) 
    170166         CALL trd_dyn( ztrdu, ztrdv, jpdyn_spg, kt ) 
    171          DEALLOCATE( ztrdu , ztrdv )  
    172167      ENDIF 
    173168      !                                      ! print mean trends (used for debugging) 
Note: See TracChangeset for help on using the changeset viewer.