Changeset 1135 for LMDZ4/branches/LMDZ4-dev/libf/phylmd/read_pstoke0.F
- Timestamp:
- Mar 26, 2009, 1:30:23 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/libf/phylmd/read_pstoke0.F
r940 r1135 17 17 C****************************************************************************** 18 18 19 19 use netcdf 20 20 USE dimphy 21 21 IMPLICIT NONE … … 121 121 if (irec .eq. 0) then 122 122 123 ncidp=NCOPN('phystoke.nc',NCNOWRIT,rcode)124 125 varidps=NCVID(ncidp,'phis',rcode)123 rcode=nf90_open('phystoke.nc',nf90_nowrite,ncidp) 124 125 rcode = nf90_inq_varid(ncidp, 'phis', varidps) 126 126 print*,'ncidp,varidps',ncidp,varidps 127 127 128 varidpl=NCVID(ncidp,'sig_s',rcode)128 rcode = nf90_inq_varid(ncidp, 'sig_s', varidpl) 129 129 print*,'ncidp,varidpl',ncidp,varidpl 130 130 131 varidai=NCVID(ncidp,'aire',rcode)131 rcode = nf90_inq_varid(ncidp, 'aire', varidai) 132 132 print*,'ncidp,varidai',ncidp,varidai 133 133 134 varidt=NCVID(ncidp,'t',rcode)134 rcode = nf90_inq_varid(ncidp, 't', varidt) 135 135 print*,'ncidp,varidt',ncidp,varidt 136 136 137 varidmfu=NCVID(ncidp,'mfu',rcode)137 rcode = nf90_inq_varid(ncidp, 'mfu', varidmfu) 138 138 print*,'ncidp,varidmfu',ncidp,varidmfu 139 139 140 varidmfd=NCVID(ncidp,'mfd',rcode)140 rcode = nf90_inq_varid(ncidp, 'mfd', varidmfd) 141 141 print*,'ncidp,varidmfd',ncidp,varidmfd 142 142 143 varidenu=NCVID(ncidp,'en_u',rcode)143 rcode = nf90_inq_varid(ncidp, 'en_u', varidenu) 144 144 print*,'ncidp,varidenu',ncidp,varidenu 145 145 146 variddeu=NCVID(ncidp,'de_u',rcode)146 rcode = nf90_inq_varid(ncidp, 'de_u', variddeu) 147 147 print*,'ncidp,variddeu',ncidp,variddeu 148 148 149 varidend=NCVID(ncidp,'en_d',rcode)149 rcode = nf90_inq_varid(ncidp, 'en_d', varidend) 150 150 print*,'ncidp,varidend',ncidp,varidend 151 151 152 varidded=NCVID(ncidp,'de_d',rcode)152 rcode = nf90_inq_varid(ncidp, 'de_d', varidded) 153 153 print*,'ncidp,varidded',ncidp,varidded 154 154 155 varidch=NCVID(ncidp,'coefh',rcode)155 rcode = nf90_inq_varid(ncidp, 'coefh', varidch) 156 156 print*,'ncidp,varidch',ncidp,varidch 157 157 158 158 c Thermiques 159 varidfmth=NCVID(ncidp,'fm_th',rcode)159 rcode = nf90_inq_varid(ncidp, 'fm_th', varidfmth) 160 160 print*,'ncidp,varidfmth',ncidp,varidfmth 161 161 162 varidenth=NCVID(ncidp,'en_th',rcode)162 rcode = nf90_inq_varid(ncidp, 'en_th', varidenth) 163 163 print*,'ncidp,varidenth',ncidp,varidenth 164 164 165 varidfi=NCVID(ncidp,'frac_impa',rcode)165 rcode = nf90_inq_varid(ncidp, 'frac_impa', varidfi) 166 166 print*,'ncidp,varidfi',ncidp,varidfi 167 167 168 varidfn=NCVID(ncidp,'frac_nucl',rcode)168 rcode = nf90_inq_varid(ncidp, 'frac_nucl', varidfn) 169 169 print*,'ncidp,varidfn',ncidp,varidfn 170 170 171 varidyu1=NCVID(ncidp,'pyu1',rcode)171 rcode = nf90_inq_varid(ncidp, 'pyu1', varidyu1) 172 172 print*,'ncidp,varidyu1',ncidp,varidyu1 173 173 174 varidyv1=NCVID(ncidp,'pyv1',rcode)174 rcode = nf90_inq_varid(ncidp, 'pyv1', varidyv1) 175 175 print*,'ncidp,varidyv1',ncidp,varidyv1 176 176 177 varidfts1=NCVID(ncidp,'ftsol1',rcode)177 rcode = nf90_inq_varid(ncidp, 'ftsol1', varidfts1) 178 178 print*,'ncidp,varidfts1',ncidp,varidfts1 179 179 180 varidfts2=NCVID(ncidp,'ftsol2',rcode)180 rcode = nf90_inq_varid(ncidp, 'ftsol2', varidfts2) 181 181 print*,'ncidp,varidfts2',ncidp,varidfts2 182 182 183 varidfts3=NCVID(ncidp,'ftsol3',rcode)183 rcode = nf90_inq_varid(ncidp, 'ftsol3', varidfts3) 184 184 print*,'ncidp,varidfts3',ncidp,varidfts3 185 185 186 varidfts4=NCVID(ncidp,'ftsol4',rcode)186 rcode = nf90_inq_varid(ncidp, 'ftsol4', varidfts4) 187 187 print*,'ncidp,varidfts4',ncidp,varidfts4 188 188 189 varidpsr1=NCVID(ncidp,'psrf1',rcode)189 rcode = nf90_inq_varid(ncidp, 'psrf1', varidpsr1) 190 190 print*,'ncidp,varidpsr1',ncidp,varidpsr1 191 191 192 varidpsr2=NCVID(ncidp,'psrf2',rcode)192 rcode = nf90_inq_varid(ncidp, 'psrf2', varidpsr2) 193 193 print*,'ncidp,varidpsr2',ncidp,varidpsr2 194 194 195 varidpsr3=NCVID(ncidp,'psrf3',rcode)195 rcode = nf90_inq_varid(ncidp, 'psrf3', varidpsr3) 196 196 print*,'ncidp,varidpsr3',ncidp,varidpsr3 197 197 198 varidpsr4=NCVID(ncidp,'psrf4',rcode)198 rcode = nf90_inq_varid(ncidp, 'psrf4', varidpsr4) 199 199 print*,'ncidp,varidpsr4',ncidp,varidpsr4 200 200
Note: See TracChangeset
for help on using the changeset viewer.