Changeset 625
- Timestamp:
- 2007-02-27T16:51:00+01:00 (18 years ago)
- Location:
- trunk/NEMO/TOP_SRC
- Files:
-
- 8 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/SMS/trcini.cfc.h90
r433 r625 4 4 !! * Module variables 5 5 INTEGER :: & 6 inum = 16! unit number6 inum ! unit number 7 7 CHARACTER (len=34) :: & 8 8 clname = 'cfc1112.atm' ! ??? … … 78 78 IF(lwp) WRITE(numout,*) 'read of formatted file cfc1112atm' 79 79 80 OPEN( UNIT = inum, FORM = 'FORMATTED',FILE = clname, STATUS = 'OLD') 80 CALL ctlopn( inum, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', & 81 & 1, numout, .FALSE., 1 ) 81 82 REWIND(inum) 82 83 -
trunk/NEMO/TOP_SRC/SMS/trcini.pisces.h90
r617 r625 54 54 55 55 INTEGER :: numriv,numdust,numbath,numdep 56 INTEGER :: numlight = 4956 INTEGER :: numlight 57 57 58 58 #if defined key_trc_kriest … … 598 598 !! A PRECOMPUTED ARRAY CORRESPONDING TO THE ATTENUATION COEFFICIENT 599 599 600 OPEN( numlight, file = 'kRGB61.txt', form = 'formatted') 600 CALL ctlopn( numlight, 'kRGB61.txt', 'OLD', 'FORMATTED', 'SEQUENTIAL', & 601 & 1, numout, .TRUE., 1 ) 601 602 DO ichl = 1,61 602 603 READ(numlight,*) ztoto,(xkrgb(iband,ichl),iband = 1,3) -
trunk/NEMO/TOP_SRC/SMS/trclsm.cfc.h90
r348 r625 26 26 clname = 'namelist.trc.sms' 27 27 INTEGER :: & 28 numnat = 8028 numnat 29 29 30 30 NAMELIST/namdates/ndate_beg, nyear_res … … 41 41 ENDIF 42 42 43 OPEN( UNIT=numnat, FILE=clname, FORM='FORMATTED', STATUS='OLD' ) 43 CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', & 44 & 1, numout, .FALSE., 1 ) 44 45 READ( numnat , namdates ) 45 46 -
trunk/NEMO/TOP_SRC/SMS/trclsm.lobster1.h90
r433 r625 58 58 59 59 60 numnat=8061 60 clname ='namelist.trc.sms' 62 OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old') 61 CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', & 62 & 1, numout, .FALSE., 1 ) 63 63 64 64 ! 1.4 namelist natbio : biological parameters -
trunk/NEMO/TOP_SRC/SMS/trclsm.pisces.h90
r617 r625 66 66 ENDIF 67 67 68 numnat=8069 68 clname ='namelist.trc.sms' 70 OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old')71 69 CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', & 70 & 1, numout, .FALSE., 1 ) 72 71 73 72 ! 1 Namelist natext : -
trunk/NEMO/TOP_SRC/prtctl_trc.F90
r426 r625 15 15 16 16 !! * Module declaration 17 INTEGER, DIMENSION(:), ALLOCATABLE :: numid_trc ! logical unit 17 18 INTEGER, DIMENSION(:), ALLOCATABLE, SAVE :: & !: 18 19 nlditl , nldjtl , & !: first, last indoor index for each i-domain … … 81 82 82 83 !! * Local declarations 83 INTEGER :: overlap, numid, jn, js, sind, eind, kdir84 INTEGER :: overlap, jn, js, sind, eind, kdir, j_id 84 85 REAL(wp) :: zsum, zvctl 85 86 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zmask, ztab3d … … 119 120 DO js = sind, eind 120 121 121 numid = 90 + js122 122 ! Set logical unit 123 j_id = js - narea + 1 123 124 ! Set indices for the SUM control 124 125 IF( .NOT. lsp_area ) THEN … … 152 153 zsum = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) & 153 154 & *zmask(nictls:nictle,njctls:njctle,1:kdir) ) 154 WRITE(numid ,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl155 WRITE(numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum-zvctl 155 156 tra_ctl(jn,js) = zsum 156 157 ENDDO … … 160 161 zsum = SUM( ztab3d(nictls:nictle,njctls:njctle,1:kdir) & 161 162 & * zmask(nictls:nictle,njctls:njctle,1:kdir) ) 162 WRITE(numid ,FMT="(3x,a,' : ',D23.16)") cl(jn), zsum163 WRITE(numid_trc(j_id),FMT="(3x,a,' : ',D23.16)") cl(jn), zsum 163 164 END DO 164 165 ENDIF … … 185 186 186 187 !! * Local declarations 187 INTEGER :: numid, js, sind, eind188 INTEGER :: js, sind, eind, j_id 188 189 !!---------------------------------------------------------------------- 189 190 … … 200 201 ! Loop over each sub-domain, i.e. number of processors ijsplt 201 202 DO js = sind, eind 202 numid = 90 + js203 WRITE(numid ,*)clinfo203 j_id = js - narea + 1 204 WRITE(numid_trc(j_id),*)clinfo 204 205 ENDDO 205 206 … … 218 219 !!---------------------------------------------------------------------- 219 220 !! * Local declarations 220 INTEGER :: js, numid, sind, eind221 INTEGER :: js, sind, eind, j_id 221 222 CHARACTER (len=31) :: clfile_out 222 223 CHARACTER (len=27) :: clb_name … … 268 269 ENDIF 269 270 271 ALLOCATE(numid_trc(eind-sind+1)) 272 270 273 DO js = sind, eind 271 numid = 90 + js272 274 WRITE(clfile_out,FMT=clb_name) js-1 273 OPEN ( UNIT=numid, FILE=TRIM(clfile_out),FORM='FORMATTED' ) 274 WRITE(numid,*) 275 WRITE(numid,*) ' L O D Y C - I P S L' 276 WRITE(numid,*) ' O P A model' 277 WRITE(numid,*) ' Ocean General Circulation Model' 278 WRITE(numid,*) ' version OPA 9.0 (2005) ' 279 WRITE(numid,*) 280 WRITE(numid,*) ' PROC number: ', js 281 WRITE(numid,*) 282 WRITE(numid,FMT="(19x,a20)")cl_run 275 j_id = js -narea + 1 276 CALL ctlopn( numid_trc(j_id), clfile_out, 'UNKNOWN', 'FORMATTED', 'SEQUENTIAL', & 277 & 1, numout, .FALSE., 1 ) 278 WRITE(numid_trc(j_id),*) 279 WRITE(numid_trc(j_id),*) ' L O D Y C - I P S L' 280 WRITE(numid_trc(j_id),*) ' O P A model' 281 WRITE(numid_trc(j_id),*) ' Ocean General Circulation Model' 282 WRITE(numid_trc(j_id),*) ' version OPA 9.0 (2005) ' 283 WRITE(numid_trc(j_id),*) 284 WRITE(numid_trc(j_id),*) ' PROC number: ', js 285 WRITE(numid_trc(j_id),*) 286 WRITE(numid_trc(j_id),FMT="(19x,a20)")cl_run 283 287 284 288 ! Print the SUM control indices … … 296 300 ENDIF 297 301 ENDIF 298 WRITE(numid ,*)299 WRITE(numid ,*) 'prt_tra_ctl : Sum control indices'300 WRITE(numid ,*) '~~~~~~~'301 WRITE(numid ,*)302 WRITE(numid ,9000)' nlej = ', nlejtl(js), ' '303 WRITE(numid ,9000)' ------------- njctle = ', njctle, ' -------------'304 WRITE(numid ,9001)' | |'305 WRITE(numid ,9001)' | |'306 WRITE(numid ,9001)' | |'307 WRITE(numid ,9002)' nictls = ', nictls, ' nictle = ', nictle308 WRITE(numid ,9002)' nldi = ', nlditl(js), ' nlei = ', nleitl(js)309 WRITE(numid ,9001)' | |'310 WRITE(numid ,9001)' | |'311 WRITE(numid ,9001)' | |'312 WRITE(numid ,9004)' njmpp = ',njmpptl(js),' ------------- njctls = ', njctls, ' -------------'313 WRITE(numid ,9003)' nimpp = ', nimpptl(js), ' nldj = ', nldjtl(js), ' '314 WRITE(numid ,*)315 WRITE(numid ,*)302 WRITE(numid_trc(j_id),*) 303 WRITE(numid_trc(j_id),*) 'prt_tra_ctl : Sum control indices' 304 WRITE(numid_trc(j_id),*) '~~~~~~~' 305 WRITE(numid_trc(j_id),*) 306 WRITE(numid_trc(j_id),9000)' nlej = ', nlejtl(js), ' ' 307 WRITE(numid_trc(j_id),9000)' ------------- njctle = ', njctle, ' -------------' 308 WRITE(numid_trc(j_id),9001)' | |' 309 WRITE(numid_trc(j_id),9001)' | |' 310 WRITE(numid_trc(j_id),9001)' | |' 311 WRITE(numid_trc(j_id),9002)' nictls = ', nictls, ' nictle = ', nictle 312 WRITE(numid_trc(j_id),9002)' nldi = ', nlditl(js), ' nlei = ', nleitl(js) 313 WRITE(numid_trc(j_id),9001)' | |' 314 WRITE(numid_trc(j_id),9001)' | |' 315 WRITE(numid_trc(j_id),9001)' | |' 316 WRITE(numid_trc(j_id),9004)' njmpp = ',njmpptl(js),' ------------- njctls = ', njctls, ' -------------' 317 WRITE(numid_trc(j_id),9003)' nimpp = ', nimpptl(js), ' nldj = ', nldjtl(js), ' ' 318 WRITE(numid_trc(j_id),*) 319 WRITE(numid_trc(j_id),*) 316 320 317 321 9000 FORMAT(a41,i4.4,a14) -
trunk/NEMO/TOP_SRC/trcctl.F90
r340 r625 55 55 istop = 0 56 56 57 ! 1. LOGICAL UNIT initialization for specifi! files for passive tracer 58 ! -------------------------------------------------------------------- 59 ! nutwrs : OUTPUT for passive tracer restart UNIT (always used) 60 ! nutrst : restart FILE INPUT UNIT (lrsttr=.TRUE.) 61 ! nutini(jptra) : UNIT for initial FILE for tracer 62 63 nutwrs = 72 64 nutrst = 73 65 66 ! 2. restart for passive tracer (input) 57 ! 1. restart for passive tracer (input) 67 58 ! ----------------------------- 68 59 … … 83 74 IF(lwp) WRITE(numout,*) ' ' 84 75 85 ! 3. OPEN FILES for initial tracer value76 ! 2. OPEN FILES for initial tracer value 86 77 ! -------------------------------------- 87 78 DO jn=1,jptra … … 101 92 ENDIF 102 93 103 ! 4. Don't USE non penetrative convective mixing option94 ! 3. Don't USE non penetrative convective mixing option 104 95 ! it's not implemented for passive tracer 105 96 ! ----------------------------------------------------- … … 114 105 ENDIF 115 106 116 ! 5. transport scheme option107 ! 4. transport scheme option 117 108 ! -------------------------- 118 109 … … 121 112 122 113 123 ! 6. SMS model114 ! 5. SMS model 124 115 ! --------------------------------------------- 125 116 -
trunk/NEMO/TOP_SRC/trclec.F90
r433 r625 89 89 ENDIF 90 90 91 numnat=8092 REWIND (numnat)93 94 91 clname='namelist.passivetrc' 95 OPEN( numnat, FILE= clname, FORM='formatted', STATUS = 'old') 92 CALL ctlopn( numnat, clname, 'OLD', 'FORMATTED', 'SEQUENTIAL', & 93 & 1, numout, .FALSE., 1 ) 94 96 95 97 96 !! 1., 2. & 3. initialization with namelist files
Note: See TracChangeset
for help on using the changeset viewer.