Changeset 36 for branches/petra/src/writeheader_txt.f90
- Timestamp:
- Feb 16, 2015, 6:21:09 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/petra/src/writeheader_txt.f90
r20 r36 1 1 !********************************************************************** 2 ! Copyright 1998 ,1999,2000,2001,2002,2005,2007,2008,2009,2010*2 ! Copyright 1998-2015 * 3 3 ! Andreas Stohl, Petra Seibert, A. Frank, Gerhard Wotawa, * 4 4 ! Caroline Forster, Sabine Eckhardt, John Burkhart, Harald Sodemann * … … 32 32 ! * 33 33 ! 7 August 2002 * 34 ! modified IP 2013 for text output * 34 ! * 35 ! 2013, IP: IP, text output * 36 ! 2/2015, PS: version string length written, enclose version string in ' 37 ! write out ldep_incr 38 ! 35 39 !***************************************************************************** 36 40 ! * … … 59 63 !************************ 60 64 61 open(unitheader,file=path(2)(1:length(2))//'header _txt', &65 open(unitheader,file=path(2)(1:length(2))//'header.txt', & 62 66 form='formatted',err=998) 63 open(unitheader_ txt,file=path(2)(1:length(2))//'header_txt_releases', &67 open(unitheader_rel,file=path(2)(1:length(2))//'header_releases.txt', & 64 68 form='formatted',err=998) 65 69 … … 68 72 !***************************** 69 73 70 write(unitheader,*) '# ibdate,ibtime, iedate, ietime, flexversion'71 write(unitheader,*) ibdate, ibtime, iedate, ietime, trim(flexversion)! 'FLEXPART V9.0'74 write(unitheader,*) '# ibdate,ibtime, iedate, ietime, len(flexversion), flexversion' 75 write(unitheader,*) ibdate, ibtime, iedate, ietime, len_trim(flexversion), "'"//trim(flexversion)//"'" ! 'FLEXPART V9.0' 72 76 !if (ldirect.eq.1) then 73 77 ! write(unitheader,*) ibdate,ibtime,trim(flexversion) ! 'FLEXPART V9.0' … … 87 91 write(unitheader,*) '# information on grid setup ' 88 92 write(unitheader,*) '#outlon0,outlat0,numxgrid,numygrid,dxout,dyout' 89 write(unitheader,*) outlon0,outlat0,numxgrid,numygrid, & 90 dxout,dyout 93 write(unitheader,*) outlon0,outlat0,numxgrid,numygrid, dxout,dyout 91 94 write(unitheader,*) '# numzgrid, outheight(.) ' 92 95 write(unitheader,*) numzgrid,(outheight(i),i=1,numzgrid) … … 119 122 120 123 121 write(unitheader_ txt,*) '# information on release points'122 write(unitheader_ txt,*) '# numpoint'123 write(unitheader_ txt,*) numpoint124 write(unitheader_ txt,*) '# for numpoint:'124 write(unitheader_rel,*) '# information on release points' 125 write(unitheader_rel,*) '# numpoint' 126 write(unitheader_rel,*) numpoint 127 write(unitheader_rel,*) '# for numpoint:' 125 128 do i=1,numpoint 126 write(unitheader_ txt,*) ireleasestart(i),ireleaseend(i),kindz(i)129 write(unitheader_rel,*) ireleasestart(i),ireleaseend(i),kindz(i) 127 130 xp1=xpoint1(i)*dx+xlon0 128 131 yp1=ypoint1(i)*dy+ylat0 129 132 xp2=xpoint2(i)*dx+xlon0 130 133 yp2=ypoint2(i)*dy+ylat0 131 write(unitheader_ txt,*) xp1,yp1,xp2,yp2,zpoint1(i),zpoint2(i)132 write(unitheader_ txt,*) npart(i),1134 write(unitheader_rel,*) xp1,yp1,xp2,yp2,zpoint1(i),zpoint2(i) 135 write(unitheader_rel,*) npart(i),1 133 136 if (numpoint.le.1000) then 134 write(unitheader_ txt,*) compoint(i)137 write(unitheader_rel,*) compoint(i) 135 138 else 136 write(unitheader_ txt,*) compoint(1001)139 write(unitheader_rel,*) compoint(1001) 137 140 endif 138 141 do j=1,nspec 139 write(unitheader_ txt,*) xmass(i,j)140 write(unitheader_ txt,*) xmass(i,j)141 write(unitheader_ txt,*) xmass(i,j)142 write(unitheader_rel,*) xmass(i,j) 143 write(unitheader_rel,*) xmass(i,j) 144 write(unitheader_rel,*) xmass(i,j) 142 145 end do 143 146 end do … … 148 151 write(unitheader,*) '# information on model switches' 149 152 write(unitheader,*) '# method,lsubgrid,lconvection, ind_source,ind_receptor' 150 write(unitheader,*) method,lsubgrid,lconvection, & 151 ind_source,ind_receptor 153 write(unitheader,*) method,lsubgrid,lconvection,ind_source,ind_receptor 152 154 153 155 ! Write age class information … … 160 162 !Do not write topography to text output file. Keep it on the binary one 161 163 !******************************** 162 163 164 !do ix=0,numxgrid-1 164 165 ! write(unitheader,*) (oroout(ix,jy),jy=0,numygrid-1) 165 166 !end do 166 167 167 168 ! Write deposition type 169 !*********************** 168 170 169 171 write(unitheader,*) '# information on incremental / accum. deposition' 172 write(unitheader,*) ldep_incr 170 173 171 174 close(unitheader) 172 close(unitheader_txt) 173 175 close(unitheader_rel) 174 176 175 177 ! open(unitheader,file=path(2)(1:length(2))//'header_nml', & … … 181 183 182 184 183 998 write(*,*) ' #### FLEXPART MODEL ERROR! THE FILE #### ' 185 998 continue 186 write(*,*) ' #### FLEXPART MODEL ERROR! THE FILE #### ' 184 187 write(*,*) ' #### '//path(2)(1:length(2))//'header_txt'//' #### ' 185 188 write(*,*) ' #### CANNOT BE OPENED. IF A FILE WITH THIS #### '
Note: See TracChangeset
for help on using the changeset viewer.