Changeset 61 for trunk


Ignore:
Timestamp:
10/26/12 05:08:09 (11 years ago)
Author:
rblod
Message:

Ensure correct restart with agrif activated (need to allocate gcx)

Location:
trunk/NEMOGCM/NEMO/OPA_SRC
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90

    r7 r61  
    2424   !!   istate_uvg    : initial velocity in geostropic balance 
    2525   !!---------------------------------------------------------------------- 
     26   USE sol_oce 
    2627   USE oce             ! ocean dynamics and active tracers  
    2728   USE dom_oce         ! ocean space and time domain  
     
    138139      IF( lk_agrif ) THEN                  ! read free surface arrays in restart file 
    139140         IF( ln_rstart ) THEN 
    140             IF( lk_dynspg_flt )   CALL flt_rst( nit000, 'READ' )      ! read or initialize the following fields 
    141             !                                                         ! gcx, gcxb for agrif_opa_init 
    142          ENDIF                                                        ! explicit case not coded yet with AGRIF 
     141            IF( lk_dynspg_flt )  THEN      ! read or initialize the following fields 
     142               !                           ! gcx, gcxb for agrif_opa_init 
     143               IF( sol_oce_alloc()  > 0 )   CALL ctl_stop('agrif sol_oce_alloc: allocation of arrays failed') 
     144               CALL flt_rst( nit000, 'READ' ) 
     145            ENDIF   
     146         ENDIF                             ! explicit case not coded yet with AGRIF 
    143147      ENDIF 
    144148      ! 
  • trunk/NEMOGCM/NEMO/OPA_SRC/SOL/solver.F90

    r7 r61  
    8383 
    8484      !                              ! allocate solver arrays 
    85       IF( sol_oce_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'solver_init : unable to allocate sol_oce arrays' ) 
     85      IF( .NOT. lk_agrif .OR. .NOT. ln_rstart) THEN 
     86         IF( sol_oce_alloc() /= 0 )   CALL ctl_stop( 'STOP', 'solver_init : unable to allocate sol_oce arrays' ) 
     87      ENDIF  
    8688 
    8789      SELECT CASE( nn_solv )          !* parameter check 
Note: See TracChangeset for help on using the changeset viewer.