Changeset 498 for trunk/NEMO/OFF_SRC/TRA/traqsr.F90
- Timestamp:
- 2006-09-12T13:07:10+02:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OFF_SRC/TRA/traqsr.F90
r343 r498 39 39 !! ** Method : The profile of solar radiation within the ocean is set 40 40 !! from two length scale of penetration (xsr1,xsr2) and a ratio 41 !! (rabs). These parameters are read in the nam qsr namelist. The41 !! (rabs). These parameters are read in the nam_qsr namelist. The 42 42 !! default values correspond to clear water (type I in Jerlov' 43 43 !! (1968) classification. … … 54 54 !! * Local declarations 55 55 56 NAMELIST/nam qsr/ ln_traqsr, xsi156 NAMELIST/nam_qsr/ ln_traqsr, xsi1 57 57 !!---------------------------------------------------------------------- 58 58 59 ! Read Namelist nam qsr : ratio and length of penetration59 ! Read Namelist nam_qsr : ratio and length of penetration 60 60 ! -------------------- 61 61 REWIND ( numnam ) 62 READ ( numnam, nam qsr )62 READ ( numnam, nam_qsr ) 63 63 64 64 ! Parameter control and print … … 69 69 WRITE(numout,*) 'tra_qsr_init : penetration of the surface solar radiation' 70 70 WRITE(numout,*) '~~~~~~~~~~~~' 71 WRITE(numout,*) ' Namelist nam qsr : set the parameter of penetration'71 WRITE(numout,*) ' Namelist nam_qsr : set the parameter of penetration' 72 72 WRITE(numout,*) ' first depth of extinction xsi1 = ',xsi1 73 73 WRITE(numout,*) ' ' … … 81 81 ENDIF 82 82 83 IF( xsi1 < 0.e0 ) THEN 84 IF(lwp) WRITE(numout,cform_err) 85 IF(lwp) WRITE(numout,*) ' 0<xsi1 not satisfied' 86 nstop = nstop + 1 87 ENDIF 83 IF( xsi1 < 0.e0 ) & 84 CALL ctl_stop( ' 0<xsi1 not satisfied' ) 88 85 89 86
Note: See TracChangeset
for help on using the changeset viewer.