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