Changeset 389 for trunk/NEMO/OPA_SRC/stpctl.F90
- Timestamp:
- 2006-03-09T18:22:04+01:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/stpctl.F90
r367 r389 56 56 INTEGER, DIMENSION(3) :: ilocu ! 57 57 INTEGER, DIMENSION(2) :: ilocs ! 58 CHARACTER(len=80) :: clname 58 59 !!---------------------------------------------------------------------- 59 60 !! OPA 9.0 , LOCEAN-IPSL (2005) … … 67 68 WRITE(numout,*) '~~~~~~~' 68 69 ! 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 ) 70 72 ENDIF 71 73 … … 80 82 ! -------------------------- 81 83 IF( lk_dynspg_flt .OR. lk_dynspg_rl ) THEN 82 83 84 ! Solver 85 IF(lwp) WRITE(numsol,9200) kt, niter, res, SQRT(epsr)/eps 84 86 85 86 87 ! Islands (if exist) 88 IF( lk_isl ) CALL isl_stp_ctl( kt, kindic ) 87 89 88 90 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,*) ' ====== *******************************' 103 103 ENDIF 104 CALL dia_wri( kt, kindic ) 105 ENDIF 104 106 ENDIF 105 107 … … 140 142 ENDIF 141 143 kindic = -3 144 142 145 CALL dia_wri( kt, kindic ) 143 146 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.