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 389 for trunk/NEMO/OPA_SRC/stpctl.F90 – NEMO

Ignore:
Timestamp:
2006-03-09T18:22:04+01:00 (18 years ago)
Author:
opalod
Message:

RB:nemo_v1_update_038: first integration of Agrif :

  • configuration parameters are just integer when agrif is used
  • add call to agrif routines with key_agrif
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/stpctl.F90

    r367 r389  
    5656      INTEGER, DIMENSION(3) ::   ilocu      !  
    5757      INTEGER, DIMENSION(2) ::   ilocs      !  
     58      CHARACTER(len=80) :: clname 
    5859      !!---------------------------------------------------------------------- 
    5960      !!  OPA 9.0 , LOCEAN-IPSL (2005)  
     
    6768         WRITE(numout,*) '~~~~~~~' 
    6869         ! open time.step file 
    69          CALL ctlopn( numstp, 'time.step', 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', 1, numout, lwp, 1 ) 
     70         clname = 'time.step' 
     71         CALL ctlopn( numstp, clname, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', 1, numout, lwp, 1 ) 
    7072      ENDIF 
    7173 
     
    8082      ! -------------------------- 
    8183      IF( lk_dynspg_flt .OR. lk_dynspg_rl ) THEN 
    82          ! Solver 
    83          IF(lwp) WRITE(numsol,9200) kt, niter, res, SQRT(epsr)/eps 
     84      ! Solver 
     85      IF(lwp) WRITE(numsol,9200) kt, niter, res, SQRT(epsr)/eps 
    8486 
    85          ! Islands (if exist) 
    86          IF( lk_isl )   CALL isl_stp_ctl( kt, kindic ) 
     87      ! Islands (if exist) 
     88      IF( lk_isl )   CALL isl_stp_ctl( kt, kindic ) 
    8789 
    8890 
    89          ! Output in numwso and numwvo IF kindic<0 
    90          ! --------------------------------------- 
    91          !    (i.e. problem for the solver) 
    92          IF( kindic < 0 ) THEN 
    93             IF(lwp) THEN 
    94                WRITE(numout,*) ' stpctl: the elliptic solver DO not converge or explode' 
    95                WRITE(numout,*) ' ====== ' 
    96                WRITE(numout,9200) kt, niter, res, sqrt(epsr)/eps 
    97                WRITE(numout,*) 
    98                WRITE(numout,*) ' stpctl: output of last fields in numwso' 
    99                WRITE(numout,*) '                                  numwvo' 
    100                WRITE(numout,*) ' ======  *******************************' 
    101             ENDIF 
    102             CALL dia_wri( kt, kindic ) 
     91      ! Output in numwso and numwvo IF kindic<0 
     92      ! --------------------------------------- 
     93      !    (i.e. problem for the solver) 
     94      IF( kindic < 0 ) THEN 
     95         IF(lwp) THEN 
     96            WRITE(numout,*) ' stpctl: the elliptic solver DO not converge or explode' 
     97            WRITE(numout,*) ' ====== ' 
     98            WRITE(numout,9200) kt, niter, res, sqrt(epsr)/eps 
     99            WRITE(numout,*) 
     100            WRITE(numout,*) ' stpctl: output of last fields in numwso' 
     101            WRITE(numout,*) '                                  numwvo' 
     102            WRITE(numout,*) ' ======  *******************************' 
    103103         ENDIF 
     104         CALL dia_wri( kt, kindic ) 
     105      ENDIF 
    104106      ENDIF 
    105107 
     
    140142         ENDIF 
    141143         kindic  = -3 
     144 
    142145         CALL dia_wri( kt, kindic ) 
    143146      ENDIF 
Note: See TracChangeset for help on using the changeset viewer.