| 1 | call WRITEDIAGFI(ngrid,"emis","Surface emissivity","w.m-1",2, |
|---|
| 2 | . emis) |
|---|
| 3 | call WRITEDIAGFI(ngrid,"tsurf","Surface temperature","K",2, |
|---|
| 4 | . tsurf) |
|---|
| 5 | call WRITEDIAGFI(ngrid,"ps","surface pressure","K",2,ps) |
|---|
| 6 | call WRITEDIAGFI(ngrid,"co2ice","co2 ice thickness","kg.m-2",2, |
|---|
| 7 | . co2ice) |
|---|
| 8 | call WRITEDIAGFI(ngrid,"temp","temperature","K",3,zt) |
|---|
| 9 | call WRITEDIAGFI(ngrid,"u","Zonal wind","m.s-1",3,zu) |
|---|
| 10 | call WRITEDIAGFI(ngrid,"v","Meridional wind","m.s-1",3,zv) |
|---|
| 11 | call WRITEDIAGFI(ngrid,"tsoil","soil temperature","K",3,tsoil) |
|---|
| 12 | if (tracer) then |
|---|
| 13 | do iq=1,nqmx |
|---|
| 14 | write(str2(1:2),'(i2.2)') iq |
|---|
| 15 | call WRITEDIAGFI(ngridmx,'q'//str2,noms(iq), |
|---|
| 16 | & 'kg/kg',3,zq(1,1,iq)) |
|---|
| 17 | call WRITEDIAGFI(ngridmx,'qsurf'//str2,noms(iq), |
|---|
| 18 | & 'kg.m-2',2,qsurf(1,iq)) |
|---|
| 19 | end do |
|---|
| 20 | end if |
|---|
| 21 | |
|---|