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