Changeset 75

Timestamp:

08/02/12 22:32:22 (12 years ago)

Author:

ymipsl

Message:

write file apbp.nc parameters for vertical discretisation

YM

Location:

codes/icosagcm/trunk/src

Files:

• disvert.f90 (modified) (1 diff)
• icosa_gcm.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/disvert.f90

@@ -47 +47 @@
   END SUBROUTINE init_disvert  
   
+  SUBROUTINE write_apbp
+  USE icosa
+  USE netcdf_mod
+  IMPLICIT NONE
+    REAL(rstd) :: val(llm)
+    INTEGER :: status
+    INTEGER :: lev,ilev
+    INTEGER :: ncid,levid,ilevid,hyaiid,hybiid,hyamid,hybmid,P0id
+    INTEGER :: l
+    
+    status = NF90_CREATE('apbp.nc', NF90_CLOBBER, ncid)
+    status = NF90_DEF_DIM(ncid,'lev',llm,lev)
+    status = NF90_DEF_DIM(ncid,'ilev',llm+1,ilev)
+    
+    status = NF90_DEF_VAR(ncid,'lev',NF90_DOUBLE,(/ lev /),levid)
+    status = NF90_PUT_ATT(ncid,levid,"long_name","hybrid level at midpoints (1000*(A+B))")
+    status = NF90_PUT_ATT(ncid,levid,"units","Pa")
+    status = NF90_PUT_ATT(ncid,levid,"positive","up")
+    status = NF90_PUT_ATT(ncid,levid,"standard_name","atmosphere_hybrid_sigma_pressure_coordinate")
+    status = NF90_PUT_ATT(ncid,levid,"formula_terms","a: hyam b: hybm p0: P0 ps: PS")
+
+    status = NF90_DEF_VAR(ncid,'ilev',NF90_DOUBLE,(/ ilev /),ilevid)
+    status = NF90_PUT_ATT(ncid,ilevid,"long_name","hybrid level at interfaces (1000*(A+B))")
+    status = NF90_PUT_ATT(ncid,ilevid,"units","Pa")
+    status = NF90_PUT_ATT(ncid,ilevid,"positive","up")
+    status = NF90_PUT_ATT(ncid,ilevid,"standard_name","atmosphere_hybrid_sigma_pressure_coordinate")
+    status = NF90_PUT_ATT(ncid,ilevid,"formula_terms","a: hyai b: hybi p0: P0 ps: PS")
+   
+    status = NF90_DEF_VAR(ncid,'hyai',NF90_DOUBLE,(/ ilev /),hyaiid)
+    status = NF90_PUT_ATT(ncid,hyaiid,"long_name","hybrid A coefficient at layer interfaces")
+
+    status = NF90_DEF_VAR(ncid,'hybi',NF90_DOUBLE,(/ ilev /),hybiid)
+    status = NF90_PUT_ATT(ncid,hybiid,"long_name","hybrid B coefficient at layer interfaces")
+
+    status = NF90_DEF_VAR(ncid,'hyam',NF90_DOUBLE,(/ lev /),hyamid)
+    status = NF90_PUT_ATT(ncid,hyamid,"long_name","hybrid A coefficient at midpoint interfaces")
+
+    status = NF90_DEF_VAR(ncid,'hybm',NF90_DOUBLE,(/ lev /),hybmid)
+    status = NF90_PUT_ATT(ncid,hybmid,"long_name","hybrid B coefficient at midpoint interfaces")
+    
+    status = NF90_DEF_VAR(ncid,'P0',NF90_DOUBLE,varid=P0id)
+
+    status = NF90_ENDDEF(ncid)   
+    
+    PRINT*,ap
+    PRINT*,bp
+    status=NF90_PUT_VAR(ncid,ilevid, ap(:)+bp(:)*Preff)
+    
+    DO l=1,llm
+      val(l)= 0.5*(ap(l+1)+ap(l)+Preff*(bp(l)+bp(l+1)))
+    ENDDO
+    
+    status=NF90_PUT_VAR(ncid,levid, val)
+
+    status=NF90_PUT_VAR(ncid,hyaiid, ap(:)/Preff)
+    status=NF90_PUT_VAR(ncid,hybiid, bp(:)/Preff)
+    
+    DO l=1,llm
+      val(l)= 0.5*(ap(l+1)+ap(l))
+    ENDDO
+    status=NF90_PUT_VAR(ncid,hyamid, val(:)/Preff)
+   
+     DO l=1,llm
+      val(l)= 0.5*(bp(l+1)+bp(l))
+    ENDDO
+    status=NF90_PUT_VAR(ncid,hybmid, val(:))
+ 
+    status=NF90_PUT_VAR(ncid,P0id, Preff)
+   
+   status=NF90_CLOSE(ncid)
+  END SUBROUTINE write_apbp
+
 END MODULE disvert_mod

codes/icosagcm/trunk/src/icosa_gcm.f90

@@ -77 +77 @@
   CALL WriteField("Ai",geom%Ai)
 !  CALL WriteField("sum_ne",sum_ne)
+  CALL write_apbp
   CALL timeloop