Changeset 5100 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common
- Timestamp:
- Jul 23, 2024, 7:00:20 AM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/grilles_gcm_netcdf_sub.F90
r5099 r5100 122 122 !------------------------- 123 123 ! Longitudes on "u" dynamical grid 124 status= NF90_DEF_VAR(ncid_out,'lonu',NF90_FLOAT,out_lonudim, out_lonuid)124 status=nf90_def_var(ncid_out,'lonu',nf90_float,out_lonudim, out_lonuid) 125 125 CALL handle_err(status) 126 126 status=nf90_put_att(ncid_out,out_lonuid,'units','degrees_east') 127 127 status=nf90_put_att(ncid_out,out_lonuid,'long_name','Longitude on u grid') 128 128 ! Longitudes on "v" dynamical grid 129 status= NF90_DEF_VAR(ncid_out,'lonv',NF90_FLOAT,out_lonvdim, out_lonvid)129 status=nf90_def_var(ncid_out,'lonv',nf90_float,out_lonvdim, out_lonvid) 130 130 CALL handle_err(status) 131 131 status=nf90_put_att(ncid_out,out_lonvid,'units','degrees_east') 132 132 status=nf90_put_att(ncid_out,out_lonvid,'long_name','Longitude on v grid') 133 133 ! Latitudes on "u" dynamical grid 134 status= NF90_DEF_VAR(ncid_out,'latu',NF90_FLOAT,out_latudim, out_latuid)134 status=nf90_def_var(ncid_out,'latu',nf90_float,out_latudim, out_latuid) 135 135 CALL handle_err(status) 136 136 status=nf90_put_att(ncid_out,out_latuid,'units','degrees_north') 137 137 status=nf90_put_att(ncid_out,out_latuid,'long_name','Latitude on u grid') 138 138 ! Latitudes on "v" dynamical grid 139 status= NF90_DEF_VAR(ncid_out,'latv',NF90_FLOAT,out_latvdim, out_latvid)139 status=nf90_def_var(ncid_out,'latv',nf90_float,out_latvdim, out_latvid) 140 140 CALL handle_err(status) 141 141 status=nf90_put_att(ncid_out,out_latvid,'units','degrees_north') … … 144 144 out_dim(1)=out_lonudim 145 145 out_dim(2)=out_latudim 146 status= NF90_DEF_VAR(ncid_out,'grille_u',NF90_FLOAT,out_dim, out_uid)146 status=nf90_def_var(ncid_out,'grille_u',nf90_float,out_dim, out_uid) 147 147 CALL handle_err(status) 148 148 status=nf90_put_att(ncid_out,out_uid,'units','m/s') … … 151 151 out_dim(1)=out_lonvdim 152 152 out_dim(2)=out_latvdim 153 status= NF90_DEF_VAR(ncid_out,'grille_v',NF90_FLOAT,out_dim, out_vid)153 status=nf90_def_var(ncid_out,'grille_v',nf90_float,out_dim, out_vid) 154 154 CALL handle_err(status) 155 155 status=nf90_put_att(ncid_out,out_vid,'units','m/s') … … 158 158 out_dim(1)=out_lonvdim 159 159 out_dim(2)=out_latudim 160 status= NF90_DEF_VAR(ncid_out,'grille_s',NF90_FLOAT,out_dim, out_tempid)160 status=nf90_def_var(ncid_out,'grille_s',nf90_float,out_dim, out_tempid) 161 161 CALL handle_err(status) 162 162 status=nf90_put_att(ncid_out,out_tempid,'units','Kelvin') … … 165 165 ! for INCA : 166 166 ! vertical levels "presnivs" 167 status= NF90_DEF_VAR(ncid_out,'presnivs',NF90_FLOAT,out_levdim, presnivs_id)167 status=nf90_def_var(ncid_out,'presnivs',nf90_float,out_levdim, presnivs_id) 168 168 CALL handle_err(status) 169 169 status=nf90_put_att(ncid_out,presnivs_id,'units','Pa') … … 172 172 out_dim(1)=out_lonvdim 173 173 out_dim(2)=out_latudim 174 status = nf90_def_var(ncid_out,'phis', NF90_FLOAT,out_dim,phis_id)174 status = nf90_def_var(ncid_out,'phis',nf90_float,out_dim,phis_id) 175 175 CALL handle_err(status) 176 176 status=nf90_put_att(ncid_out,phis_id,'units','m') 177 177 status=nf90_put_att(ncid_out,phis_id,'long_name','surface geopotential height') 178 178 ! gridcell area 179 status = nf90_def_var(ncid_out,'aire', NF90_FLOAT,out_dim,area_id)179 status = nf90_def_var(ncid_out,'aire',nf90_float,out_dim,area_id) 180 180 CALL handle_err(status) 181 181 status=nf90_put_att(ncid_out,area_id,'units','m2') … … 193 193 !------------------------- 194 194 ! 1D : lonu, lonv,latu,latv ; INCA : presnivs 195 status= NF90_PUT_VAR(ncid_out,out_lonuid,rlonudeg,[1],[iip1])196 CALL handle_err(status) 197 status= NF90_PUT_VAR(ncid_out,out_lonvid,rlonvdeg,[1],[iip1])198 CALL handle_err(status) 199 status= NF90_PUT_VAR(ncid_out,out_latuid,rlatudeg,[1],[jjp1])200 CALL handle_err(status) 201 status= NF90_PUT_VAR(ncid_out,out_latvid,rlatvdeg,[1],[jjm])202 CALL handle_err(status) 203 status= NF90_PUT_VAR(ncid_out,presnivs_id,rlev,[1],[llm])195 status=nf90_put_var(ncid_out,out_lonuid,rlonudeg,[1],[iip1]) 196 CALL handle_err(status) 197 status=nf90_put_var(ncid_out,out_lonvid,rlonvdeg,[1],[iip1]) 198 CALL handle_err(status) 199 status=nf90_put_var(ncid_out,out_latuid,rlatudeg,[1],[jjp1]) 200 CALL handle_err(status) 201 status=nf90_put_var(ncid_out,out_latvid,rlatvdeg,[1],[jjm]) 202 CALL handle_err(status) 203 status=nf90_put_var(ncid_out,presnivs_id,rlev,[1],[llm]) 204 204 CALL handle_err(status) 205 205 … … 209 209 210 210 COUNT(2)=jjp1 ! for "u" and "s" grids 211 status= NF90_PUT_VAR(ncid_out,out_uid,uwnd,start, count)211 status=nf90_put_var(ncid_out,out_uid,uwnd,start, count) 212 212 CALL handle_err(status) 213 213 COUNT(2)=jjm ! for "v" grid 214 status= NF90_PUT_VAR(ncid_out,out_vid,vwnd,start, count)214 status=nf90_put_var(ncid_out,out_vid,vwnd,start, count) 215 215 CALL handle_err(status) 216 216 COUNT(2)=jjp1 ! as "s" grid, for all the following vars 217 status= NF90_PUT_VAR(ncid_out,out_tempid,temp,start, count)217 status=nf90_put_var(ncid_out,out_tempid,temp,start, count) 218 218 CALL handle_err(status) 219 219 status = nf90_put_var(ncid_out, phis_id, phis_loc,start,count)
Note: See TracChangeset
for help on using the changeset viewer.