/[MITgcm]/MITgcm_contrib/darwin2/pkg/darwin/wavebands_init_vari.F
ViewVC logotype

Diff of /MITgcm_contrib/darwin2/pkg/darwin/wavebands_init_vari.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by jahn, Wed Apr 13 18:56:25 2011 UTC revision 1.2 by jahn, Wed Dec 4 21:19:12 2013 UTC
# Line 36  CEOP Line 36  CEOP
36    
37  #ifdef WAVEBANDS  #ifdef WAVEBANDS
38    
39    C     !FUNCTIONS:
40          LOGICAL MASTER_CPU_THREAD
41          EXTERNAL MASTER_CPU_THREAD
42    
43  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
44  C     == Local variables ==  C     == Local variables ==
45         _RL outfile(npmax,tlam)         _RL outfile(npmax,tlam)
# Line 86  c ANNA create output files for WAVEBANDS Line 90  c ANNA create output files for WAVEBANDS
90  c ANNA see format_helps folder for more info on file structures  c ANNA see format_helps folder for more info on file structures
91    
92  c file for ap_types assigned via coin flips  c file for ap_types assigned via coin flips
93            IF (MASTER_CPU_THREAD(myThid)) THEN
94           CALL MDSFINDUNIT( IniUnit1, mythid )           CALL MDSFINDUNIT( IniUnit1, mythid )
95           open(IniUnit1,file='p-ini-char-aptype.dat',           open(IniUnit1,file='p-ini-char-aptype.dat',
96       &        status='unknown')       &        status='unknown')
# Line 100  c file for ap_types assigned via coin fl Line 105  c file for ap_types assigned via coin fl
105           close(IniUnit1)           close(IniUnit1)
106           close(IniUnit2)           close(IniUnit2)
107  120      format(2i5)  120      format(2i5)
108            ENDIF
109    
110  c file of total absorption spectra  c file of total absorption spectra
111  c rows = pwaves, columns = np        c rows = pwaves, columns = np      
# Line 108  c rows = pwaves, columns = np Line 114  c rows = pwaves, columns = np
114               outfile(np,ilam) = aphy_chl(np,ilam)               outfile(np,ilam) = aphy_chl(np,ilam)
115             enddo             enddo
116            enddo            enddo
117            IF (MASTER_CPU_THREAD(myThid)) THEN
118           CALL MDSFINDUNIT( IniUnit3, mythid )           CALL MDSFINDUNIT( IniUnit3, mythid )
119           open(IniUnit3,file='p-ini-char-apspec.dat',           open(IniUnit3,file='p-ini-char-apspec.dat',
120       &        status='unknown')       &        status='unknown')
# Line 125  c rows = pwaves, columns = np Line 132  c rows = pwaves, columns = np
132  c make sure outfile is defined above with the correct size  c make sure outfile is defined above with the correct size
133           close(IniUnit3)           close(IniUnit3)
134           close(IniUnit4)           close(IniUnit4)
135            ENDIF
136    
137  c file for absorption spectra of PS's only  c file for absorption spectra of PS's only
138  c rows = pwaves, columns = np  c rows = pwaves, columns = np
# Line 133  c rows = pwaves, columns = np Line 141  c rows = pwaves, columns = np
141               outfile_ps(np,ilam) = aphy_chl_ps(np,ilam)               outfile_ps(np,ilam) = aphy_chl_ps(np,ilam)
142             enddo             enddo
143            enddo            enddo
144            IF (MASTER_CPU_THREAD(myThid)) THEN
145           CALL MDSFINDUNIT( IniUnit5, mythid )           CALL MDSFINDUNIT( IniUnit5, mythid )
146           open( IniUnit5,file='p-ini-char-apspec-psc.dat',           open( IniUnit5,file='p-ini-char-apspec-psc.dat',
147       &         status='unknown')       &         status='unknown')
# Line 150  c rows = pwaves, columns = np Line 159  c rows = pwaves, columns = np
159           close(IniUnit5)           close(IniUnit5)
160           close(IniUnit6)           close(IniUnit6)
161  130      format(9999f10.4)  130      format(9999f10.4)
162            ENDIF
163    
164    
165  c file for wavebands used  c file for wavebands used
166            IF (MASTER_CPU_THREAD(myThid)) THEN
167  c        open(23,file='pwaves-check.dat',status='new')  c        open(23,file='pwaves-check.dat',status='new')
168           CALL MDSFINDUNIT( IniUnit7, mythid )           CALL MDSFINDUNIT( IniUnit7, mythid )
169           open( IniUnit7,file='pwaves-check.dat',           open( IniUnit7,file='pwaves-check.dat',
# Line 160  c        open(23,file='pwaves-check.dat' Line 171  c        open(23,file='pwaves-check.dat'
171           write(IniUnit7,140)pwaves           write(IniUnit7,140)pwaves
172           close(IniUnit7)           close(IniUnit7)
173  140      format(i5)  140      format(i5)
174            ENDIF
175    
176    #endif /* WAVEBANDS */
 #endif  
177    
178           return           return
179           end           end
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22